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ABSTRACT 


We  have  developed  a simple  but  general  model  in  terms  of  which  the  electronic 
structure  and  spectra  of  Si02.  Ge02  and  the  ABO^-type  oxides  are  studied  in  a systematic 
way.  Methods  have  also  been  developed  in  terms  of  which  we  calculate  the  energy  levels  of 
impurities  in  bulk  Si02  and  at  the  Si>Si02  interface,  and  the  electronic  structure  of  free 
surfaces  and  interfaces  between  crystalline  materials.  The  main  objectives  of  the  work 
reported  here  were  to  obtain  a theoretical  description  of  the  electronic  structure  and  properties 
of  the  oxides  of  the  tetrahedral  semiconductors  (in  particular,  Si02.  Ge02,  and  the  ABO^-type 
oxides,  where  AB  is  a tetrahedral  semiconductor)  and  their  interfaces  with  other  materials  and 
vacuum.  The  results  of  the  work  are  summarized  below.  Details  are  contained  in  Appendices 
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1.  BULK  ELECTRONIC  PROPERTIES 

We  have  developed  a simple  tight>binding  model,  based  on  ideas  previously  developed 
and  tested  for  the  simple  tetrahedral  semiconductors,*'^  in  terms  of  which  three  distinct  kinds 
of  properties  were  studied:  (a)  Average  Properties:  These  are  properties  which  correspond  to 
an  average  over  all  the  one-electron  states.  Examples  are  effective  charges  on  atoms,  dielectric 
constants,  cohesive  energies,  etc.  (b)  Dispersive  Properties:  These  are  properties  which 
depend  on  the  details  of  the  distribution  of  one-electron  states.  Examples  are  energy  bands, 
densities  of  states,  photoemission  spectra.  X-ray  emission  spectra,  etc.  (c)  Local  Properties, 
such  as  impurity  and  defect  levels  and  excitons. 

The  main  achievements  of  our  work  using  this  model  are  the  following:  The  first  sets 
of  energy  bands  for  Si02,  Ge02  and  the  ABO^-type  oxides  were  obtained  (Ref.  3 and 
Appendix  A).  The  photoemission  spectrum  and  the  X-ray  emission  spectra  of  Si02  were 
calculated  explicitly  for  the  first  time  and  compared  with  experiment,  as  opposed  to  previous 
work  which  arrived  at  interpretations  of  the  spectra  in  terms  of  electronic  transitions  (Appendix 
A).  The  photoemission  spectra  of  Si^Gej,^02  alloys  were  also  calculated  and  found  to  agree 
well  with  experimental  dau  (Appendix  B).  Two  different  complementary  model  calculations 
of  the  optical  absorption  spectrum  of  Si02  were  carried  out  which  established  that  the 
observed  peaks  are  excitonic  in  nature  (Appendix  C).  Finally,  by  calculating  the  dielectric 
constant  of  Si02  explicitly  as  a function  of  the  Si-O-Si  angle,  it  was  possible  to  deduce 
information  about  the  structure  of  various  pol3rmorphs  (Appendices  A and  D). 
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2.  IMPURITIES  IN  BULK  SiO, 

The  above  model  for  the  bulk  properties  of  Si02  has  been  extended  to  provide  a 
description  for  substitutional  impurities  in  the  bulk  material.  In  particular,  we  have  studied  P 
and  As  at  an  O site  and  predicted  their  energy  levels  in  the  band  gap.  By  using  these  results 
we  were  able  to  interpret  experimental  dau  on  ion-implanud  Si02  and  identify  the  observed 
centers  as  substitutional  P and  Aa  at  O sites  (Appendix  E). 
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3.  THE  STOICHIOMETRY  OF  THE  Si-SiO,  INTERFACE 

A large  number  of  experiments  have  been  carried  out  which  have  attempted  to 
determine  whether  there  exists  an  SiO^  (0  < x < 2)  layer  at  the  Si-Si02  interface,  and,  if  so, 
determine  its  width.  In  our  study  (Appendix  F),  we  pointed  out  that  the  definition  of  a width 
is  not  unique  and  suggested  several  possibilities  that  would  give  different  numbers.  We  also 
constructed  continuous-random-network  (CRN)  models  for  the  (100)  interface  and  obtained 
the  following:  (a)  It  is  possible  to  construct  models  with  a totally  abrupt  interface  which  have 
distortions  that  are  visually  comparable  to  distortions  present  in  models  with  a finite  SiO^ 
layer,  (b)  A simple  force  model  was  used  to  perform  computer  calculations  which  showed  that 
the  net  strain  energy  was  reduced  as  the  finite  SiO^  layer  was  gradually  reduced  in  width, 
suggesting  that  an  abrupt  interface  may  in  fact  be  preferred.  (See  Appendix  F for  details.) 

4.  SURFACES  AND  INTERFACES 

In  this  work  the  main  goal  was  to  develop  efficient  methods  to  calculate  the  electronic 
structure  of  surfaces  and  interfaces.  At  this  stage,  the  work  was  restricted  to  crystalline 
materials  for  which  two-dimensional  periodicity  could  be  exploited.  We  made  use  of  the 
Koster-Slater**  idea  for  localized  perturbations  which  was  first  introduced  for  point  defects^ 
and  was  later  extended  to  surfaces  by  Koutecky^  and  others.  The  advantage  of  the  method  is 
that  it  builds  in  the  properties  of  the  perfect  bulk  crystal  from  the  start  and  then  directly 
calculates  changes  in  the  electronic  structure  produced  by  the  perturbation.  In  contrast, 
alternative  methods,  such  as  cluster  and  slab  methods,^  rely  on  a finite  number  of  atoms  or 
layers  of  atoms  to  give  an  adequate  description  of  both  bulk  properties  (e.g.,  bandwidths, 
bandgaps)  and  localized  states.  In  the  case  of  surfaces,  we  made  use  of  empirical  tight-binding 
Hamiltonians  and  introduced  a novel  way  to  define  the  perturbation  that  creates  free  surfaces. 
Instead  of  "cutting  bonds"  to  divide  an  infinite  solid  into  two  decoupled  semi-infinite  solids, 
we  "remove"  layers  of  atoms  by  setting  their  on-site  Hamiltonian  matrix  elemenu  to  infinity,  a 
procedure  first  used  to  describe  a single  vacancy.^.  This  choice  of  pertnrbation  makes  the 


calculations  straightforward  and  fast.  Calculations  have  been  carried  out  for  test  purposes  and 
for  cubic  SiOj.  The  details  of  the  method  and  the  results  are  discussed  in  Appendix  G. 

We  also  extended  the  Koster-Slater  approach  to  describe  interfaces  between  crystalline 
materials  by  assuming  the  unperturbed  system  to  consist  of  two  noninteracting  infinite  solids 
and  then  constructing  an  appropriate  perturbation  matrix.  Calculations  so  far  have  only  been 
done  on  the  Ge-GaAs  and  Ge*ZnSe  systems  for  which  other  theoretical  work  is  also  available. 
The  method  and  results  are  discussed  in  Appendix  H. 

The  perturbation  matrices  for  free  surfaces  and  for  interfaces  may  be  combined  so  that 
calculations  can  be  performed  for  monolayers  or  multiple  layers  of  a given  material  on  a 
crystalline  substrate.  This  technique  is  therefore  capable  of  providing  information  about 
gradual  oxidation  of  materials. 


5.  IMPURITIES  AT  INTERFACES 

Theoretical  work  has  been  done  in  two  different  areas; 

(a)  It  has  been  observed^  that  some  of  the  traps  introduced  into  Si02  near  a Si 

interface  by  ionizing  radiation  and  by  electron-beam  irradiation  (procedures  used  in  the 

I 

''  fabrication  of  some  devices)  do  not  anneal  out.  There  are  two  types  of  such  traps,  one  with  a 

I 

j cross  section  of  about  10*'^  cm^,  and  the  other  with  a cross  section  of  about  10'**  cm^.  The 

latter  was  not  produced  in  samples  irradiated  at  intensities  below  200  W/cm^  at  any  tempera- 
I ture.  The  rate  of  production  of  the  former  was  suppressed  by  lowering  the  irradiation 

I temperature  down  to  90K.  These  data  have  been  explained  in  terms  of  two  different  process- 

es, namely  the  formation  of  point-defect  bound  complexes  by  recombination-enhanced 

I diffusion  and  by  dislocation  migration.  Details  are  given  in  Appendix  I. 

1 

1 • 

(b)  We  have  developed  a method  to  calculate  binding  energies  of  shallow  impurities  at 
or  near  the  Si-SiOj  interface.  The  method  is  based  on  effective-mass  theory,  which  is  adapted 
for  the  cylindricai  symmetry  of  the  interface  problem.  Results  and  comparison  with  experi- 
mental data  are  contained  in  Appendix  J. 
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Some  properties  of  the  oxides  of  the  tetrahedral 
semiconductors  and  the  oxide-semiconductor  interfaces*^ 

Sokrates  T.  Pantelides 

IBM  Thomas  J.  IValson  Rtsianh  Ctnttr,  Yorklown  Htights,  Sew  York  10598 
(Received  10  Febniary  1977;  accepCed  4 April  1977) 

Continuous>nndoni-network  models  have  been  constructed  for  the  Si-SiOj  interface.  It  is 
found  that  an  abrupt  interface  with  no  SiO,  layer  is  posible.  A simple  tight-binding  model  is 
described  that  is  applicable  fo'  the  calculation  of  the  electronic  properties  of  the  bulk  oxides 
and  the  oxide-semiconductor  interfaces.  Results  are  given  only  for  selected  bulk  properties, 
namely  the  photoemission  and  x-ray  emission  spectra,  and  the  dielectric  constants  of  SiO„ 
OeO],  and  various  AB04-type  oxides. 

PACS  numbers:  71.43.Nt,  79.60.Eq,  73.40.(}v,  71.20. -l-c 


I.  INTRODUCTION 

In  order  to  develop  a theory  for  the  electronic  properties 
of  a given  nuterial,  it  is  essential  that  the  stoichiometry  and 
the  atomic  surangement  be  known.  In  this  paper,  a step  is 
taken  toward  understanding  the  stoichiometry,  atomic  ar- 
rangement, and  electronic  properties  of  the  lutive  oxides  of 
the  tetrahedral  semiconductors  and  the  oxide-semiconductor 
interfaces. 

II.  STOICHIOMETRY  AND  ATOMIC 
ARRANGEMENT  OF  INTERFACES 

The  native  oxide  of  silicon,  SiOz,  exists  in  a variety  of 
crystalline  forms  (quartz,  tridymite,  etc.i),  but  the  filnu  grown 
themudly  on  Si  for  the  manufacture  of  devices  are  anwrphous. 
As  for  the  Si-SiOz  interface,  the  situation  is  complicated  even 
moreby  the  question  whether  an  SiOx  layer,  withz  ^ 2,  exists 
between  Si  and  SiOz-  Answers  to  this  question  have  been 
sought  by  various  experimental  techniques,^  but  evidence 
thus  far  is  inconclusive.  Estimates  of  the  width  of  the  SiO* 
layer  include  values  of  essentially  zero^  (abrupt  interface), 
about  4 A (Ref.  3),  12-15  A (Ref.  4),  and  even  more.’-^  In  an 
attempt  to  obtain  some  insight  into  this  question,  a number 
of  continuous-random-network  (CRN)  models  have  been 
built  Detaib  of  the  constructian  will  be  given  elsewhere.^  The 
main  idea,  however,  was  to  start  with  a crystalline  Si  substrate, 
expose  a well-defined  surface,  and  then  attempt  to  build  up 


Fic.  1.  Pbatograph  of  the  CRN  model  deecribed  in  the  text. 

N8  J.  Vae.  Set  TeehnoL,  Vet  14,  No.  4,  July/Aug.  1977 


an  amorphous  oxide  layer  bond  by  bond.^  Clearly,  an  SiOi- 
type  interface  layer  of  arbitrary  width  can  be  built  before  full 
SiOz  coordination  is  reached,  and  therefore  the  aim  was  to 
strive  for  a minimum  interface  width.  The  first  modeb  that 
were  built  indicated  that  a finite  minimum  width  existed,  but 
finally  a model  was  constructed  with  zero  interface  width 
(abrupt  interface).  No  dangling  bonds  were  allowed  and  the 
local  distortions  are  comparable  to  those  present  in  bulk 
amorphous  CRN  modeb.  Thb  model  b shown  in  Fig.  1. 

The  construction  of  the  model  indicates  that  an  atomically 
sharp  interface  may  in  fact  be  energetically  possible,  but  does 
not  rule  out  the  possibility  that  some  real  interfaces  have  an 
SiOz  layer  of  finite  width,  which  may  be  a function  of  the 
growth  conditions.  An  attempt  b now  under  way  to  quantify 
the  amount  of  local  distortions  in  the  various  modeb  in  order 
to  extract  preferred  configurations.  Work  b also  in  progress 
to  calculate  the  electronic  properties  of  the  Si-SiOz  interface, 
assuming  the  atomic  anangements  of  these  models.  The  na- 


Flc.2.  'nwbaeioihilih  lot  the  mcdel  dmeribed  In  the  text  («)  orfattilnaed 
in  the  ideal  cubic  (d-crMobalite)  stmchire;  (b)  orbitak  uad  (or  the  mote 
fenenl  caee.  The  Si-O-Si  angle  b 144*. 
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ture  of  interface  states  may  suggest  preferred  configurations. 
The  basic  aspects  of  the  model  used  in  this  calculations  and 
results  on  bulk  properties  are  described  in  the  next  section. 

It  should  be  noted  that  the  CRN  models  described  above 
apply  to  the  compound  semiconductors  as  well,  assuming  that 
the  native  oxide  is  of  the  ABO4  type.  Other  possibilities,  of 
course,  exist.  For  example,  any  of  three  oxides  may  be  formed 
on  CaAs,  namely  CagOs,  AlgOs,  and  GaAsOa. 

III.  BULK  ELECTRONIC  PROPERTIES 

SiOg  and  GeOg  are  special  cases  of  the  general  AB04-type 
oxide,  where  AB  is  a tetrahedral  semiconductor,  with  A « B 
* Si  and  A * B ■ Ge,  respectively.  The  properties  of  all  these 
materials  can  therefore  be  described  in  a systematic  manner. 
A simple  parametrized  tight-binding  model  has  previously^ 
been  introduced  for  SiOg  and  GeOg  and  has  now  been  ex- 
tended to  treat  the  AB04-type  oxides.  The  model  is  similar 
to  that  used  for  the  tetrahedral  semiconductors.  In  this 
paper,  the  model  will  be  discussed  briefly  and  a few  new 
applications  will  be  givea 

Central  to  the  understanding  of  the  model  is  the  concept 
of  a bonding  unit.  In  the  case  of  the  tetrahedral  semicon- 
ductors, one  constructs  the  familiar  sp^  hybrids  on  each  atom. 
A bonding  unit  is  then  defined  to  consist  of  two  hybrids,  h^ 
and  hg.  on  nearest  neighbors,  directed  toward  each  other.  One 
then  constructs  bonding  and  antibonding  orbitals  in  each 
bonding  unit.  By  focusing  on  the  bonding  orbitals  alone 
(bond-orbital  approximation^-**’),  two  distinct  kinds  of  cal- 
culations can  be  carried  out:  (a)  Calculations  of  the  valence 
energy  bands  by  setting  up  an  LCAO  secular  matrix.  The 
form  of  the  bands  is  determined  by  the  off-diagonal  matrix 
elements,  Le.,  interactions  between  bonding  orbitals  on  dif- 
ferent sites,  (b)  Calculations  of  the  total  electronic  energy  both 
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Fig.  S.  The  UPS  spectrum  of  SiOj  (dotted  curve)  compered  with  the 
broadened  DOS  (solid  curve).  The  experimental  curve  is  from  B.  Fischer. 
R.  A.  Poliak,  T.  H.  DiStefano,and  W.  D.  Crobman,  Phys.  Rev.  B 15, 3193 

(1977). 

in  the  absence  and  presence  of  external  perturbations  and  thus 
of  explicit  analytical  expressions  for  response  functions.  For 
this  kind  of  calculation,  one  needs  only  the  diagonal  matrix 
elements  and  no  matrix  diagonalization  or  Brillouin  zone  sums 
are  necessary. 

The  bonding  unit  in  an  AB04-type  oxide  is  an  sp®  hybrid 
on  an  A atom,  and  sp®  hybrid  on  a B atom  and  oxygen  orbitals 
on  the  oxygen  atom  in  between.  For  valence-band  and 
total-energy  studies  the  oxygen  2p  orbitals  are  adequate.  The 
important  point  is  that  an  A-O  or  B-O  bond  is  not  a well- 
defined  entity  and  one  must  construct  A-O-B  bonding  or- 
bitals, which  will  be  referred  to  as  extended  bonding  orbitals 
(EBOs). 


Fig.  6.  Thmetkal  (lolhl  curva)  and  mpcrUnental  (dotted  curvet)  x-ray 
emiaian  ipectra  of  SiOt  The  experimental  curvet  are  from  C.  Wiech,  in 
Soft-X-Bay  Band  Spectra,  edited  by  D. ).  Fabian  (Academic,  New  York, 
1968),  and  C.  Klein,  and  H.  U.  Chun,  Phyt.  Stat.  Solidi  B 49, 167  (1972). 
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TablE'1.  Theoretical  and  experimental  values  for  the  electronic 
(high-frequency)  dielectric  constants  of  SiO,,  GeO,,  and  various 
ABO,-type  oxides. 


Material 

A-O-B  angle  (°) 

e(th.) 

«(exp») 

SiO,  (ot-quaitz) 

144b 

2.4 

2.4 

SiO,  (d^uaitz) 

144b 

2.3 

2.3 

SiO,  (oi-cristobalite) 

144b 

2.3 

2.2 

SiO,  (i3-cristobalite) 

180b 

2.7 

2.2 

SiO,  (fi-cristobalite) 

144C 

2.2 

2.2 

SiO,  (coesite) 

120b 

2.0 

2.6 

SiO,  (coesite) 

144c 

2.6 

2.6 

SiO,  (vitreous) 

144d 

2.1 

2.1 

GeO,  (quartzlike) 

130b 

2.9 

2.9 

BPO. 

135b 

2.0 

2.6 

BAsO. 

125b 

2.3 

2.8 

AIPC. 

145b 

2.4 

2.3 

AlAsO, 

146b 

2.9 

GaPO. 

135b 

2.6 

^ Landolt-Botnsteitt,  ZMenwerte  und  Funktionen  (Springer,  Berlin, 
1962).  Vol.  2,  pt.  8. 
b Reference  1. 
cSee  text, 
d .Assumed. 


The  internal  structure  of  the  EBO’s  is  rather  intriguing.  Let 

and  hg  denote  the  hybrids  in  a bonding  unit  and  px,  Py, 
Pt,  the  p orbitals  on  the  oxygen  (Fig.  2).  Instead  of  and  hg, 
are  given  elsewhere. In  this  paper  a few  new  results  will 
be  given  to  demonstrate  the  power  of  the  method.  Figure  3 
shows  the  valence  energy  bands  and  densities  of  states  (DOS) 
one  can  work  with  their  linear  combinations  + fts  and 
a = Ha  ~ ha-  Note  that  these  would  be  the  bonding  and  an- 
dbonding  orbitals,  respectively  in  an  elemental  semiconductor 
like  Si.  In  Si02.  the  main  EBO  is  formed  from  pt  and  the 
antibonding  combination,  a,  so  that  it  is  odd  about  the  mid- 
point between  the  two  Si’s.  The  structure  of  some  of  the  va- 
lence-band orbitals  of  Si02  is  thus  similar  to  the  structure  of 
the  conduction-band  orbitals  of  Si  and  vice-versa. 

In  constructing  the  EBO’s  for  SiOg,  two  quantities,  dpi  and 
dpi,  enter,  which  describe  the  charge  transfer  between  the 
oxygen  and  its  two  neighboring  Si’s,  in  constructing  the  EBO’s 
for  a typical  AfiOrtype  oxide,  b and  a defined  above  are  not 
adequate,  just  as  they  are  not  adequate  for  the  description  of 
zinc-blende-type  compounds. For  the  latter,  modified  h’s 
and  a’s  are  introduced,  given  by  h ■ (1  + ap)*^®hA  + (I  ~ 
«p)‘/*hBanda  * (1  - ap)*/*hA  “ (I  + «p)*^*bB,  where  Op, 
the  polarity,  describes  the  charge  transfer  between  atom  A 
and  B.  Similar  quantities  enter  the  construction  of  the  EBO’s 
for  the  ABOi-type  oxides.  They  are  denoted  by  apx  and  apt 
and  are  referred  to  as  external  polarities.  The  quantities  dpi 
and  dpi  are  referred  to  as  internal  polarities. 

Mathematical  detail  and  the  method  of  parametrization 


for  the  ideal  cubic  form  of  Si02  and  a typical  ABOi-type 
oxide.  Note  the  opening  of  an  additional  gap  in  ABO4,  which 
is  analogous  to  the  so-called  antisymmetric  gap  in  the  bands 
of  the  AB-type  semiconductors.*®  Figure  4 shows  the  DOS  for 
a more  general  Si02  structure,  which  is  a simulation®  of 
amorphous  Si02  and  has  the  Si-O-S^  angle  to  be  the  observed 
144°  instead  of  the  ideal  180°.  Figure  5 shows  the  same  DOS 
broadened  by  0.8  eV  and  compares  it  with  the  measured  ul- 
traviolet photoemission  sp^trum  (UPS).  The  agreement  be- 
tween theory  and  experiment  is  substantially  better  than  that 
obtained  previously**  with  only  0.5  eV  of  broadening.  Figure 
6 shows  the  calculated  and  experimental  x-ray-emission 
spectra  of  Si02.  This  is  the  first  such  calculation  of  these 
spectra  and  the  agreement  with  experiment  is  seen  to  be  very 
satisfactory.  Comparable  spectra  have  been  calculated  for 
ABOi-type  oxides  but  no  experiments  are  available.  One  of 
the  response  functions  that  can  be  calculated  without  needing 
the  energy  bands,  as  mentioned  earlier,  is  the  electronic 
(high-frequency)  dielectric  constant  <,  for  which  an  explicit 
analytical  expression  is  obtained.®  *^  It  turns  out  that,  in  ad- 
dition to  the  density  of  electrons,  t depends  on  the  angle  at  the 
oxygen  atom  in  an  important  way.  In  Table  I,  a number  of 
theoretical  and  experimental  values  of  < are  given.  In  the  case 
of  the  d-cristobalite  and  coesite  forms  of  Si02,  t was  calculated 
with  the  Si-O-Si  angles  listed  by  Wyckoff  * and  also  with  the 
quartz  value  of  144°.  The  180°  value  for  d-cristobalite  has 
been  questioned  by  many  authors.  The  comparison  seems  to 
suggest  that  the  two  extreme  angles  may  in  fact  be  unrealistic 
and  that  all  forms  of  Si02  have  angles  in  the  neighborhood  of 
144°*.  The  values  of  e listed  for  the  boron  compounds  are  not 
very  reliable  due  to  complications  with  first-row  elements.*® 
A more  extensive  treatment  of  this  and  other  questions  is  given 
in  Ref.  12. 

••Work  supported  in  pert  by  the  Office  of  Naval  Research  under  contract  Na 
N00014-76-C-0934. 
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THE  ELECTRONIC  STRUCTURE  OF  SiO,,  GeO,  AND  INTERMEDIATE  ! 

Si.Ge,..0,  COMPOSITIONS:  EXPERIMENT  AND  THEORY 

Sokrates  T.  Pantelkles,  Bernhard  Fischer.t  Roger  A.  PoUak,  and  Thomas  H.  Di  Stefano 
IBM  Thomas  J.  Watson  Research  Center,  Yorktown  Heights,  New  York  10598,  USA 
(Received  27  December  1976  by  E.  Burstein) 

Theoretical  calculations  are  reported  for  SiOj,  GeO,  and  the  intermediate  compos-  j 

itkx  SijGejO,  which  reproduce  the  main  observed  features  and  trends  in  experimen- 
tal photoemissioii  spectra.  The  agreement  between  the  two  establishes  the  impor-  < 

tance  of  band  theory  in  understanding  the  electronic  structure  of  these  materials,  and  { 

demonstrates  that  detailed  quantitative  predictions  are  feasible  for  such  complex 

materials  in  terms  of  the  empirical  tight-binding  method.  The  calculations  further  I 

establish  that  the  structure  in  the  valence  bands  is  determined  mainly  by  nearest-  ! 

neighbor  oxygen-oxygen  interactions. 

I 

i 

I 

The  electronic  structure  of  SiO,  and  GeO,  has  been  and  much  of  the  internal  structure  are  largely  deter- 

the  subject  of  considerable  attention  in  the  literature.  mined  by  the  nearest-neighbor  oxygen-oxygen  p-orbital 

Most  attempts  have,  however,  mainly  focused  on  interactions,  a factor  which  should  be  central  in  under- 

MwrprafiMg  experimental  data  (of^cal  absorption,  x-ray  standing  the  electronic  propeties  of  all  the  silicate  glass- 

absorption,  x-ray  emission,  photoemission.  etc.)  in  es. 

terms  of  simple  "molecular-orbital"  energy  levels.  Such  The  ultraviolet  photoemission  spectroscopy  (UPS) 

assignments  have  worked  well  for  certain  experimental  measurements  were  performed  using  a cylindiical  mir- 

spectra,  but  at  the  same  time  were  incompatible  with  ror  electrosutic-deflection  type  electron  energy  analyz- 

other  spectra.''^  Furthermore,  several  assignment  er  (resolution  0.2S  eV)  and  40.8  eV  photons  from  a 

schemes,  based  on  similar  ideas,  differ  substantially  He-discharge  lamp.  The  x-ray  photoemission  spectros- 

from  each  other.  One  difficulty  arises  from  the  fact  copy  (XPS)  measuremenu  were  performed  with  a 

that  the  crystal  structures  of  these  materials  are  com-  Hewlett-Packard  x-ray  photoelectron  spectrometer 

plex,  with  low  symmetry  and  many  atoms  in  a primitive  which  has  a monochromatic  A1  Ka, , x-ray  source 

unit  cell.  This  complexity  has  inUbited  the  successful  (1486.6  eV)  and  a resolution  of  0.6  eV!  First,  Si,Ge,., 

application  of  band-theoretic  techniques  which  have  in  films  several  thousand  Angstroms  thick  were  deposit^ 

the  past  proven  very  powerful  for  simpler  materials.  on  sapphire  substrates  by  DC  sputtering  from  targets 

Recently,  Pantelkles  and  Harrison^  (PH)  and  Schneider  formed  by  melting  high-purity  polycrystalline  silicon 

and  Fowler^  (SF)  performed  energy  band  calculations  and  germanium  onto  a molybdenun  substrate.  The 

for  SiO,  and  independently  pointed  out  the  importance  Si,Ge,.,0,  films  which  were  studied  with  UPS  were 

of  energy-band  theory  for  a systematic  interpretation  of  oxidized  in  jim  in  the  spectrometer  by  heating  to  4S0" 

the  data.  PH  went  beyond  SF  in  calculating  densities  of  C for  IS  minutes  in  10"’  torr  of  oxygen.  The  Aims 

states  (DOS)  but  both  used  the  band  results  only  as  studied  with  XPS  were  oxidized  by  heating  in  oxygen  to 

fMdn  to  interprete  the  dau,  leaving  a host  of  questions  6S0°  C for  two  hours  in  a tube  furnace  and  then  imme- 

unresolved.  The  nature  of  the  electronic  sutes  in  these  diately  transferred  to  the  spectrometer  vacuum.  The 

materials  has  not  therefore  thus  far  been  conclusively  compositions  of  the  films  u^  in  the  UPS  experiments 

determined  due  to  the  inability  of  theory  to  directly  and  were  determined  by  electron  microprobe  analysis  and 

unambiguously  reproduce  the  observed  spectra.  the  compositions  of  the  films  used  in  the  XPS  experi- 

In  the  present  paper  we  present  theoretical  calcula-  ments  were  estimated  from  the  relative  intemsities  of 

tkMS  of  the  photoemission  spectra  of  SiO„  GeO„  and  XPS  spectra  from  the  Ge  3d.  Si  2p,  and  O Is  core  lev- 

the  intermediate  composition  Si,Ge  ,0,,  which  repro-  els.  The  oxide  films  were  thicker  than  the  escape  depth 

duce  the  important  features  and  trends  in  the  observed  of  the  photoelectrons,  as  evidenced  by  the  absence  of 

spectra.  The  results  show  that  a molecular-orbital  pic-  unoxidized  Si  or  Ge  core  levels  in  the  spectra.  More 

ture  about  a single  oxygen  atom  is  inadequate.  They  experimental  details  are  given  in  Ref.  4.  The  expen- 

also  establish  the  importance  of  band  theory  for  these  mental  spectra  are  shown  in  Fig.  1. 

materials  and  show  that  quantitative  predictions  are  The  theoretical  calculations  were  carried  out  in 

feasible  for  materials  of  such  high  complexity  and  low  terms  of  the  empirical  tight-binding  model  employed  by 

symmetry.  Finally,  they  esublish  that  the  bandwidth  PH.  The  basis  orbitals  are  the  tetrahedrally-oriented  sp^ 

* Work  supported  in  pan  by  the  Office  of  Naval  Research,  Washington  D.C.  under  contract  No.  N 00014-76- 

C-0934  and  in  pen  by  the  Advanced  Research  Product  Agency  and  monitored  by  AFCRL  under  contract  No. 

Fm28-76-C-0249 

’ Present  address:  Max-Planck  Institut  fur  Festkorperforschung.  Stuttgart,  West  Germany. 
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Fig. 


Experimental  UPS  and  XPS 


spectra  for  Si,Ge,.,0;  films. 


hybrids  on  the  Si/Ge  atoms  and  the  p orbitals  on  the 
oxygens.  These  orbitals  (Fig.  2)  are  allowed  to  interact, 
giving  rise  to  five  types  of  orbitals:  (1 ) a bonding  orbital 
B,  comprized  of  p,  and  the  odd  combination  of  the 
hybrids,  (2)  an  orbital  B,  comprized  of  p,  and  the  even 
combination  of  the  hybrids  (B,«p,  for  the  ideal  cubic 
structure.  Fig.  2a),  (3)  a non-bonding  orbital  B,>p,, 
and  (4,S)  two  corresponding  antibonding  orbiuls  A, 
and  A,.  The  B’s  are  occupied,  giving  rise  to  the  valence 
bands  while  the  A’s  are  empty,  giving  rise  to  the  con- 
duction bands.  For  the  intermediate  Si,Ge  ,0,  compos- 
ition, the  main  additional’  effect  is  a mixing  between  B, 
and  B,. 


Fig.  2.  The  tight-binding  orbitals  employed  in  the  cal- 
culations. (a)  For  the  ideal  cubtc  struaure,  and  (b)  for 
the  general  case. 


The  precise  composition  of  the  various  A's  and  B’s, 
as  well  as  their  energy  positions  are  entirely  determined 
in  terms  of  two  parameters,  the  hybrid— p-orbital  inter- 
action Wj  and  the  hybrid-p-orbital  energy  separation 
W,.  These  have  been  fixed  by  making  use  of  the  optical 
spectrum  of  SiO,  (see  Ref.  2).  The  parameters  for 
GeO,  and  the  intermediate  compositions  are  then  deter- 
mined by  making  use  of  the  relative  hybrid  energies  of 
Si  and  Ge  and  the  d~^  rule,  where  d is  the  bond  length.’ 
With  the  basis  orbitals  just  defined,  a tight-binding  ener- 
gy band  calculation  is  carried  out  by  retaining  only 
nearest-neighbor  interactions.  This  introduces  two  new 
parameters,  a hybrid-hybrid  interaction  V,  (carried  over 
from  previous  work  on  the  tetrahedral  semiconductors) 
and  the  oxygen-oxygen  ppa  interaction  V,.  The  latter  is 
adjusted  to  the  observed  total  width  of  the  valence 
bands  of  SiO,.  For  GeO,  and  the  intermediate  compos- 
itions V,  is  scaled  with  the  d'’  rule.  The  valence  bands 
obtained  for  cubic  SiO,  with  these  choices  are  in  very 
good  accord  with  the  first-principles  bands  of  SF.  We 
are  now  in  position,  however,  to  bend  the  Si-O-Si 
chains  to  the  observed  angles  (144°)  and  repeat  the 
calculations  and  obtain  the  valence  bands  and  corre- 
ponding  DOS  for  the  various  cases  of  interest.  The 
latter  are  then  broadened  by  convoluting  with  a gaus- 
sian  of  width  O.S  eV.  This  broadening  turns  out  to  be 
extremely  important  as  it  eliminates  most  of  the  sharp 
spiked  structures  in  the  DOS  (see  Fig.  3)  which  pre- 
cluded PH  from  carrying  out  more  than  a qualitative 
analysis  of  the  observ^  spectra. 

The  final  broadened  DOS  curves  for  SiO,,  GeO, 
and  Si,,Ge,,02  (or  SiGeO.),  are  shown  in  Fig.  4 for 
direct  comparison  vvith  Fig.  I . In  view  of  the  fact  that 
the  calculation  did  not  include  photoemission  matrix 
elements,  the  comparison  must  be  limited  to  the  overall 
trends  among  the  various  compositions  and  to  the  posi- 
tions of  dips  and  peaks  in  the  individual  spectra.  No 
comparison  can  be  made  of  the  relative  heights  of 
peaks. 

It  is  seen  that  the  overall  trends  present  in  the  data 
are  reproduced  by  the  theory.  In  particular,  the  total 
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width  decreases  as  one  goes  from  SiO,  to  GeO,.  Also 
the  main  dip  at  about  4 eV  fills  up  slowly  and  moves  to 
higher  energies  in  both  the  experimental  and  theoretical 
curves,  as  shown  by  the  dashed  lines  in  Figs.  I and  4. 


Fig.  3.  Density  of  States  (DOS)  and  broadened  DOS 
for  the  valence  bands  of  SiO,. 


Fig.  4.  Theoretical  photoemission  spectra  for 
Si.Cei  .O,  films.  Compare  with  Fig.  I. 

For  more  detail  companson  between  theory  and 
experiment,  the  positions  of  the  peaks  and  dips  in  the 
theoretical  and  experimental  curves  of  SiO.  and  GeO. 
are  listed  in  Table  I.  For  the  purposes  of  this  Table,  the 
zero  of  energy  in  the  theoretical  curves  was  redefined 
by  the  technique  employed  in  obtaining  the  zero  of 
energy  in  the  expenmental  spectra,  namely  by  linear 
extrapolation  of  the  descending  curve.  The  main  discre- 
pancy in  both  materials  it  peak  1 which  appears  as  a 
double  peak  in  the  theoretical  curves.  The  agreement 
between  experiment*  and  theory  for  the  remaining  fea- 
tures if  better  than  13%. 


TABLE  1 

Experimental  and  theoretical  positions  (in  eV)  of 
peaks  and  dips  in  the  photoemission  spectra  of  SiO,  and 
GeO,.  The  numbering  refers  to  Fig.  I in  order  of  in- 
creasing binding  energy. 


Feature 

XPS 

UPS 

Theory 

sto, 
dip  1 

5.3 

5.3 

4.7 

dip  2 

7.6 

8.1 

8.5 

peak  1 

2.8 

2.6 

1. 2,2.5 

peak  2 

6.5 

6.4 

5.7 

peak  3 

9.8 

9.8 

9.9 

total  width 

11.5 

11.2 

11.7 

GcO. 


dip  1 

4.0 

4.0 

3.7 

dip  2 

6.6 

6.9 

7.0 

peak  1 

1.9 

1.8 

1. 2,2.9 

peak  2 

5.1 

5.0 

4.5 

peak  3 

8.7 

9.3 

8.3 

total  width 

10.4 

10.1 

10.0 

The  empirical  tight-binding  method  employed  in 
the  calculations  presented  above  allows  a detailed  anal- 
ysis of  the  results  in  terms  of  simple  physical  concepts. 
For  example,  it  has  long  been  suggested  that  the  top  of 
the  valence  bands,  in  particular  the  broad  feature  from 
0 to  about  4 eV,  arises  from  the  non-bonding  oxygen  p 
orbiuls.  whereas  the  remainder  of  the  bands  arises  from 
the  bonding  orbitals.  The  present  calculations  show  this 
separation  to  be  an  oversimplification  of  reality.  First 
note  that  the  Si-O  interaction  (W,)  is  responsible  for 
the  formation  of  the  bonding  orbiuls.  It  is  the  0-0 
interaction  V,  alone  that  causes  the  banding  of  the  non- 
bonding orbitals,  and  it  is  together  with  V,  that  cause 
the  banding  of  the  bonding  orbiuls.  The  same  V^,  how- 
ever, is  also  responsible  for  the  admixture  between  the 
bonding  and  non-bonding  orbitals,  the  net  result  being  a 
strong  hybridization  between  the  bonding  and  non- 
bonding bands.  The  nature  of  the  hybridization  is  par- 
ticularly interesting.  If  it  is  turned  off,  the  pure  non- 
bonding bands  overlap  considerably  with  the  pure  bond- 
ing ones.  When  the  hybridization  is  turned  on,  the  two 
sets  of  bands  push  strongly  on  each  other,  opening  up  a 
gap  at  4 eV  (Fig.3).  TTiis  is  opposite  to  what  a 
molecular-orbiul  picture  would  suggest,  since  the  pres- 
ence of  a gap  between  two  levels  would  be  viewd  as  a 
consequence  of  hybridization. 

In  conclusion,  we  have  presented  new  experimental 
spectra  which  help  esubiish  the  electronic  structure  of 
SiO,  and  GeO,.  Theoretical  calculations  in  terms  of  the 
empirical  tight-binding  method  (ETBM)  have  proved 
successful  in  reproducing  the  experimental  spectra  and 
establishing  the  validity  of  the  band  picture  for  these 
rather  complex  materials.  It  should  be  emphasized  that 
only  the  total  width  of  the  bands  of  SiO,  was  used  as  an 
input  in  performing  the  calculations.  The  internal  struc- 
ture in  the  spectrum  of  SiO,,  as  well  as  the  other  spectra 
in  their  entirety,  were  frnucttd  ana  found  to  be  in  good 
agreement  with  experiment.  This  demonstrates  the 
capability  of  the  ETBM  to  yield  quantitative  predictions 
for  materials  for  which  other  methods  may  not  be  too 
convenient  to  apply  due  to  low  symmetry  and  compli- 
cated lattices. 
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ABSTRACT;  We  show  that  inclusion  of  electron-hole  interactions  not  only 
produces  an  exciton  at  10.4  eV,  but  also  modify  the  interband  continuum  in  a 
substantial  way.  The  result  is  a spectrum  consisting  of  three  excitonic  peaks,  in 
good  agreement  with  the  observed  spectrum.  These  peaks  are  identified  in 
terms  of  the  bond-orbital  model  as  members  of  a crystal-field  split  Frenkel-type 
exciton  corresponding  to  the  free  oxygen  atom  2p  3s  transitions. 
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INTRODUCTION 

The  optical  absorption  spectrum  of  SiO,,  shown'  in  Figs,  la  and  2,  has  been  the 
subject  of  considerable  interest  and  controversy*.  Following  early  attempts  to 
interpret  the  spectrum,--^  DiStefano  and  Eastman"*  concluded  from  a series  of 
photoconductivity  and  internal  phoioemission  measurements  on  Si-SiO^  structures 
that  the  band  gap  is  about  9 eV.  These  results  raised  the  question  why  no  strong 
excitons  appear  in  the  absorption  spectrum  below  the  band  edge  as  is  the  case  in 
other  ionic  wide-gap  materials,  such  as  the  alkali  halides.  The  observed  peaks  thus 
became  the  subject  of  conflicting  interpretations  in  terms  of  excitons.  interband 
transitions  and  various  combinations  thereof  (see  e.g.  Fig.  6 in  Ref  2 for  summary). 
More  recently,  Schneider  and  Fowler^  calculated  a band  structure  for  cubic  SiO,  and 
suggested  that  excitons  are  not  formed  because  of  a symmetry-forbidden  edge. 
Chelikowsky  and  Schluter  (CS)*’  actually  calculated  an  interband  spectrum  (Fig.  Ic) 
and  found  an  indirect  bandgap  of  9.2  eV,  which  is  consistent  with  the  measurements 
of  Ref.  4.  They  also  found  that  matrix  elements  suppress  absorption  below  ~10.4 
eV.  They  therefore  suggested  that  the  first  observed  peak  corresponds  to  an  exciton. 
as  previously  proposed  by  others,**^  while  the  rest  of  the  spectrum  is  due  to  inter- 
band transitions.  In  support  of  the  latter  assignment,  they  displayed  the  experimental 
spectrum  shifted  by  0.7  eV  and  noted  the  resulting  good  agreement  between  theory 

^Supported  in  part  by  the  Office  of  Naval  Research  under  Contract  No.  N00014-76-C- 
0934. 
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and  experiment  (Fig.  lb  and  Ic).  More  recently,  Mott’  investigated  the  subject  and 
suggested  that  there  exists  an  a'lowed  exciton  at  9 eV  which  has  a dispersion  of  ~2.4 
eV  due  to  hopping,  resulting  in  the  10.4  eV  peak  (see  also  the  paper  by  Mott  in 
these  Proceedings). 


Fig.  1 (a)  The  optical  absorption  spectrum  of  SiO,  meas- 

ured by  Klein  and  Chun  [Ref.  1 ; the  spectrum 
shown  is  actually  E“e,(E)].  (b)  The  same  spectrum 
shifted  by  0.7  eV  in  Ref.  4 in  order  to  align  it  with 
the  theoretical  interband  spectrum,  curve  (c). 


In  this  paper  we  present  two  completely  independent  treatments  of  electron-hole 
(e-h)  interactions  and  conclude  that  the  observed  peaks  are  a series  of  excitonic 
resonances.  This  interpretation  is  consistent  with  but  goes  beyond  a band-theoretic 
approach  and  is  also  consistent  with  the  known  data. 


BAND-THEORETIC  APPROACH  AND  THE  EFFECT  OF  E-H 
INTERACTIONS 
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In  most  semiconductors  which  have  large  dielectric  constants,  e-h  interactions  tend  to 
be  weak  and  produce  excitons  with  very  small  binding  energies  and  oscillator 
strengths.  As  a result,  the  interband  spectrum  remains  essentially  unmodified  and  is 
in  general  well  reproduced  by  band-structure  calculations.  In  the  case  of  widc-gap 
insulators,  however,  dielectric  constants  are  small,  making  e-h  interactions  very 
strong.  The  resulting  excitons  are  then  localized  and  carry  considerable  oscillator 
strength.  This  oscillator  strength  is  pulled  out  of  the  interband  continuum  (the  total 
oscillator  strength  remains  constant)  in  such  a way  that  the  final  spectrum  is  substan- 
tially different  from  the  independent-particle  spectrum.  This  assertion  has  been 
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supported  by  theoretical  calculations*-’  and  by  a study  of  39  x-ray  absorption  spectra 
of  alkali  halides.'" 

In  this  section  we  start  with  the  interband  spectrum  of  SiO,  calculated  by  CS^',  which 
indicated  that  the  10.4-eV  peak  should  not  be  attributed  to  interband  transitions. 
We  will  then  incorporate  e-h  interactions  by  means  of  a model  calculation,  which 
indicates  that  e-h  interactions  strong  enough  to  produce  an  excitonic  peak  at  10.4  eV 
also  modify  the  interband  spectrum  very  strongly,  so  that  all  the  observed  peaks  are 
excitonic  in  nature. 

The  model  calculation  we  intend  to  carry  out  is  similar  to  that  done  previously’  for 
LiF.  We  start  with  the  interband  spectrum  calculated  by  CS  and  assume  it  can  be 
simulated  by  a two-band  model  and  momentum-independent  transition  matrix 
elements.  In  such  a model,  the  absorption  coefficient  for  interband  transitions  is 
given  by  the  imaginary  part  of  a two-particle  Green’s  function  G°  describing  inde- 
pendently propagating  electrons  and  holes.  We  therefore  construct  an  analytical 
form  for  G°(E)  as  in  Ref.  9,  chosen  so  that  its  imaginary  part,  when  broadened  by 
0.5  eV,  reproduces  CS’s  interband  spectrum.  E-h  interactions  are  then  introduced  in 
the  contact  approximation,*  ’ i.e.  as  a single  on-site  matrix  element  V in  the  Wannier 
representation.  The  new  absorption  spectrum  is  given  by  the  imaginary  part  of  the 
new  Green’s  function  G,  which  satisfies  Dyson’s  equation: 


G - G° -I- G°VG.  (1) 


Eq.  ( 1 ) is  solved  directly  to  yield 


G(E)  = G»(E)/[1  - G‘»(E)V]  (2) 


so  that  a bound  state  occurs  when  G°(E)  = 1/V.  The  value  of  V was  chosen  to 
yield  a bound  state  at  10.4  eV. 

The  results  of  the  calculation  are  shown  in  Fig.  2.  The  dashed  curve  is  — Im  G"(E), 
with  E having  a 0.5  eV  imaginary  part  for  broadening,  and  corresponds  to  CS’s 
interband  spectrum  (cf.  with  Fig.  Ic).  The  solid  curve  marked  theory  is  — Im  G(E), 
also  broadened  by  0.5  eV.  These  results  clearly  show  that  if  the  10.4-eV  peak  is 
indeed  an  exciton,  e-h  interactions  also  modify  the  interband  spectrum  strongly,  so 
that  the  final  spectrum  is  excitonic  in  nature,  i.e.  consists  of  a series  of  excitonic 
resonances  (sometimes  called  metastable  excitons).  Note  that  whereas  the  two  peaks 
in  the  interband  spectrum  are  separated  by  ~2.3  eV,  the  two  peaks  in  the  final 
spectrum  are  separated  by  1 .4  eV,  in  excellent  agreement  with  experiment.  The 
spectrum  at  higher  energies  is  not  reproduced  very  well  as  one  might  expect,  since 
oscillator  strength  at  those  energies  would  have  to  be  brought  down  by  e-h  interac- 
tions from  even  higher  energies,  which  are  left  out  by  the  present  model.  Neverthe- 


less,  the  calculated  peak  at  15  eV  compares  rather  well  with  the  observed  peak  at  14 


eV. 
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Fig.  2 Theoretical  absorption  spectrum  compared  with  the 
optical  absorption  spectrum  [(^(E)]  measured  by 
Philipp  (Ref.  1).  The  dashed  line  is  interband 
spectrum  before  e-h  interactions.  See  text. 

Apart  from  the  success  of  the  above  model  calculation  in  predicting  the  energy 
separations  and  relative  heights  of  the  peaks  in  the  optical  spectrum,  the  main 
conclusion  is  that  in  wide-gap  insulators  one  cannot  ignore  the  effect  of  e-h  interac- 
tions on  interband  transitions,  as  done  in  semiconductors,  in  contrast,  e-h  interac- 
tions have  a rather  pervasive  effect,  so  that  excitonic  peaks  dominate  the  entire 
spectrum.  One  should  not,  therefore,  hope  to  extract  a band  gap  from  such  a 
spectrum. 


BOND-ORBITAL  APPROACH 

The  above  calculation  indicates  that  the  absorption  spectrum  of  SiO.,  can  be  inter- 
preted as  a series  of  excitonic  peaks  arising  from  localized  excitations,  but  dues  not 
provide  information  about  the  nature  of  states  contributing  to  the  two  ob.served 
peaks.  As  was  the  case  for  alkali  halides,'"  an  atomistic  point  of  view  can  he 
complementary  and  useful.  For  this  purpose,  we  make  use  of  the  tight-binding 
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bond-orbital  model"  for  SiOj,  according  to  which  the  valence  band  states  arise  from 
three  different  bond  orbitals,  associated  with  a Si-O-Si  bonding  unit:  which  is  a 

pure  O 2p  orbital  perpendicular  to  the  Si-O-Si  plane  (lone-pair  or  non-bonding  p 
orbital);  which  is  ~85%  an  O 2p  orbital  in  the  Si-O-Si  plane  perpendicular  to 
the  Si-Si  axis  and  -15%  Si  sp^  hybrid  (partially  bonding  orbital);  and  B^,  which  is 
-65%  an  O 2p  orbital  along  the  Si-Si  axis  and  -35 ‘K)  Si  sp^  hybrid  (bonding 
orbital).  These  three  orbitals  at  each  bonding  unit  have  been  found  adequate  for  a 
calculation  of  the  valence  bands  and  the  related  photoemission  and  x-ray-emission 
spectra." 

According  to  the  same  bond-orbital  scheme,  the  lowest  unoccupied  band  states  arise 
from  a bond  orbital  denoted  in  Ref.  10  by  A^,  consisting  mainly  of  O 3s  and  Si 
hybrid  admixtures.  We  suggest  that  during  photoabsorption  the  electron  and  the  hole 
do  not  occupy  propagating  band  states.  Instead,  e-h  interactions  tend  to  localize  the 
electron  and  the  hole  in  the  same  bonding  unit,  so  that  excitonic  peaks  are  expected 
at  the  By  -*•  A^,  -►  A^,  and  B^  A^  energies,  in  that  order.  Note  that  these 

transitions  are  essentially  localized  on  an  O atom  and  correspond  to  the  atomic 
transition  2p  ■*  3s.  In  the  free  atom  the  transition  energy  is  about  9.2  eV.  In  the 
solid,  By,  Bjj,  and  B^  may  be  viewed  as  the  "crystal-field"  split  members  of  the  O 2p 
orbital,  resulting  in  a "crystal-field"  split  Frenkel-type  exciton  with  peaks  at  10.4  eV. 
11.7  eV,  and  -14  eV.  These  energies  cannot  be  predicted  by  some  simple  quantita- 
tive calculation,  but  the  interpretation  may  be  further  supported  by  a variety  of 
arguments.  First,  CS’s  calculation  showed  that  the  10.4-eV  should  not  be  attributed 
to  interband  transitions.  Second,  the  calculation  of  the  previous  section  demonstrated 
that  if  the  10.4-eV  peak  is  indeed  an  exciton,  the  other  peaks  are  also  predominantly 
excitonic  in  nature.  Third,  using  the  identification  presented  above,  one  can  infer 
that  the  three  peaks  ought  to  be  successively  shorter,  because  of  the  smaller  weight 
of  O 2p  in  each  successive  B,  and  also  broader,  because  of  the  shorter  lifetime  of  a 
hole  in  each  succesive  B.  Both  these  observations  are  in  agreement  with  experiment 
and  consistent  with  the  calculation  of  the  previous  section. 
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Recent  Advances  in  the  Theory  of  Electronic  Structure  of  SiO: 


Silicon  dioxide  (SiO;)  is  a mutcriul  that  has  attracted  the  attention  of  scientists 
for  a long  time.  It  is  available  in  nature  in  a variety  of  forms  and  has  been  used 
extensively  by  the  construction  industry,  by  tiic  glass  industry,  and.  during  the 
last  fifteen  years,  in  technological  applications,  m particular  the  manufacture 
of  metal-oxide  semiconductor  (MOS)  devices.  There  exists  a vast  experimental 
and  theoretical  literature  on  its  physical,  chemical,  electronic,  optical,  thermal, 
transport,  magnetic  and  other  properties. 

Theoretical  understanding  of  the  electronic  properties  of  SiO;  began  vvitii 
qualitative  molecular-orbital  pictures  aimed  at  the  interpretation  of  particular 
experimental  spectra,  such  as  the  optical  absorption  st;cct.’'a.  the  x-ray  emission 
spectra,  etc.  These  early  models  were  reviewed  by  Ruffa.'  who  identified  '.heir 
relative  merits  and  shortcomings.  New  insights  were  gamed  by  subsequent  work 
by  several  workers’  who  quantified  such  studies  with  scIf-consistcnt  eaicui.itious 
on  small  clusters,  but  the  ambiguities  clouding  the  connection  between  the 
discrete  energy  levels  of  the  cluster  with  the  continuous  excitation  spectra  of 
tile  solid  were  a persistent  problem.  In  particular,  clusters  ol  inc.'casing  sbe  were 
not  found  to  yield  a converging  picture  of  i-:vci  densities.  Furthermore,  ciusters 
centred  about  different  atoms  (Si  or  0)  were  needed  to  interpret  the  x-rcy 
emission  spectra  associated  with  core  levcis  on  different  atoms,  so  that  a con- 
sistent and  universal  level  sciiciuc  was  elusive. 

During  tltc  last  two  years  several  advances  liavc  been  made  in  our  understand- 
ing of  the  electronic  properties  of  SiO;,  in  paiiicniar  from  applications  of  banc 
theory.  In  view  of  the  technological  importance  of  Sit);,  it  is  rather  surpri.\ing 
that  band  calculations  have  been  lackirc.  It  is  therefore  interesting  first  to 
explore  the  reasons  for  this.  Wo  will  then  examine  the  recent  accomplishments 
of  band  theory  and  assess  the  current  and  future  objectives  of  theoretical  work 
in  SiOj. 

The  reasons  why  energy  band  calculations  have  not  been  available  for  SiO; 
lie  in  its  crystal  structure,  or,  rather  crystal  siruciiircs.  SiO;  e.xists  in  natu.-c  in  a 
variety  of  crystalline  forms,  cailcd  polymorphs,  ali  of  which  have  low  symmetry 
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and  large  unit  cells.  For  example,  a-quartz,  the  best  known  polymorph,  has  nine 
atoms  in  the  primitive  unit  cell.  By  comparison,  some  of  t.he  materials  that  are 
studied  most  oricn  have  one  or  two  atoms  per  unit  cell  (simple  and  trenstticn 
metals,  alkali  halides,  tetraiiedral  semiconductors,  etc.).  The  consequence  of 
many  atoms  in  the  primitive  unit  ceil  is  that  many  more  bands  are  .needed  for  a 
complete  description.  Fur  example,  a-quartz  has  24  valence  bands,  compared 
with  1 , 3 or  4 in  the  materials  mentioned  above.  Tlie  large  number  of  bands  calls 
fur  very'  large  basts  sets  for  proper  convergence,  which  might  strain  even  large 
computers  or  cumputei  ''-idgets.  Apart  from  the  problem  of  computer  time, 
oii'cr  Iturdics  liavc  ais«>  been  present.  Fur  example,  the  construction  of  a puten* 
tial  for  widc-gap  insulators,  especially  compounds  involving  a first-row  element, 
such  as  oxides,  niindcs.  fiuoridcs,  etc.,  is  known  to  involve  large  uncertainties. 
Calculations  on  such  materials  using  lucal-dcnsity  approximations  to  the  ex- 
cltange  and  correlation  potentials  in  general  procuce  bandwidths  arc  bandgaps 
which  are  too  narrow;  pure  Hartree-Fock  calculations  give  the  apposite:  on  the 
other  hand,  local-density  approximations  have  been  quits  successful  for  the  less 
ionic  compound  tetrahedral  semiconductors  so  t.hat  SiOj  might  ce  though:  to  bs 
a borderline  material. 

How  w'ere  these  problems  circumvented?  The  firs:  two  calculations^*'^*  were 
done  on  the  ideal  cubic  fonn  of  SiO;.  which  has,  in  the  past,  been  tnough:  to  ba 
the  structure  of  the  polymorph  known  as  y-cristoba!ite.  This  structure  has  six 
atoms  per  primitive  unit  cell,  but  the  symmetry  is  rather  high,  namely  the  same 
as  that  of  diamond  (space  group  0),’ ).  The  atomic  arrangement  is  actually  quite 
interesting,  in  that  tlie  Si's  form  a diamond  lattice  (just  like  Si  itself,  except  for 
a somewhat  larger  lattice  spacing),  and  the  oxygens  lie  at  the  mid-points  of  the 
lines  connecting  nca.''cst-nctgfib'jr  St's,  the  so-called  bond  sites.  Tims,  each  Si 
is  surrounded  by  four  oxygen  atoms  in  tetrahedral  directions,  and  each  oxygen 
has  two  Si  nearest  neighbors,  just  like  in  most  real  forms  of  SiO} . The  only 
exception  is  that  the  Si-O-Si  angle  is  180“’  instead  of  the  observed  I-i0“-I50°. 
Parenthetically,  we  may  note  that  the  o-cristobalite  structure  can  be  arrived  at  by 
starting  with  the  cubic  structure  just  described  and  bending  the  Si-O-Si  chains 
in  a systematic  way.  A tetragonal  distortion  then  occurs.  In  a similar  fashion, 
the  structure  of  tridymite  can  be  arrived  at  by  starting  with  a wurtzite  lattice  of 
Si's  (instead  of  a diamond  lattice  of  Si's)  and  then  inserting  the  o.xygens  between 
the  Si's  and  bending  the  Si-O-Si  chains. 

The  first  band  calculation  was  donc^  for  p-cnstobalitc  in  the  ‘empirical  tight- 
binding  method'  (ETBM)  winch  employs  a tight-binding  basis  set  and  treats  the 
matrix  elements  as  parameters.  This  technique  resembles  the  often  used  pro- 
cedure of  obtaining  a tight-binding  yi;  to  energy  bands  calculated  by  a sophisti- 
cated firs- •principles  method.  The  fundamental  difference  is  that  in  the  UTBM 
one  proceeds  in  the  absence  of  any  band  calculation  whatever,  and  determines 
the  values  of  the  parameters  by  scaling  parameters  from  other  materials  and/or 
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experimental  data  for  the  material  of  interest.  It  is  thus  analogous  to  the  empirical 
pseudoporential  method  (EP.M)  which  is  also  occasionally  used  as  an  interpolation 
scheme.  A particular  virtue  of  both  the  EPM  and  the  ETBM  is  that  tlicy  are  tests 
of  the  single-particle  model  to  interpret  the  excitation  spectra  of  solids  without 
depending  upon  the  ability  of  contemporary  theory  to  construct  from  first- 
principles  the  one-electron  Hamiltonian,  in  particular  the  self-energy  opciator 
which  describes  the  complicated  electron-electron  interaction.  On  the  other 
hand,  a serious  drawback  of  the  ETBM  is  that  it  produces  only  valence  bands. 
Later,  more  laborious  calculations,  using  first-pnnciples  potentials,  by  Schneider 
and  Fowler”  and  by  Ciraci  and  Batra*  produced  similar  valence  bands  which 
confirmed  the  validity  of  the  assumptions  made  ir.  the  ETBM.^  The  latter 
calculations  were  more  complete  in  that  they  also  produced  conduction  hands. 
Calculations  on  cubic  SiOj,  however,  do  not  necessarily  rciaie  directly  with 
expcriinentai  data  because  ail  the  .-eal  polymorphs  of  SiOj  have  hcni  Si-O-.Si 
chains,  with  angle  at  the  0 atoms  about  144*.  compared  with  the  180"  for  cubic 
SiOj.  Nevertheless,  the  energy  bands  of  cubic  SiOj  help  draw  interesting 
analogies  with  the  well-known  bands  ofSi,^  and  suggest  qualitative  interpre- 
tations*’* of  the  spectra*’*  of  SiOj. 

Experimental  optical  spectra  available  for  amorphous  SiOj  and  for  cry'stallinc 
quartz  We'Ve  known  to  be  identical,''  which  led  to  the  conclusion  that  the  local 
coordination  determined  all  the  structure,  with  long-range  order  being  immaterial. 
Pantelides  and  Harrison^  then  assinncd  bent  Si-O-Si  cliains  in  a hypothetical 
diamond  lattice,  a model  which  simulates  the  local  topology  of  amorphous  SiO;, 
and  obtained  valence  bands  and  density  of  states.  Pantclidcs^'’^^  later  extended 
this  to  the  calculation  of  tlic  various  ,\-ray  emission  spectra.  After  appropriate 
broadening  was  inciuded.  the  calculated  photocmission  ami  x-ray  emission 
spectra  agreed  in  a quantitative  way  with  experiment.  Similar  calculation  also 
reproduced  the  trends  in  the  cxpcnincntal  spectra  of  SivGc,.^0:  alloys.* 

More  recently,  .Schlutcr  and  Chclikowsky*  (SC)  and  Calabrese  and  Fowler* 
(CF)  reported  energy  band  calculations  for  a-quarf/..  SC  employed  parametrized 
pseudopotcntials  fur  the  core  potentials  and  carried  out  a self-consistent  calcu- 
lation for  the  valence  electrons  in  a plane-wave  representation.  CF  used  a super- 
position of  atomic  potentials  and  a mixed  basis  set.  Similar  results  were  obtained, 
but  only  SC  were  able  to  calculate  enough  points  in  tiic  Briilouin  zone  to 
obtain  explicit  spectra 

Schlutcr  and  Chclikowsky  calculated  detailed  photocmission  and  x-ray 
emission  spectra  that  agreed  very  well  with  ex|H’riincnt.  Furthcrinorc.  llicy 
obtained  detailed  charge  density  plots  which  cunfirincd  the  basic  assumptions 
made  earlier' about  the  nature  of  the  bunding  and  non-bonding  orbitals.  For 
the  first  time,  however.  SC  also  calculated  a theoretical  absorption  spcct.-um 
for  SiO].  The  calculation  was  limited  to  interband  transitions,  but  it  immediately 
revealed  that,  even  though  absorption  begins  at  tlic  indirect  gap  of  about  9.2  cV. 
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it  rcnuiiiN  essentially  zero  I'ur  aliout  I cV  and  finally  peaks  at  1 1 cV  ami  then 
again  at  LV.<  eV.  SC  iJcntifieil  the  first  experimental  peak  at  10.4  cV  as  an 
cxciton.  in  agreement  with  previous  assignments,  ami  attnbuted  the  structure  at 
higlicr  enereies  to  intcrbaml  transitions.  T!ie  interpretation  v^as  suggested  by 
the  separation  of  the  tncorctical  peaks  being  the  same  as  the  separation  of  the 
experimental  peaks  at  1 1.7  and  14.0  cV.  One-to-one  correspondence  between  the 
experimental  and  theoretical  spectra  was  obtained  by  siiifting  the  two  spectra 
relative  to  each  other  by  0.7  cV. 

Another  property  that  has  recently  been  calculated  in  addition  to  the  spectra 
mentioned  above  is  the  dielectric  constant.  This  was  done  by  1‘anlclides  and 
Harrison^  who  obtained  a rather  simple  analytical  expression,  which  depended 
explicitly  on  the  value  of  the  Si-()-Si  angle  and  other  known  quantities.  A new 
unknown  parameter  of  order  unity  was  introduced,  but  once  its  value  was  fixed 
for  a-quart/.  the  dielectric  constanis  for  the  other  SiO;  polymorphs  could  be 
predicted,  and  agreement  with  experiment  was  good.  The  most  interesting 
result**^  was  for/J-cristobalite.  if  the  ideal  (cubic)  value  of  ISO”  for  the  Si-O-Si 
angle  was  used,  the  calculated  dielectric  constant  t was  rixa  large  by  When 
an  angle  of  144“  (the  Q-quartz  angle)  was  n.sed,  c became  identical  to  the  experi- 
mental value.  Similarly,  for  cocsitc.  if  the  value  of  I Z0“  was  used  (as  listed  by 
Wyckoff'®),  e was  too  small  by  25ri.  Again,  use  of  144“  brought  e in  line  with 
experiment.  Since  /^cri$tobalitc  and  cocsitc  are  the  only  polymorphs  thougiit 
to  have  angles  significantly  different  from  I44“,  these  results  suggest  that  all 
SiOj  polymorphs  may  have  Si-O-Si  angles  of  about  144“.  These  results  also 
reveal  the  amount  of  detailed  information  that  can  be  extracted  from  a para- 
metrized method.  Studies  of  the  dependence  of  spectra  and  properties  on  the 
• Si-O-Si  angle  would  otherwise  be  very  involved  calculations. 

\Vli.n  tlien  have  band  structure  calculations  accomplished'’  For  the  first  time, 
a consistent  picture  of  bonding  in  SiO]  has  been  obtained.  It  has  been  shown 
that  the  details  of  the  photocmission  spectra  and  the  various  x-ray  emission 
sjKctra  can  be  systematically  reproduced  only  via  an  energy,'  band  calculation, 
and  througlit  the  calculation  of  densities  of  states  with  and  without  matrix 
elements,  a procedure  that  requires  knowledge  of  the  bands  at  many  points  in 
the  Brilloiiin  /one.  It  has  also  been  shown  that  dielectric  constants  can  be 
understood  in  detail.  Finally,  progress  has  been  made  In  elucidating  the 
optical  absorption  spectrum  of  SiO],  but  that  has  been  limited  by  the  fact 
that  electron-hole  interactions  cannot  yet  be  conveniently  included. 

Where  do  we  go  from  here?  Well,  much  remains  to  be  done.  To  begin  with, 
the  available  energy  bands  don’t  agree  with  each  miter  in  detail.  For  example, 
Calabrese  and  Fowler  obtained  u direct  gap  for  a-quartz,  wheras  Schluter  and 
CItclikiiwsky  find  an  indirect  gap.  Furthermore,  effective  masses  have  not  been 
calculated  accurately.  (Sclilutcr  and  Chcli’xowsky  obtain  0.3mo  for  the  con* 
duction  band  mass,  whereas  Schneider  and  Fowler  estimated  O.Smo-)  These 
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qujntitics,  as  well  as  the  interaction  of  electrons  and  liuies  with  the  lattice 
for  the  formation  of  polarons.  are  important  in  theories  of  transport  pro- 
perties and  tunneling,  whicli  are  still  not  very  well  understood. 

Puzzling  questions  still  remain  about  the  optical  spectra.  What  is  the  nature 
of  the  10.4  cV  exciton,  in  view  of  the  fact  that  it  is  above  the  threshold  for 
interband  transitions?  Is  it  possible  that  the  formation  of  a strong  exciton  leaves 
the  interband  spectrum  unchanged  in  its  overall  shape,  as  suggested  by  SC's 
analysis?  Similar  suggestions  have  been  made  previously  in  the  case  of  core 
cxcilons  in  alkali  halides,  but  subsequent  analysis*  showed  that  electron-hole 
interactions  in  malcriais  with  low  dielectric  constants  not  only  give  rise  to 
liighly  localized  cxcitons,  but  also  perturb  tltc  interband  spectrum  in  a sub- 
stantial way. 

Some  uncertainties  also  still  remain  in  the  photoemission  and  X-ray 
emission  spectra.  Both  Schluter  and  Chclikowsky’  and  the  present  author^** 
do  very  well  for  the  Si  K spectrum.  Sciiluter  and  Clielikowsky  also  do  very 
well  for  the  0 K spectrum  but  fail  to  reproduce  one  of  the  two  peaks  in  the 
Si  Lj,j  spectrum,  which  they  attribute  to  possible  sample  contamination.  The 
present  author,  on  the  other  hand,  obtains  a good  Si  spectrum  but  a less 
satisfactory  0 K spectrum.  As  for  the  x-ray  absorption  spectra,  they  arc 
virtually  not  understood  at  all.  Similarly,  theories  of  defect  and  impurity 
states  are  in  primitive  stages." 

Detailed  understanding  of  SiO-  can  also  lead  to  understanding  of  a large  class 
of  materials  for  which  SiO:  is  a prototype.  We  recall  that  the  structure  of  SiO: 
may  be  understood  by  starting  with  a lattice  of  Si  atoms  and  inserting  o.xygen 
atoms.  One  could  also  insert  oxygens  in  a zince-blcndc-typc  lattice,  such  as 
CaAs  or  CaP.  The  result  is  then  Ga.AsO*,  GaP04,  etc.  This  class  of  materials 
may  be  described  by  the  general  formula  ABO4,  where  AB  Is  a zince-blcnde-iypc 
semiconductor.  SiOj  is  then  a special  case  of  the  AB04-type  oxides  with 
A ■ B ■ Si  (SiOj  ■ SiSi04,  or  silicon  silicate!)  and  is  the  prototype  of  the  whole 
class.  Just  as  Si  is  the  prototype  of  the  tetrahedral  AB-type  semiconductors.’' 

Finally,  an  area  where  detailed  tlicurcticai  calculations  are  needed  is  the 
Interface  between  SiO:  -"J  S*.  as  well  as  the  interfaces  between  SiO;  and 
various  iiielals.  Tlie  properties  of  these  interfaces  are  cfieial  to  ilic  '.K'iiaviour  of 
MOS  devices  and  have  long  been  studied  by  experimental  •.cclmiques.  One  of 
the  outstanding  problems  is  the  thickness  of  such  interfaces.  For  example, 
the  question  whether  an  SiO.rCO  <x<2)  layer  c.xists  between  Si  and  SiO; 
has  attracted  considerable  attention.  Experimental  estimates  of  the  thickness  cf 
the  SiO,  layer  range  from  less  than  4 A to  1 5-20  Panteiidcs  and  long”  recently 
constructed  a continuous  random  network  model  which  shows  that  a zero- 
width  interface  is  possible  for  the  (100)  Si-SiO;  interface,  but  the  question 
still  remains  open.  This  uncertainty  complicates  even  further  the  diffieult  task  of 
performing  calculations  of  interface  state.  Thus  far,  oniy  qualitative  modcis 
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have  been  developed. '*  Some  quantitative  calculations  were  recently  reported  by 
Ciraci  and  Baira,”  but  the  field  is  wide  open,  and  .d  indications  are  t-tat 
rigorous  calculations  should  become  possible  in  t-ne  immediate  future. 

In  conclusion,  considerable  progress  has  been  achieved  i.".  our  understanding 
of  the  electronic  properties  of  SiOj . Band-theoretic  techniques  have  now  been 
applied  and  have  already  resolved  many  outstanding  problems.  Continued  work 
along  the  same  lines  ought  to  supply  even  more  answers.  Finally,  mere  speciaiiaei 
techniques  are  still  needed  for  the  understanding  of  the  eicctronic  proaerties 
of  the  interfaces  of  SiO- . The  tuture  looks  promising  for  this  important 
material. 
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ABSTRACT 

Ion  implantation  has  been  used  to  incorporate  As  and  P into  the  thermally  grown  SiO^ 
layer  of  metal-silicon  dioxide-silicon  structures.  These  impurities  increase  the  electron  trap 
density  in  the  oxide  layer  proportionally  to  the  ion  fluence  (1-3  x 10'^  cm'^).  Avalanche 
injection  from  the  silicon  substrate  was  used  to  populate  the  electron  trapping  sites.  It  is 
shown  that  the  negative  charge  is  removable  from  the  trapping  centers  under  illumination  with 
photons  between  3 and  5.6  eV.  From  the  detrapping  experiments,  we  detemined  a photoioni- 
zation threshold  of  % 4 eV  both  for  As  and  P related  centers,  and  a spectrally  resolved 
effective  photoionization  cross  section;  the  latter  quantity  is  defined  as  the  convolution  of  the 
photoionization  cross  section  with  the  optically  accessible  trap  distribution  in  the  energy  gap  of 
the  Si02.  A simple  model  is  presented  which  suggests  that  the  observed  detrapping  originates 
from  levels  corresponding  to  substitutional  P and  As  at  O sites  in  the  Si02. 
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INTRODUCTION 

Photodepopulation  spectroscopy  in  metal-insulator-semiconductor  (MIS)  structures  has 
recently  been  used  to  study  both  the  spatial  and  energy  distribution  of  optically  accessible 
charges  in  the  insulator  film  [1],  Several  variations  of  the  technique  exist,  all  relying  upon  the 
optical  stimulation  of  electrons,  trapped  in  deep  levels  in  the  insulator,  into  cqnducting  states 
from  which  they  drift  towards  an  electrode. 

For  MIS  structures  with  Si02  as  the  insulator,  literature  is  available  on  a 2.4-2. S eV  deep 
electron  trap,  possibly  related  to  sodium  incorporated  at  elevated  temperatures  in  an  immobile 
configuration  in  the  oxide  film  [2].  Few  studies,  however,  have  been  carried  out  to  character- 
ize other  charge  trapping  centers.  Several  ions  such  as  phosphorus  and  arsenic  were  recently 
demonstrated  to  exhibit  electron  trapping  characteristics  when  implanted  into  the  Si02  layer 
[3].  It  was  determined  that  P implantation  results  in  a dominant  electron  trap  with  a capture 
cross  section  of  3x10'*^  cm^  [3],  whereas  for  As  related  centers  the  dominant  trap  has  a 
capture  cross  section  of  » 10''^  cm^  [4].  The  trapping  sites  are  in  a net  neutral  state  before 
electron  capture.  In  both  cases,  the  integrated  trap  density  increases  with  the  ion  fluence. 
Good  agreement  was  found  between  the  trapped  electron  distribution  centroid  and  both 
theoretical  predictions  and  experimental  observations  of  the  ion  distribution  centroid  [3]. 

In  this  work,  we  report  on  spectrally  resolved  detrapping  measurements  for  charged  MOS 
structures,  from  which  a photoionization  threshold  and  an  effective  photoionization  cross 
section  spectrum  were  determined.  We  also  present  a simple  theoretical  model  which  suggests 
that  the  observed  levels  are  those  of  substitutional  P and  As  at  O sites. 
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EXPERIMENTAL 

Sample  Fabrication 

The  starting  material  was  p-type  <100>  silicon  with  a resistivity  of  0. 1-0.2  Qcm.  The 
wafers  were  oxidized  at  1000‘’C  in  a "dry"  oxygen  ambient  to  oxide  thicknesses  ranging  from 
560  to  1430  A as  determined  by  ellipsometry.  Then  P'*'  or  As'*'  implantation  wa*s  performed  at 
if  room  temperature  with  energies  of  20  to  80  keV  and  fluences  of  I to  3 x 10'^  cm'^.  The  ion 

I current  at  target  during  implantation  was  of  the  order  of  1 x 10*^  A.  All  wafers  were  then 

cleaned  and  annealed  in  nitrogen  at  lOGO^C  for  30  min.  Using  a shadow  mask,  semitranspar- 
ent aluminum  electrodes  (lOO-ISO  A thick)  with  an  area  of  0.0052  cm^  were  deposited  to  form 
MOS  capacitors.  Finally,  all  devices  were  given  a post-metallization  annealing  treatment  at 
400°C  in  forming  gas  for  20  min. 

Electron  Trapping  and  Photodetrapping 

1 

j 

I The  oxide  traps  were  charged  by  avalanche  injection  of  hot  electrons  from  the  silicon 

I substrate  driven  into  deep  depletion  [S].  The  amplitude  of  the  50  kHz  ramp  wave  used  for  this 

I purpose  was  constantly  adjusted  in  order  to  keep  the  average  dc  injection  current  constant. 

The  currents  ranged  from  5 x 10''®  to  3 x 10"’  A.  Some  of  the  injected  electrons  were 

I 

trapped  in  the  Si02,  causing  a shift  in  flat-band  voltage,  determined  from  1 MHz  differential 
capacitance-voltage  characteristics. 

After  reaching  a given  charge  level,  the  avalanche  injection  was  stopped  and  the  sample 
mounted  in  a set-up  for  photoelectric  measurements,  consisting  of  a 900  W xenon  arc  lamp  in 
combination  with  a 500  mm  grating  monochromator  (Bausch  and  Lomb).  The  incident  photon 
flux  at  the  sample  position  was  measured  over  the  spectrum  using  a thermopile  in  combination 
with  an  electronic  chopper.  An  electric  shutter  was  used  to  control  the  illumination  time  which 
was  usually  5 min.  The  change  in  flat-band  voltage  induced  by  this  illumination  at  room 
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temperature  was  monitored  to  within  1 mV  with  an  automatic  tracking  system.  The  spectrum 
was  scanned  step-by-step  in  the  direction  of  increasing  photon  energies. 

The  detrapping  experiments  are  performed  with  zero  gate  bias  in  order  to  avoid  injection 
from  the  contacts  at  photon  energies  greater  than  the  Si-SiOj  or  Al-Si02  energy  barrier,  by 
using  the  internal  fields  in  the  Si02  layer  due  to  the  trapped  negative  charges  as  a potential 
barrier  against  this  electron  injection.  This  situation  is  preserved  throughout  the  entire 
detrapping  experiment,  since  only  small  portions  of  the  total  charge  are  removed  at  each  step. 

PHOTODETRAPPING  ANALYSIS 

If  first  order  kinetics  (neglecting  charge  retrapping)  is  assumed,  the  local  depopulation  of 
occupied  traps  under  illumination  is  governed  by  the  following  equation: 

^n^  lx,E.t)/8t  - -Fp(x,hw)  <7p{x.E.hu)  n,{x,E,t),  (1) 

where  n,(x,E,t)  is  the  trapped  electron  concentration  per  unit  energy,  Fp(x,hu>)  is  the  local 
photon  flux  in  the  Si02  layer  and  Op(x.E,hw)  is  the  trap  photoionization  cross  section.  The 
photon  flux  Fp  is  a function  of  both  the  photon  energy  hw  and  the  position  x in  the  SiOi  layer 
due  to  the  optical  interference  phenomenon  [6].  This  interference  gives  rise  to  a standing 
wave  pattern  in  the  Si02  layer,  which  depends  upon  photon  energy. 

The  pbotoionization  cross  section  Op  is  a function  of  the  trap  energy  level  E in  the  Si02 
band  gap  (E^O  and  Eg  at  the  top  of  the  Si02  valence  band  and  the  bottom  of  the  conduction 
band,  respectively)  and  of  the  photon  energy  hu  since  it  includes  the  transition  probability  to  a 
final  state  E-f  hw.  and  may  be  position  dependent  through  variation  of  the  electric  field  due  to 
the  presence  of  charge  in  the  Si02  layer.  The  field  dependence  of  <7^  was  experimentally 
found  to  be  weak  an  ' therefore,  neglected  here. 

It  was  concluded  from  discharging  experiments  [4],  that  the  charge  centroid  x (measured 
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from  the  Al-Si02  interface)  is  constant  in  time  and  that,  if  an  energy  spectrum  of  trapping 
centers  is  present,  they  all  have  the  same  spatial  charge  distribution  centroid.  These  two 
observations  allow  us  to  separate  the  variables  determining  n,,  i.e. 

n,(x,E,t)  » n,°(t)  n,'(E)  n,2(x),  (2) 

where  n,'  and  n,^  are  normalized  distributions.  For  the  spatial  distribution  n(^(x)  a Gaussian 
is  used  with  the  charge  centroid  as  median  value  and  the  same  standard  deviation  as  for  the 
ion  distribution. 

Equation  1 is  integrated  over  the  oxide  thickness  L and  over  the  Si02  energy  gap,  using 
the  expression  for  the  flat-band  voltage  shift: 

AVpB(t)  - q(x/e)^  * n,(x,E,t)dEdx,  (3) 

where  e is  the  static  permittivity  of  Si02  and  q is  the  electron  charge.  At  a particular  photon 
energy  (hw),  only  centers  between  EvEg-h<j  and  E^Eg  can  be  depopulated,  if  thermal 
broadening  of  the  trapping  levels  and  two-photon  processes  are  disregarded.  We  thus  get: 


d(AVFB(t)l/dt  - 


-AVFB(t)£  Fp(x,h«)n,2(x)dxJ^  *^0j,(E,h«)n,'(E)dE. 

(4) 

r*- 

Let 

X (h«)  * Fp  (x,hw)  n,2  (x)  dx 

(5a) 

and  2 (hw)  - 1 o-(E,hw)  n,'  (E)  dE, 

•'E.-hu 

(5b) 

where  >>  the  convolution  of  the  photon  flux  with  the  spatial  distribution  of  the  trapping 

centers  and  2(h<d)  is  the  convolution  of  the  photoionization  cross  section  with  the  optically 
accessible  trap  distribution  over  energy.  The  latter  quantity  will  be  viewed  as  an  effective 
photoionization  cross  section.  If  a single  monoenergetic  trap  is  involved,  then  Zfhw)  reduces 
to  0p(hw). 


r- 
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The  solution  to  eq.  4 can  be  approximated  as: 


[AVFB(0)-AVpB(t)]/AVFB(0)  = x(hu)2(hw)t  (6) 


if 


t «[x(h<o)2(h<u)]->. 


(7) 


If  long  discharging  times  were  used,  the  charge  centroid  would  ultimately  be  determined 
by  the  minimum  in  the  standing  wave  pattern  of  the  light.  However,  since  the  illumination 
interval  t was  kept  small  compared  to  the  discharging  time  constant,  we  mdy  disregard  the 
effect  of  light  interference  upon  the  charge  centroid. 


Since  the  variation  in  flat-band  voltage  shift  due  to  illumination  for  a period  t can  be 
measured  and  the  quantity  x(h>^)  can  be  calculated  knowing  the  sample  geometry,  the  optical 
constants  of  the  various  materials  of  the  multilayer  structure,  the  incident  photon  flux  and  the 
light  energy  used,  the  effective  photoionization  cross  section  Sfhu)  can  be  determined  from 
Eq.  6. 


RESULTS 


The  photo  I-V  technique  which  is  a sensitive  method  of  determining  both  the  density  and 
the  centroid  of  oxide  charges  [7],  was  initially  used  to  ascertain  that  the  negative  oxide  charge 
is  removed  by  exposure  to  light,  and  not  compensated  by  positive  charge  [4].  It  was  also 
experimentally  verified  that  the  discharging  phenomenon  obeys  first  order  kinetics,  and  that 
the  discharging  time  constant  is  long  compared  to  the  illumination  time  [4].  Figures  1 and  2 
display  the  effective  photoionization  cross  section  spectra,  respectively  for  P and  As  implanted 
samples,  for  which  the  experimental  parameters  are  summarized  in  Table  1.  All  the  samples 
were  implanted  with  a fluence  of  1 x 10'^  cm'^,  except  for  wafer  P8C,  which  was  implanted 
with  3 X 10‘^  P*/cm^.  In  the  latter  case,  partial  penetration  of  the  P ions  into  the  Si  substrate 
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was  taken  into  account.  In  view  of  the  low  doping  level  of  the  Si02  introduced  by  the  ion 
implantation,  we  used  the  optical  constants  for  unimplanted  Si02  in  the  calculation  of  the 
standing  wave  pattern.  Incorporating  the  finite  band-pass  of  the  monochromator  (SO  A full 
width  at  half  maximum)  into  the  analysis  reduces  the  magnitude  of  Ilfhu)  in  Figs.  1 and  2 by 
25-40  %,  the  larger  correction  being  for  higher  photon  energies,  without  changing  the 
observed  structure  in  the  spectra. 

DISCUSSION 

The  nature  of  the  traps  that  exist  in  the  ion-implanted  samples  cannot  be  determined  in  a 
direct  way.  Because  of  the  high  annealing  temperature  used  after  implantation,  it  is  very  likely 
that  all  major  structural  damage  in  the  Si02  is  healed.  Therefore,  we  consider  trapping  levels 
related  to  the  implanted  ions  themselves  rather  than  due  to  the  damage. 

At  first,  one  might  expect  that  P would  occupy  Si  sites  and  form  phosphate-type  clusters, 
as  is  the  case  when  P is  diffused-in  instead  of  implanted.  That  is  possible,  but  P and  As  at  Si 
sites  would  act  as  donors  and  not  as  electron  traps  and  would,  therefore,  go  undetected  in  our 
experiments.  On  the  other  hand,  a simple  tight-binding  model  for  substitutional  impurities  in 
Si02,  based  on  the  bond-orbital  description  of  the  bulk  material  [8],  reveals  that  P and  As  are 
very  likely  to  be  substitutional  at  O sites.  An  energy  level  diagram  describing  P at  an  O site  is 
shown  in  Fig.  3.  The  energy  separation  between  the  O 2p  and  P 3p  levels  is  about  S.8  eV  [9]. 
The  P 3p  orbitals  bond  with  the  neighboring  Si  hybrid  orbitals  and,  as  a result,  two  levels  are 
expected  at  midgap.  The  tower  level  contains  two  electrons,  whereas  the  upper  level  contains 
only  one,  and  would,  therefore,  act  as  an  electron  trap.  Photoionization  of  this  electron  would 
have  a threshold  at  approximately  3.7  eV,  which  is  close  to  that  observed.  Above  threshold,  a 
resonance  corresponding  to  the  atomic  P 3p  -»  4s  transition  is  expected,  in  analogy  with  the 
ootical  soectnim  of  the  pure  material  [20].  Such  a resonance  may  be  present  in  the  data  of 
Fig.  I.  According  to  the  same  model.  As  at  an  O site  would  behave  just  like  P with  the  only 
exception  that  the  level  in  the  gap  will  be  approximately  0.4  eV  lower  in  energy.  The  data  of 
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Fig.  2 show  that  the  As  threshold  is  0.2-0.3  eV  higher  than  that  of  P,  in  agreement  with  the 
theoretical  prediction.  We  feel  that  these  results  suggest  that  the  observed  traps  are  P and  As 
at  O sites,  but  we  cannot  at  present  exclude  other  possibilities. 

If  our  model  is  correct,  we  expect  the  photoionization  cross  section  above  threshold  to  be 
insensitive  to  implantation  energy  and  oxide  thickness.  As  observed  in  Figs.  1 and  2,  the 
analysis  indeed  compensates  for  implantation  energy,  but  not  completely  for  oxide  thickness. 
Our  experimental  oxide  thickness  range  allowed  us  to  perform  accurate  flat-band . voltage 
measurements  on  MOS  devices,  but  on  the  other  hand  was  shown  to  introduce  a strong 
sensitivity  due  to  minute  thickness  variations  in  the  interference  calculation.  In  view  of  this, 
the  structure  in  the  spectra  has  to  be  considered  preliminary.  Experiments  are  in  progress  with 
thinner  oxide  layers  (200-300  A),  in  which  case  the  influence  of  oxide  thickness  variations 
upon  the  photoionization  cross  section  spectrum  should  be  minimized. 
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Table  I.  Experimental  Parameters  for  the  Samples. 


Sample 

Ion 

Oxide 

Aluminum 

Ion 

Charge 

Ion  standard 

thickness 

thickness 

energy 

centroid 

deviation 

(A) 

(A) 

(keV) 

(A) 

(A) 

PE 

P 

1283 

135 

30 

420 

142 

PG 

P 

1277 

140 

40 

545 

181 

P8C 

P 

583 

106 

25 

290  * 

121 

AS7D 

As 

1433 

120 

60 

470 

118 

AS7G 

As 

1415 

120 

80 

570 

149 

AS4X 

As 

1272 

150 

60 

465 

118 

AS8D 

As 

559 

134 

20 

185 

51 
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PHOTON  ENERGY  (eV) 


Fig.  1.  Effective  photoionization  cross  section  spectra  for  P implanted  Si02. 
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Fig.  2.  Effective  photoionization  cross  section  spectra  for  As  implanted  Si02. 
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Fig.  3.  Bond-orbitai  energy-level  structure  for  SiOt:?.  The  gap  marked  ‘~8  eV  cone- 


sponds  to  the  onset  of  absorption  (not  the  band  gap),  so  that  3 7 eV  is  the  predicted 


onset  of  absorption  of  the  impurity 
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It  was  found  that  models  could  be  constructed  with  or  without  an  SiO^  layer 
between  the  crystalline  Si  and  the  amorphous  stoichiometric  SiO,.  The  atomic 
coordinates  were  relaxed  using  a computer  program  and  a simple  force  model. 
In  order  to  probe  the  width  w of  the  SiO^  layer,  additional  O atoms  were 
gradually  inserted,  thus  reducing  w,  and  a normalized  elastic  energy  was 
monitored.  Using  a wide  range  of  choices  for  the  force  constants,  it  was  found 
that  the  elastic  energy  was  always  lowered  as  w - 0. 
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The  stoichiometry  of  the  Si-SiOj  interface  is  a subject  that  has  attracted  considerable 
attention.  The  interface  has  been  probed  experimentally  by  a wide  variety  of 
techniques,  many  of  which  are  discussed  in  a series  of  papers  in  these  Proceedings. 
The  main  objective  of  such  experiments  is  to  determine  the  width  of  the  SiO^ 
(x  9*  2)  layer  that  may  exist  at  the  interface.  A variety  of  results  has  been  report- 
ed, ranging  from  less  than  4 A (Refs.  1 and  2)  to  8 A (Ref.  3),  12-15  A (Ref.  4), 
lS-20  A (Ref.  5),  and  so  on.  Much  of  this  variation  is  due  to  the  differences  in 
experimental  techniques  and  interpretation  of  results.  For  example,  some  of  the 
techniques  may  alter  the  composition  of  the  sample  during  measurement.  Also,  some 
techniques  probe  the  scattering  properties  of  atoms,  whereas  others  probe  the 
behavior  of  electrons,  which  may  see  a different  effective  width.  This  observation 
brings  up  the  subject  of  how  to  define  an  appropriate  width.  One  possibility  is  to 
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define  the  width  strictly  on  the  basis  of  stoichiometry.  Even  then  a number  of 
alternatives  exist.  Take,  for  example,  an  abrupt  (100)  interface.  One  would 
normally  assume  that  such  an  interface  has  no  SiO^  layer.  However,  as  Stern'*  has 
pointed  out,  the  Si  atoms  in  the  last  crystallographic  plane  of  atoms  on  the  Si  side 
are  bonded  to  two  Si  atoms  and  to  two  O atoms,  corresponding  to  SiO  stoichiometry. 
[For  a (111)  interface  this  plane  of  atoms  corresponds  to  SiO^j  stoichiometry.] 
Since  it  is  only  one  plane  of  atoms,  one  might  define  the  width  to  be  the  distance 
between  atoms  on  either  side  of  that  plane,  which  is  about  3 A. 

Another  possibility  is  to  define  the  width  having  in  mind  the  behavior  of  electrons. 
Qearly,  even  if  one  had  an  abrupt  interface,  electrons  in  the  vicinity  of  the  first  few 
layers,  including  core  electrons,  would  not  behave  as  bulk  electrons  do.  One  meas- 
ure of  the  width  could  be  taken  to  be  the  range  over  which  the  self-consistent 
potential  is  different  from  bulk  values.  Calculations  by  Baraff,  Appelbaum,  and 
Hamann  have  found  this  width  to  be  very  small  for  the  Ge-GaAs  interface.  No 
calculations  are  available  for  the  Si-SiO^  interface.  Other  properties,  such  as  the 
dielectric  constant,  which  depend  on  excited  states,  may  in  fact  exhibit  differences 
from  bulk  values  over  an  even  wider  range. 

In  this  paper,  we  attempt  to  shed  some  light  on  the  question  of  the  interface  stoi- 
chiometry by  constructing  continuous-random-network  (CRN)  models.  This  techni- 
que has  been  widely  used  in  modeling  the  atomic  structure  of  amorphous  solids.  In 
fact,  pioneering  work  using  CRN  models  was  done  for  amorphous  SiOj  by  Evans  and 
King’  and  by  Bell  and  Oean.^  Later,  the  technique  was  used  to  study  amorphous  Si, 
Ge  and  other  materials. 

In  carrying  out  the  present  study  we  first  built  several  mo(^els  of  amorphous  SiO,. 
We  used  the  "Framework  Molecular  Models"  manufactured  by  Prentice-Hall,  Inc.,  of 
Englewood  Cliffs,  New  Jersey.  They  consist  of  sturdy  plastic  tubing  which  can  easily 
be  cut  to  size  and  which  fits  snugly  over  the  prongs  of  small  metal  units,  representing 
atoms.  The  units  used  to  represent  Si  atoms  have  four  metal  prongs  pointing  in 
tetrahedral  directions.  Oxygen  atoms,  which  ought  to  have  two  prongs,  were 
represented  by  small  pieces  of  bare  copper  wire  bent  at  an  angle  of  144°,  which 
corresponds  to  the  observed  average  angle  for  amorphous  SiO,.'^  Since,  however, 
the  plastic  tubing  was  slightly  bendable,  the  resulting  Si-O-Si  angles  actually  have  a 
distribution  about  144°.  The  Si-O  bonds  were  represented  by  blue  tubing  cut  to  1.6 
inches,  so  that  one  inch  corresponds  to  1 a.  Every  Si  atom  was  always  connected  to 
four  O atoms  in  the  tetrahedral  directions,  and  every  O atom  was  always  connected 
to  two  Si  atoms  in  order  to  maintain  proper  stoichiometry.  No  excessive  bending  of 
the  plastic  tubing  connecting  atoms  was  allowed.  O atoms  were  not  allowed  to 
approach  each  other  by  less  than  about  2 inches  (O  nearest-neighbor  atoms  in 
crystalline  SiO,  are  separated  by  about  2.5  A).  While  in  crystalline  SiO,  one  finds 
only  six-fold  rings  (counting  an  Si-O-Si  unit  as  a "fold"),  we  allowed  from  four-fold 
to  eight-fold  rings.^ 

For  the  purpose  of  modeling  the  interface,  we  first  had  to  choose  an  orientation  for 
the  crystalline  Si  substrate.  We  chose  the  (100)  orientation  because  it  is  the  one 
used  most  often  in  devices.  We  therefore  constructed  a substrate  of  several  (100) 
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planes  of  crystalline  Si.  Red  tubing  cut  to  2.35  inches  was  used  for  the  Si-Si  bonds 
(1  inch  » 1 A).  The  dimensions  of  the  exposed  surface  were  approximately  2 ft  x 3 
ft.  The  objective  of  the  study  was  to  continue  this  structure  and  build  several  layers 
of  amorphous  Si02. 


Fig.  1 A photograph  of  the  Si-SiOj  CRN  Fig.  2 A histogram  of  the 

model  showing  a stoichiometrically  abrupt  Si-O-Si  angle  distribution, 

interface. 


At  the  first  stage  of  our  interface  studies  we  relied  totally  on  visual  tests.  Clearly, 
we  could  use  both  Si-Si  bonds  (red)  and  Si-O  bonds  (blue)  and  construct  an  SiO^ 
layer  of  arbitrary  thickness.  We  could  also  allow  an  arbitrary  amount  of 
roughness.^- '•*  Our  objective  was  therefore  to  attempt  to  construct  a model  with  an 
interface  region  which  is  as  smooth  as  possible  and  as  thin  as  possible,  without 
excessive  distortions  of  the  bonds  and  without  dangling  bonds.  Such  an  exercise 
would  thus  not  answer  the  question  of  how  thick  the  interface  region  actually  is  in 
real  systems  (it  probably  varies  from  sample  to  sample),  but  it  would  provide  an 
indication  as  to  how  thin  it  can  be  and  perhaps  set  a lower  limit.  Our  initial  work 
indicated  that  an  SiO^  region  of  about  4-6  A was  necessary  to  avoid  excessive 
distortions.  We  finally  succeeded,  however,  in  constructing  a model  with  a complete- 
ly abrupt  interface,  i.e.,  without  an  intermediate  layer  containing  both  red  and  blue 
bonds.  The  distortions  of  the  bonds  were  similar  to  those  of  the  other  models.  A 
photograph  of  this  model  is  shown  in  Fig.  1.  The  conclusion  that  one  reaches  is  that 
an  abrupt  interface  in  real  samples  is  not  ruled  out  by  geometrical  considerations  and 
should  be  achievable,  unless  energetic  considerations  favor  an  interface  with  an  SiO^ 
layer  several  Angstroms  thick. 

In  the  second  stage  of  our  work  we  attempted  to  calculate  which  configuration  would 
actually  be  favored  by  energetic  considerations.  For  this  purpose  we  used  one  of  the 
CRN  models  which  had  an  SiO^  interface  region  of  about  6 A.  The  coordinates  of 
all  the  atoms  were  measured  with  a ruler  to  an  accuracy  of  better  than  0.1  inches. 
These  coordinates  were  then  used  as  a starting  point  in  a computer  rela.xation 
program  similar  to  the  one  used  in  Ref.  12.  The  total  force  on  each  atom  was 
calculated  as  the  sum  of  three  types  of  forces: 
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(a)  a bond-stretching  force  with  a potential 

U, --to  2 [(R  - Ri)2  - 

^ i 

where  R is  the  position  vector  of  the  atom  of  interest,  R|  is  the  position  vector  of  a 
nearest  neighbor.  R^  is  the  measured  crystalline  bond  length,  and  a is  the  bond- 
stretching force  constant. 

(b)  a bond-bending  force  with  a potential 

Ub  - Y ^ [(R  - Rj)  • (R  “ Rk)  - c]' 

where  Rj  and  R^  are  the  position  vectors  of  the  two  nearest  neighbors,  and  C is  a 
constant  such  that  the  potential  is  minimized  for  a bond  angle  of  144”,  if  the  atom  at 
R is  an  oxygen,  and  109.6”  if  the  atom  at  R is  a silicon.  )3  is  the  bond-bending  force 
constant  and  is  different  for  Si  and  O atoms  (see  below). 

(c)  a van  der  Waal’s  force  with  a Lennard-Jones  potential  with  a cutoff: 


A ^ B 
I R-Rj  1 ^ I R-Rj  I '2 


R-RJ 

< R, 

R-RJ 

>R. 

The  main  role  of  this  force  is  to  keep  the  O atoms  from  going  too  close  to  one 
another.  The  above  forces  are  certainly  not  the  totality  of  forces  on  atoms  in  solids. 
They  are,  however,  a minimal  set  needed  for  a stable  solid  and  have  been  found  to 
work  well  within  the  context  of  CRN  models. The  choice  of  force  constants  will 
be  discussed  later. 


In  the  relaxation  procedure,  each  atom  was  displaced  in  turn  toward  its  equilibrium 
position  determined  by  the  total  potential  U - U,  -f-  + U^,.  This  was  done  by 

computing  the  local  gradient  vector  and  moving  the  atom  in  the  direction  of  the  force 
by  small  steps.  The  cycle  was  repeated  until  the  atoms  approached  a good  approxi- 
mation to  equilibrium  in  both  energy  and  position. 

Once  the  coordinates  of  the  original  model  were  relaxed,  a search  was  made  by  the 
computer  to  locate  all  the  Si-Si  bonds  lying  between  6 and  5 \ from  the  top  layer  of 
the  crystalline  Si  substrate.  Oxygen  atoms  were  then  inserted  at  the  midpoints  of  all 
such  bonds  and  the  relaxation  program  was  run  again.  Fig.  2 shows  a typical 
histogram  for  the  Si-O-Si  angle  after  relaxation.  Once  equilibrium  was  reached,  the 
final  elastic  energy  was  recorded.  Si-Si  bonds  lying  between  S and  4 k were  then 
eliminated  by  inserting  oxygen  atoms  and  relaxing.  The  procedure  was  repeated  for 
each  successive  inch  until  all  Si-Si  bonds  above  the  top  layer  of  the  substrate  were 
eliminated.  The  result  for  our  initial  choice  of  force  constants  indicated  that  the 
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elastic  strain  energy  decreased  as  the  width  of  the  SiO^  region  was  decreased.  The 
same  result  was  obtained  for  the  elastic  distortion  energy  of  the  entire  system  and 
also  for  a bond-normalized  elastic  distortion  energy.'^  We  therefore  repeated  the 
entire  procedure  for  a wide  range  of  force  constants  withtn  physically  acceptable 
limits.  In  all  cases  an  abrupt  interface  was  preferred. 

The  final  conclusion  of  this  work  is  that  an  abrupt  interface  is  not  ruled  out  by  either 
geometric  or  energetic  considerations  at  T = 0”  K.  The  model  we  used  is  not 
reliable  enough  to  allow  us  to  conclude  that  the  abrupt  interface  would  in  fact  be 
formed  under  realistic  growing  conditions  at  finite  temperatures.  The  effect  of 
impurities,  defects,  such  as  dangling  bonds,  finite  temperature,  and  other  factors 
would  have  to  be  included  before  a more  definite  conclusion  can  be  reached. 
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ABSTRACT 

We  report  the  development  of  a method  for  calculating  the  electronic  structure  o(  semiconduc- 
tor surfaces.  It  is  based  on  the  Koster-Slater  idea  for  treating  localized  perturbations,  which 
was  later  extended  to  describe  surfaces.  The  present  method  makes  use  of  semi-empirical 
tight-binding  Hamiltonians  and  of  a novel  way  to  create  free  surfaces.  .A  Green's  I unction 
scattering-theoretic  formulation  is  employed.  The  properties  of  the  bulk  crystal  are  built  in 
and  preserved.  Bound-surface-state  energies  are  determined  unambiguously  and  accurately 
even  for  states  whose  wavefunctions  are  very  extended.  The  total  and  local  changes  occurring 
in  the  density  of  states  due  to  the  surface  can  be  calculated  directly  and  therefore  very 
accurately,  without  having  to  subtract  two  large  quantities.  Some  of  the  structure  in  the 
state-density  changes  is  in  the  form  of  narrow  peaks,  which  can  be  identified  as  resonances  or 
antiresonances.  In  order  to  point  out  advantages  of  the  method  and  to  compare  our  results 
with  the  results  of  slab  calculations,  we  present  applications  to  the  Si  and  Ge  (100)  free 
surfaces.  The  present  method  is  shown  to  be  very  efficient,  accurate  and  fast.  Despite  the 
fact,  that  a tight-binding  Hamiltonian  of  a truly  semi-infinite  system  is  treated  exactly,  the 
method  employs  matrices  which  are  much  smaller  than  those  arising  in  slab  calculations. 
Finally,  the  method  is  applied  to  a study  of  the  (100)  surfaces  of  the  isoelectronic  series 
Ge-GaAs-ZnSe  and  to  the  (100)  surface  of  cubic  SiO;. 
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I.  INTRODUCTION 

The  problem  of  the  electronic  structure  of  solid  surfaces  has  attracted  considerable 
attention  during  the  last  several  years. A variety  of  methods  have  been  developed  and  used 
to  study  surfaces  of  metals  and  semiconductors.  These  methods  belong  to  one  of  four  general 
categories:  (a)  methods  that  simulate  the  semi-infinite  solid  with  a finite  number  of  atoms^'*^ 
(cluster  methods);  (b)  methods  that  simulate  the  semi-inifinte  solid  with  a finite  number  of 
atomic  layers  having  two-dimensional  periodicity**'^**  (slab  methods);  (c)  methods  that  solve 
the  Hamiltonian  of  the  semi-infinite  crystal  directly  (wavefunction-matching  methbd^^'^^; 
Green's  function  matching  method^^'^*;  transfer-matrix  or  continued  fraction  method^*"*^); 
and  (d)  methods  that  treat  the  creation  of  a surface  as  a perturbation  of  an  otherwise  perfect 
infinite  solid.*^'^^  All  the  above  methods  may  be  used  with  a variety  of  one-electron  Hamilto- 
nians (semi-empirical,  Hartree-Fock,  self-consistent  pseudopotential,  etc.)  together  with  one  of 
many  possible  basis  sets  in  terms  of  which  the  wavef unctions  are  expanded  (plane  waves, 
atomic  orbitals,  muffin-tin  orbitals,  etc.) 

All  four  classes  of  methods  have  been  used  widely  in  studies  of  ideal  model  systems  as 
well  as  studies  of  real  materials.  For  semiconductors,  calculations  have  thus  far  been  done  in 
terms  of  only  the  first  three  classes  of  methods.  Methods  belonging  to  the  fourth  class,  i.e. 
those  treating  the  creation  of  a surface  as  a perturbation  of  the  otherwise  infinte  perfect 
crystal,  have  been  formulated  by  several  authors,  in  all  cases  making  use  of  Green’s  functions. 
Applications  have  been  mainly  on  ideal  systems  that  simulate  no  real  material  and  on  the 
so-called  jellium  model  used  to  simulate  simple  metals.^^-^^  More  recently,  limited  applications 
of  such  methods  have  been  reported  for  surface  states  of  some  transition  metals. 

In  this  paper  we  develop  a method  which  treats  the  creation  of  a surface  as  a perturbation 
and  IS  particularly  suited  for  semiconductors  and  insulators.  Its  foundations  are  the  same  as 
those  of  the  methods  used  or  discussed  in  Refs.  (43)-(S3),  in  that  it  is  based  on  a Green’s 
function  formulation  of  scattering  theory.  It  therefore  has  a number  of  distinct  advantages.  It 
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deals  with  a truly  semi-infinite  solid,  instead  of  a finite  cluster  of  atoms  or  finite  number  of 
layers.  As  a consequence,  the  band  continua  are  described  in  terms  of  continuous  functions, 
instead  of  being  simulated  by  a set  of  discrete  states.  Furthermore,  the  changes  produced  in 
the  electronic  structure  of  the  infinite  solid  by  the  creation  of  the  free  surface  are  obtained 
directly,  thereby  avoiding  the  subtraction  of  two  large  quantities.  From  the  computational 
point  of  view,  it  deals  with  matrices  whose  size  is  determined  by  the  range  of  the  potential 
change  due  to  the  surface.  In  contrast,  cluster  and  slab  methods  deal  with  matrices  whose  size 
is  determined  by  the  range  of  surface-state  wavef unctions,  which  is  almost  always  considerably 
larger. 

' In  our  particular  formulation,  a semi-empirical  linear-combination-of-atomic-orbitals 

(LCAO)  representation  is  used  as  in  Refs.  (43-45,  47-53)  but  a novel  way  of  creating  the 
surface  is  introduced,  which  allows  the  surface  perturbation  to  be  represented  by  very  small 
matrices.  More  specifically,  the  matrices  encountered  in  our  method  are  usually  an  order  of 
magnitude  smaller  than  corresponding  slab  matrices. 

I The  plan  of  this  paper  is  as  follows.  In  Sec.  II  we  give  the  general  results  of  scattering 

theory  in  solids,  as  it  has  been  developed  by  Callaway,^^  Garcia-Moliner^^  and  others,  and 
give  the  general  form  of  the  equations  in  an  LCAO  representation.  The  equations  are  then 
specialized  for  perturbations  which,  like  surfaces,  have  two-dimensional  periodicity.  Finally, 
the  form  of  the  perturbation  which  corresponds  to  the  creation  of  a surface  is  described  and 
the  method  of  calculation  is  discussed.  In  Sec.  Ill  we  present  applications  of  the  method  to  the 
I Si  and  Ge  (100)  surfaces  using  the  same  Hamiltonian  employed  by  Pandey  and  Phillips'^*^*  in 

I slab  calculations  in  order  to  make  a direct  comparison  and  demonstrate  the  power  of  the 

I present  method  to  describe  real  materials.  In  Sec.  IV  we  apply  the  method  to  a study  of  the 

(100)  surfaces  of  the  isoelectronic  series  Ge-GaAs-ZnSe  and  of  Si02.  We  end  with  the 
conclusions  in  Sec.  V.  Details  of  the  calculations  are  given  in  a series  of  Appendices. 
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II.  SCATTERING-THEORETIC  FORMALISM 


A.  General  Results 


Let  H*  be  a one-electron  Hamiltonian  describing  an  infinite  periodic  solid.  The  corre- 
sponding eigenvalue  problem. 


^ '^'nk  “ ^nk  V'ak 


yields  the  energy  bands  E°|^  and  the  corresponding  Bloch  functions  The  one-particle 

Green's  function  G"  of  this  system  is  defined  for  outgoing  waves  by 


G (E)  - lim  ! ^ (2) 

€-0+  E+i«-H  E+-H 

where  the  superscript  + on  E in  the  second  form  is  a short-hand  notation  for  the  full  expres- 
sion involving  the  limit  as  e->^0  from  above.  The  operator  G"  can  also  be  conveniently 

• o 

expressed  in  terms  of  the  complete  set  of  Bloch  eigensolutions  of  H as 

G“(E).  2 ' ‘ (3) 

nk  E*-E»k 

As  usual,  the  density  of  states  of  H°  is  given  by 


N (E)  - - — Im  Tr  G (E) 
ir 


Now  let  U be  any  perturbation,  and  let 


H - H 4-U 


be  the  new  one-electron  Hamiltonian  for  the  perturbed  system.  The  new  eigenvalue  problem 


(H  +U)  - E + . 


and  the  solutions  fall  into  two  categories:  States  with  energies  within  the  forbidden  gaps  of 
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the  spectrum  of  H , and  states  within  the  band  continua.  One  therefore  usually  seeks  to 
determine  the  positions  of  the  discrete  states  in  the  band  gaps  and  the  changes  produced  in  the 
density  of  states  within  the  band  continua. 


For  states  in  the  band  gaps  one  can  immediately  rewrite  (6)  in  the  form 


(1-G  U)  - 0 . 


If  ^ is  expanded  in  terms  of  an  orthonormal  basis  set  4 , i.e. 


(7) 


^ « 2 c-  , 


(8) 


eq.  (7)  can  be  transformed  to  a set  of  linear  algebraic  equations  for  c^ 


j iKn  - 2 G:,(E)  Cp  - 0 . 

where  the  Green’s  function  in  the  {4  } representation  is  directly  obtained  from  eq.  (3): 


G.^(E)  - 2 

nk 


<«.l^nk><^nkl♦^> 


E-E 


nk 


From  eq.  (9),  bound  states  correspond  to  zeros  of  the  determinant 


D(E)  - det  I 5,^  - 2 G,JE)  U « | . 
y 


(9) 


(10) 


(11) 


An  important  feature  of  this  result  is  that  the  effective  order  of  the  determinant  in  eq.  (II)  is 
equal  to  the  order  of  the  matrix  U,  so  that  the  method  is  particularly  suitable  for  short-range 
perturbations.  The  scattering-theoretic  aspect  of  this  approach  is  quite  obvious  in  that,  in 
more  formal  terms,  equation  (11)  actually  determines  the  discrete  poles  of  the  scattering 
matrix.  The  process  under  consideration  here  is,  of  course,  the  scattering  of  Bloch  waves  by 
the  perturbation  potential  creating  the  surface. 

Once  the  bound-state  energy  is  determined  from  equation  (11),  the  corresponding  wave 


function  is  calculated  from  equation  (9). 
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The  change  in  the  density  of  states  AN(E)  - N(E)  - N (E)  within  the  bands,  is  obtained 
by  first  defining  the  Green’s  function  G(E)  for  the  perturbed  system  by 


and  relating  it  to  G (E)  by  Dyson’s  equation,  i.e. 


G - G + G UG  . 


The  new  density  of  states  is  given  by 


N(E)  » - -L  Im  Tr  G(E)  . 

It 


The  change  in  the  density  of  states  AN(E)  can  now  be  conveniently  written  as 


^n(E)  . -L  ^ 

V dE 


where  the  quantity  j(E)  is  the  phase  shift,  given  by 


«(E)  - - tan-'  (Im  D(E)/Re  D(E)1 


where  0(E)  is  as  defined  by  eq.  (II).  It  should  be  noted,  that  AN(E)  depends  only  on  G and 
U,  so  that  its  evaluation  does  not  require  the  solution  of  Dyson’s  equation  (13).  Equation 
(16)  shows  that  6(E)  goes  through  an  odd  multiple  of  v/2  every  time  Re  D(E)  m 0.  If  such 
an  energy  is  labeled  £„,  one  may  then  expand  D(E)  in  the  vicinity  of  E^  to  obtain  the 
corresponding  contribution  to  the  change  in  the  density  of  states  in  a Lorentzian  form  centered 


an(E)  - L 1__  . 

2ir  (E-Eo)2  + rV4 

with  the  half-width  F given  by 


2 Im  D(E„) 
Re  D'(E„) 


Page  6 


In  (18).  the  prime  denotes  differentiation  with  respect  to  energy.  When  r>0,  AN(E)  corre- 
sponds to  a resonance  of  half-width  F,  whereas  when  r<0,  AN(E)  corresponds  to  an  antireso- 
nance with  half-width  | F | . Notice  that  the  same  expressions  apply  within  the  gaps  where  Im 
D(E)  « 0 so  that  F >>0  and  eq.  (17)  reduces  to  a 5-function  indicating  a discrete  state. 

AN(E)  satisfies  the  important  sum  rule 

J AN(E)  -0.  (19) 

00 

known  as  Levinson’s  theorem.  When  the  discrete  states  in  the  gaps  are  counted  separately, 
(19)  becomes 

f ilN(E)  - - Nb  (20) 

bands 

where  N),  is  the  number  of  discrete  states  in  the  gaps. 

The  relation  for  the  discrete  states,  eq.  (9),  as  it  applies  to  point  defects  in  semiconduc- 
tors was  first  derived  by  Koster  and  Slater. A more  general  treatment  has  been  given  by 
Callaway^^  and  by  Garcia-Moliner.^^  Callaway  employed  Wannier  functions  for  the  set 
which  are  convenient  for  general  analytical  results  but  very  cumbersome  for  actual  calcula- 
tions. More  recently,  the  problem  was  formulated  in  terms  of  LCAO  basis  sets  in  the  above 
fashion  by  Bernholc  and  Pantelides,^  who  used  the  method  to  calculate  the  electronic 
structure  of  the  vacancy  in  several  semiconductors.  Here,  we  turn  to  surface-related  problems 
and  study  an  arbitrary  perturbation  that  retains  two-dimensional  periodicity.  Before  we 
proceed  further,  however,  we  must  introduce  some  necessary  LCAO  terminology. 

B.  LCAO  Representation  ; Planar  Perturbations 

In  order  to  describe  a perfect  infinte  crystal  in  an  LCAO  representation,  we  introduce 
atomk-like  orbitals  on  each  atom  in  the  primitive  unit  cell  and  denote  them  by  Va(r-^),  where 
^ are  the  positioas  of  the  atoms  in  the  unit  cell  and  a labels  the  s,  p ...  etc.  character  of  the 
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orbitals.  For  the  purpose  of  solving  for  the  eigensolutions  of  H°  we  first  define  the  Bloch 
sums 

I ik»(R.+T  ) 

xrd)  - ^ 2 e-  -J  - 9..  (r-Rj-i,)  (21) 

where  Rj  are  the  Bravais  lattice  vectors  and  N3  is  the  number  of  Bravais  lattice  points  in  the 
3>dim.  bulk  lattice.  The  Bloch  functions  are  then  expanded  as  follows: 

•Anicd)  - 2 c;,(k)  xj'd)  • (22) 

- aff 

Upon  diagonalization  of  the  secular  matrix  <xn  I Xr  ’^  > 3t  each  k,  one  obtains  the 
eigenvalues  corresponding  eigenvectors  which  define  the  Bloch  functions.  In 

the  applications  to  be  presented  in  this  paper,  the  matrix  elements  I H°  | <p^/^.Y>  which 
are  needed  to  construct  <Xk1  | xj  *>  are  treated  as  parameters.  The  method  has  often  been 
referred  to  as  the  Empirical  Tight  Binding  Method  (ETBM). 

For  an  arbitrary  perturbation  U,  one  must  first  choose  a basis  set  in  terms  of  which  to 
expand  the  new  wavefunctions  ip,  express  the  Green’s  function  G*.  evaluate  the  matrix  of  U, 
and  proceed  with  the  formalism  described  in  the  previous  subsections.  In  early  formulations  of 
the  problem,  the  Wannier  functions  were  thought  to  be  the  natural  set  for  localized  perturba- 
tions. When  an  LCAO  representation  is  used  for  the  bulk  band  structure,  however,  the  set  of 


atomic  orbitals  has  been  shown  to  be  a more  natural  and  useful  choice  for  point 
defects.^*  Along  the  same  lines,  for  planar  perturbations,  it  is  natural  to  use  the  set  {9:),}  and 


define  layer  orbitals  in  order  to  take  advantage  of  the  fact  that  the  system  retains  two- 
dimensional  periodicity.  Because  of  this  periodicity,  a two-dimensional  3 vector  is  a good 
quantum  number  and  layer  orbitals  are  simply  Bloch  sums  in  two  dimensions.  They  are 
defined  by 


♦5'"'(r) 


^ j 


V.fr-gj-Ag^ 


where  £j  are  the  lattice  vectors  of  the  two-dimensional  Bravais  lattice  and  ^ are  the  position 
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vectors  of  the  atoms  in  the  two-dimensional  unit  cell.  The  number  of  Bravais  lattice  points  in 
the  2-dim.  lattice  is  denoted  as  N2.  The  index  m labels  different  planes  of  atoms  and  the  index 
r labels  different  atoms  in  the  same  plane.  The  above  definition  of  layer  orbitals  is  the  same 
as  that  used  for  slab  ‘calculations  in  an  LCAO  representation.  For  convenience,  we  will 
occasionally  use  a composite  index  f «(ami’)  for  the  layer  orbitals. 

States  of  the  perturbed  system  are  then  expanded  in  terms  of  the  layer  orbitals 

(r)  - 2 (r)  . (24) 

where  s numbers  the  various  states  at  each  3.  The  Green’s  function  G‘  and  the  perturbation 
matrix  U can  then  be  evaluated  in  the  representation  at  each  £,  from  which  the  determi- 
nant D(E),  the  bound  states,  the  change  in  the  state  density  AN(E),  etc.  can  be  determined. 

C Free  Surfaces 

The  next  task  is  to  identify  the  perturbation  U describing  the  creation  of  a free  surface. 
There  are  essentially  two  ways  to  accomplish  this  task.  One  possibility  is  to  "switch  off"  the 
interactions  between  orbitals  on  the  atoms  on  a number  of  adjacent  planes  (Fig.  la)  so  that 
two  non-interacting  semi-infinite  solids  are  formed.  Clearly,  the  number  of  planes  involved  in 
this  procedure  is  determined  by  the  number  of  nearest-neighbour  interactions  that  are  included 
in  the  tight-binding  Hamiltonian  H*  of  the  bulk.  For  example,  if  only  first-nearest-neighbour 
interactions  are  included  in  the  bulk  Hamiltonian,  switching  off  the  interactions  between  only 
two  adjacent  planes  will  accomplish  the  decoupling  between  the  twin  surfaces.  The  interaction 
between  orbital  a and  orbitals  j3  on  two  different  atoms  is  formally  switched  off  by  simply 
setting  « - H*d.  In  this  procedure  both  atoms  become  involved  in  the  process  of  creating 
a surface. 

An  alternative  way  to  create  a free  surface  is  to  remove  one  or  more  layers  of  atoms  so 
that  again  two  non-interacting  semi-infinite  solids  are  formed  (Fig.  Ib).  As  it  is  clear  from 
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Fig.  I,  the  number  of  atoms  involved  in  creating  the  surface  is  now  smaller  (in  most  cases 
smaller  by  a factor  of  two).  The  "removal"  of  atoms  is  accomplished  by  simply  setting  their 
diagonal  matrix  elements  to  a constant  u and  then  letting  that  constant  go  to  infinity.  This 
procedure  was  first  introduced  by  Bernholc  and  Pantelides^  in  the  case  of  the  vacancy  in  a 
bulk  semiconductor.  In  fact,  the  creation  of  a surface  is  equivalent  to  the  introduction  of  a 
"planar  vacancy"  (in  the  case  of  nearest-neighbour  coupling)  or  "planar  divacancy"  (in  the 
case  of  next-nearest-neighbour  coupling)  etc.  For  example  in  a one-dimensional  crystal  with 
only  nearest-neighbour  interactions,  a vacancy  creates  two  non-interacting  semi-infinite  solids. 
Algebraically,  the  matrix  elements  of  the  perturbation  U between  atomic  orbitals  a and  P are 
of  the  form 


(25) 

for  a and  P on  the  atoms  to  be  removed.  With  this  form  of  perturbation  and  in  the  limit  u-»9e, 
the  condition  for  bound  states,  eq.  (11),  becomes 

D(E)  - det  I G?/»'(E)  B - 0 (26) 

where  / and  ('  are  layer  orbitals  on  the  planes  of  atoms  to  be  removed.  The  change  in  the 
density  of  states  AN(E)  is  given  by  the  same  expression  (IS)  and  (16)  with  D(E)  replaced  by 
D(E)  defined  by  eq.  (26). 

The  new  procedure  of  creating  free  surfaces,  which  we  will  refer  to  as  the  "removal 
method",  has  a number  of  distinct  advantages  over  the  previously  used  method,  which  we  will 
refer  to  as  the  "switch-off  method",  (a)  In  the  "switch-off"  method  one  must  set  up  the 
Green’s  function  matrix  and  the  nondiagonal  perturbation  matrix  whereas  in  the 

"removal  method"  it  is  sufficient  to  calculate  the  Greens’  function  matrix,  as  can  be  seen  in 
equation  (26).  This  result  illustrates  that  the  bulk  Green’s  function  contains  complete 
information  about  the  electronic  structure  of  an  ideal  surface  (i.e.  no  change  of  the  bulk  matrix 
elements  up  to  the  surface),  (b)  Computationally,  the  dimension  of  the  resultant  matrix  in  the 


"removal  method"  is  usually  smaller  by  a factor  two.  (c)  The  Green’s  function  matrix  G°  is 
Hermitian  for  the  bound  states,  whereas  1-G°U,  encountered  in  the  "switch  off"  method  is 
generally  not.  As  a result,  the  labor  may  be  reduced  at  various  stages  by  at  least  a factor  of 
four  or  eight. 

In  the  above  discussion  we  compared  the  relative  merits  of  two  different  ways  of  con- 
structing the  perturbation  that  creates  a free  surface  in  the  scattering-theoretic  formalism.  In 
either  way,  the  relevant  matrices  are  substantially  smaller  than  the  secular  matrices  encoun- 
tered in  slab  calculations.  In  particular,  as  will  be  seen  in  the  next  section,  the  secular  matrices 
of  typical  slab  calculations  are  an  order  of  magnitude  larger  than  our  Green’s  function 
matrices. 

III.  THE  (100)  SURFACE  OF  Si  and  Ge  - COMPARISON  WITH  OTHER  CALCULA- 
TIONS 

In  this  section  we  present  calculations  of  the  electronic  structure  of  the  ideal  (100) 
surface  of  Si  and  Ge  in  order  to  illustrate  that  calculations  using  the  scattering-theoretic 
method  described  above  are  in  fact  rather  straightforward  for  real  materials  and  can  easily 
provide  a wealth  of  information. 

The  geometrical  arrangement  of  the  atoms  for  the  100  surface  of  Si  is  shown  in  Fig.  2. 
The  axes  shown  are  the  crystallographic  y and  z axis.  The  two-dimensional  unit  cell  is  marked 
out  by  the  dashed  line;  it  contains  one  atom  in  each  layer.  The  corresponding  two-dimensional 
Brillouin  zone  is  shown  in  Fig.  3.  We  will  use  the  same  ETBM  Hamiltonian  used  by  Pandey 
and  Phillips'^-^*  to  study  the  (100)  surface  of  Si  and  Ge  so  that  a direct  comparison  will  be 
possible.  This  Hamiltonian  uses  only  s and  p orbitals  on  every  atom  and  retains  only  first-  and 
some  second-nearest  neighbour  interactions.  The  values  of  the  parameters  in  standard 
notation  are  given  in  Table  I.  An  ideal  surface  is  then  defined  as  an  abrupt  termination  ol  the 
bulk  crystal,  keeping  all  the  intra-  and  inter-atomic  interactions  unaltered.  In  this  approxima- 
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tion,  no  new  parameters  are  necessary  for  the  calculation  of  the  electronic  structure  of  the  free 
surface.  The  creation  of  the  surface  can  be  accomplished  by  removing  two  layers,  as  shown  in 
Fig.  4. 

I 

A.  Energy  Levels  of  Bound  Surface  States 

In  the  scattering-theoretic  method,  the  first  task  is  the  evaluation  of  the  projected  band 
structure  (PBS),  i.e.  the  projection  of  the  bulk  bands  with  k - (3,  k^^)  for  each  g in  the 
surface  Brillouin  zone.  The  PBS  is  needed  for  the  calculation  of  the  G°-matrix  (Appendix  A), 
but  it  is  also  extremely  useful  in  its  own  right.  It  allows  one  to  identify  all  the  gaps  and 
"pocitets"  where  a search  must  be  made  in  order  to  determine  the  positions  of  states  that  are 
truly  bound  at  the  surface.  We  have  calculated  the  PBS  for  Si  (100)  and  Ge  (100)  and  display 
them  in  Figs.  5 and  6.  Instead  of  the  usual  practice  of  cross  hatching  the  continuum  regions 
uniformly,  we  display  the  actual  projected  eigenvalues  at  each  3 for  a fixed  set  of  kj^  values. 
These  are  precisely  the  subset  of  bulk  states  which  would  have  to  be  used  if  the  surface-state 
wavefunctions  were  to  be  expanded  in  terms  of  bulk  Bloch  functions.  The  distribution  of 
points  at  each  3 in  Figs.  S and  6,  therefore,  provides  visual  information  about  which  bands 
contribute  more  heavily  to  the  formation  of  surface  states.  In  particular  it  gives  a visual 
impression  of  the  density  of  bulk  states  at  each  3. 

Bound  states  in  the  gaps  and  pockets  of  the  PBS  are  determined  in  an  unambiguous  way 
by  searching  for  the  zeros  of  D^(E)  given  by 

Dq(E)  - det  II  1 G*  I <I>|  m'M'>  | (27) 

where,  for  the  examples  under  consideration,  the  indices  m and  m'  run  over  only  the  two  layers 
that  are  removed  for  the  creation  of  the  twin  non-interacting  semi-infinite  solids  (Fig.  4).  The 
size  of  the  relevant  Green’s  function  matrix  is  therefore  8x8  and  can  be  evaluated  with 
modest  effort  (Appendix  A).  The  search  for  bound  states  can  then  be  made  by  using  an 
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efficient  algorithm  described  in  Appendix  B.  The  results  are  shown  as  solid  lines  in  Figs.  5 
and  6.  We  will  collectively  refer  to  the  surface  bands  as  the  surface  band  structure  (SBS). 

We  turn  now  to  a comparison  of  the  above  results  with  those  obtained  by  Pandey^*  using 
the  slab  method  in  which  the  semi-infinite  solid  is  simulated  by  a slab  of  finite  thickness.  For 
the  (100)  surfaces  of  Si  and  Ge  Pandey^^  used  20  layers  of  atoms.  The  electronic  states  of 
the  slab  were  then  expanded  in  terms  of  layer  orbitals  and  the  corresponding  secular  matrix 

9 / / / 

<<(""''■  I H I <^|  "*  **  > was  set  up  and  diagonalized  at  each  £.  The  size  of  the  slab  matrix  is 
thus  80x80,  which  is  an  order  of  magnitude  larger  than  the  corresponding  scattering-theoretic 
matrix  of  eq.  (27).  The  slab  matrix  needs  of  course  be  diagonalized  only  once  at  each  g, 
instead  of  at  a mesh  of  energies  and  furthermore  the  Hamiltonian  matrix  elements  encountered 
in  the  slab  method  are  easier  to  evaluate  than  the  Green’s  function  matrix  elements.  Neverthe- 
less, the  scattering  theoretic  method  is  less  time-consuming  if  our  way  of  creating  the  surface  is 
used  together  with  the  efficient  root-finding  algorithm  employed  in  our  calculations.  As  for 
the  actual  results  in  Figs.  S and  6,  a comparison  with  the  slab  calculation  shows  excellent 
agreement.  For  highly  localized  states  the  surface  state  energies  in  the  case  of  Si  agree  within 
better  than  0.01  eV  so  that  there  was  no  reason  to  plot  the  slab  results  separately.  In  the  case 
of  Ge  the  same  agreement  is  found  for  most  of  the  states  except  for  a rigid  shift^^  of  0.06  eV. 
(States  in  the  small  pockets  and  the  empty  states  in  the  conduction  band  were  not  reported  in 
Ref.  S8.)  As  one  might  expect,  for  extended  states  there  are  appreciable  differences  in  the 
energy  positions  of  surface  states.  One  example  is  the  surface  state  at  T in  Ge.  We  find  that 
the  state  extends  over  more  than  20  layers.  Thus,  it  is  not  surprising,  that  the  slab-calculation 
result  deviates  from  our  surface  state-energy  by  0.32  eV.  Theoretically,  these  differences  are 
sigiuficant  in  that  we  try  to  solve  a given  model  problem  as  accurately  as  possible.  Experimen- 
tally, extended  states  are  difficult  to  detect,  so  that  the  slab  method  describes  the  experimen- 
tally interesting  states  satisfactorily.  We  will  now  turn  to  a discussion  of  surface  state 


wavef  unctions. 
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Surface-State  Wavef unctions 

The  calculation  of  the  wavefunction  in  our  method  is  rather  straightforward,  and  is 
described  in  detail  in  Appendix  C. 

In  Figs.  7 and  8 we  plot  orbital  sums  of  the  wavefunction  amplitude  on  each  layer  defined 
by 

f^fE)  - 2 B 12,  (28) 

a 

for  several  surface  states  of  Si  and  Ge.  In  Fig.  7 we  plot  for  comparison  a number  of  wave- 
functions  with  distinct  behavior.  Example  7a  shows  the  wavefunction  amplitude  of  a gap  state 
which  is  completely  localized  in  the  surface  plane.  An  orbital  decomposition  of  fj*  q(E)  shows, 
that  this  state  is  50%  Py-  and  50%  p^-  like.  This  type  of  bond  acts  like  a link  between 
surface  atoms  and  is  therefore  referred  to  as  the  bridge  bond.  The  second  type  of  gap  states 
usually  encountered  on  a 100-surface  in  group  IV  semiconductors  is  the  so  called  dangling 
bond.  The  wavefunction  corresponding  to  the  dangling-bond  state  at  J is  shown  in  Fig.  7b. 
The  orbital  decomposition  of  f^  ^^fE)  for  this  state  shows  that  it  has  predominantly  s-  and  p^- 
character,  resulting  in  charge  density  that  "dangles"  perpendicular  to  the  surface.  A typical 
example  for  a backbond  is  given  in  Fig.  7c.  This  stale  is  again  very  localized,  but  its  amplitude 
is  shared  mainly  by  the  first  two  layers  indicating  a charge  buildup  in  the  "back  bond". 
Example  7d  shows  the  wavefunction  of  an  empty  state  in  one  of  the  tiny  pockets  in  the 
conduction  band  part  of  the  PBS.  Such  states  lie  very  near  the  bulk  energy  continua  and  are 
therefore  very  extended.  Similar  states  in  the  valence  band  pockets  are  more  localized,  but 
interestingly  enough  they  are  localized  in  the  third  to  fifth  layer. 

When  a surface  state  band  approaches  the  bulk  band  projection,  the  states  might, 
depending  on  the  symmetry  of  their  wavefunctions,  become  more  and  more  extended.  This  is 
shown  for  the  dangling-bond  band  in  Ge  between  T and  K in  Fig.  8.  At  K (Fig.  8d)  the 
wavefunction  decays  within  2 layers.  As  the  surface  band  approaches  the  bulk  band  projection 
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(Figs.  8b  and  c)  the  wavefunction  becomes  more  extended  and  eventually  at  F (Fig.  8a)  it 
extends  over  more  than  20  layers. 

We  have  included  in  Figs.  7 and  8 a comparison  between  the  surface  state  energies 
calculated  by  our  method  and  by  the  slab  method  (energy  values  in  parenthesis)  for  the  states, 
whose  wave  functions  we  show.  For  highly  localized  states  (see  e.g.  7a,  7c,  8c,  8d)  the 
energies  are  exactly  the  same  in  the  two  calculations.  For  extended  states  (7b,  8a,  8b)  more 
or  less  pronounced  deviations  occur.  For  such  states  the  slab  method  would  have  to  employ 
thicker  and  thicker  slabs.  The  accuracy  of  the  scattering-theoretic  method,  remains  unaffected 
by  the  extent  of  the  wavefunctions,  because  it  deals  with  a truly  semi-infinite  solid.  Where  the 
results  of  the  two  methods  disagree,  the  scattering-theoretic  method  is  to  be  viewed  as  more 
accurate. 

C.  Changes  in  State  Densities  - Surface  Resonances 

In  addition  to  the  bound  surface  states  lying  within  the  gaps  and  pockets  of  the  PBS,  the 
presence  of  a surface  also  induces  changes  in  the  density  of  states  within  the  band  continua. 
These  changes  can  be  calculated  directly  by  using  eq.  (IS).  For  the  particular  problem  at  hand 

we  have 

AN-(E)  - — «(E)  (28) 

a ir  dE 

where 

S^(E)  m - tan- 1 [Im  Dj(E)/Re  6j(E))  (29) 

The  changes  given  by  eq.  (28),  however,  apply  to  the  complete  "algebraic"  system,  consisting 
of  two  semi-infinite  solids  and  the  "removed"  layers  with  states  at  infinity.  In  particular,  if  we 
define  n^  to  be  the  new  partial  density  of  states  for  the  m'**  layer  and  n*  to  be  the  partial 
density  of  states  for  any  layer  in  the  infinite  bulk  crystal,  we  have 


Af5.(E)  - I [nWE)-n*(E)]  . (30) 

* m ^ ^ 
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where  the  sum  is  over  ail  the  layers,  including  those  "removed"  by  the  perturbation.  For  the 
actual  physical  system  of  interest,  i.e.  a semi-infinite  solid,  the  relevant  quantity  should  not 
include  the  "removed"  layers  and  their  states  at  infinity.  Leaving  those  terms  out.  we  obtain 

:iNq(E)  - lANq(E)  + 2n"(E)l/2  (31) 

where  the  factor  2 inside  the  square  brackets  corresponds  to  the  number  of  rembved  layers  and 
the  division  by  2 gives  the  final  value  for  one  semi-infinite  solid. 

Both  quantities,  i.e.  AN^(E)  and  AN^fE)  were  calculated  and  found  to  satisfy  Levinson’s 
theorem:  AN^fE)  integrates  over  the  energy  range  of  the  bulk  bands  to  -8  corresponding  to 
the  8 states  at  infinity  introduced  by  the  "removal"  perturbation  and  AN^(E)  integrates  to 
zero.  The  number  of  states  is  thus  conserved. 

A typical  AN^^fE)  for  Si  at  q-(0.5,  0.5)  2>r/a  (the  J point  of  the  surface  Brillouin  zone 
see  Figs.  3 and  5),  broadened  by  0.1  eV  is  given  in  Fig.  9.  If  comparison  is  made  with  Fig.  5, 
one  observes  that  the  spikes  at  1.5  eV  and  -1  eV  correspond  to  true  surface  states  in  the  gap; 
the  remaining  structure  corresponds  to  changes  in  the  continuum  density  of  states.  Some  of 
the  structure  can  be  identified  as  resonances  and  antiresonances,  as  discussed  in  Sec.  If.  For 
example,  the  negative  spikes  at  -10  eV,  -1.2  eV.  and  4.3  eV  are  clearly  antiresonances.  The 
structure  that  appears  in  the  range  -8.5  eV  to  -5  eV  looks  rather  peculiar.  At  first  glance, 
one  might  view  it  as  two  neighbouring  resonances  giving  rise  to  a twin-peak  feature.  Such  an 
explanation  is  not  adequate,  however,  because  the  integral  of  ANg(E)  over  this  range  corre- 
sponds only  to  one  extra  state.  More  careful  examination  reveals  that  the  structure  corre- 
sponds to  two  resonances  of  about  the  same  width  plus  a wider  antiresonance  between  the 
two.  The  two  negative  dips  are  thus  the  remnants  of  the  antiresonance  tails.  Therefore,  the 
net  change  in  the  number  of  states  in  this  energy  region  should  be  one,  as  it  was  indeed  found 
to  be.  This  stiucture,  therefore,  helps  demonstrate  the  fact  that  the  change  in  the  continuum 
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density  of  states  may  be  considerably  more  complicated  than  a mere  succession  of  Lorentzian- 
type  resonances  and  antiresonances. 

If  we  integrate  ANjj(E)  of  Fig.  9 over  the  valence  bands,  i.e.  up  to  about  -1  eV  (the  top 
of  the  projected  valence  bands  at  J;  see  Fig.  S),  we  find  that  a total  of  one  state  has  been 
removed.  Similarly,  if  we  integrate  over  the  energy  range  of  the  projected  conduction  bands 
alone,  we  also  find  that  a total  of  one  state  has  been  removed.  The  two  missing  states  are  of 
course  balanced  by  the  two  states  in  the  gap,  and  Levinson’s  theorem  is  satisfied.  However, 
the  finding  also  reveals  that  the  two  bound  states  are  derived  from  both  the  valence  and  the 
conduction  bands.  This  realization  raises  some  questions  about  the  validity  of  a semi-empirical 
Hamiltonian  which  is  chosen  to  fit  the  valence  bands  very  well  while  doing  rather  poorly  for 
the  conduction  bands.  It  is  plausible  that  the  lower  of  the  two  gap  states  is  determined 
predominantly  by  the  valence  bands  and  is  thus  quite  reliable,  while  the  upper  state  is  deter- 
mined by  the  conduction  bands  (this  point  is  explored  further  in  Sec.  Ill  D below).  The  most 
crucial  test  of  a semi-empirical  surface  Hamiltonian  is  of  course  a comparison  with  the  surface 
states  obtained  by  self-consistent  calculations,  as  pursued  by  Pandey,^^  who  concluded  that  the 
results  were  quantitatively  accurate. 

In  Fig.  10  we  give  the  total  change  in  the  density  of  states  AN(E),  which  was  obtained  by 
summing  over  g in  the  surface  Brillouin  zone.  The  main  features  can  be  understood  by 
referring  to  Fig.  S.  The  states  in  the  gap  correspond  primarily  to  what  we  described  as  the 
bridge-bond  band.  The  peak  at  about  —0.8  eV  corresponds  to  the  dangling-bond  band,  also 
shown  in  Fig.  S.  Finally,  the  extra  states  at  about  -6  eV  correspond  to  the  surface  band  that 
runs  from  a point  near  K toward  J’  in  Fig.  S. 

We  turn  now  to  compare  with  the  slab  method  once  more.  By  its  nature,  the  slab  method 
is  less  suitable  for  the  description  of  changes  in  the  band  continua  because,  unlike  the  scatter- 
ing theoretic  method,  it  simulates  the  continuum  with  a set  of  discrete  states.  In  the  case  of 
the  (100)  surface  of  Si  and  Ge  discussed  here,  the  80x80  slab  secular  matrix  yields  80  energy 
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levels  at  each  g.  A few  of  these  correspond  to  the  bound  states  in  the  gaps  and  pockets 
discussed  earlier.  Some  of  the  remaining  states  which  lie  within  the  band  continua  may  be 
identified  as  resonances  if  their  wavefunctions  exhibit  localization.  For  example  in  Ref.  S8  a 
resonance  was  found  at  (he  J point  of  Si  (100)  at  —7  eV  which  is  in  the  middle  of  the  broad 
structure  we  discussed  earlier  (Fig.  9).  We  found  that  there  are  actually  two  resonances  and 
an  antiresonance  within  that  region.  It  is  clear  that  the  procedure  of  examining  slab  wavefunc- 
tions is  likely  to  miss  both  antiresonances  and  broad  resonances.  The  latter  correspond  to 
charge  localization  building  up  from  small  contributions  by  states  over  a rather  wide  energy 
range. 

One  procedure  that  has  been  used  in  conjunction  with  slab  calculations  is  to  construct  a 
histogram  out  of  the  discrete  states.  Such  histograms  for  each  g tend  to  be  rather  crude 
because  of  the  small  number  of  states.  They  become  more  meaningful  for  the  total  density  of 
states  when  a sum  over  g is  carried  out.  Note,  however,  that  this  procedure  yields  the  total 
density  of  states  for  the  perturbed  system.  One  must  subtract  large  similar  total  state  densities 
for  the  perturbed  and  unperturbed  systems  in  order  to  obtain  the  small  changes  due  to  the 
perturbation  - a clearly  error-prone  procedure. 

D.  Local  Densities  of  States 

In  the  previous  subsection  we  discussed  the  total  changes  in  the  density  of  states  within 
the  band  continua.  Those  curves  contain  resonances,  antiresonances  and  other  structure,  but 
convey  no  information  about  how  localized  these  features  are.  What  is  needed  is  something 
analogous  to  the  wavefunctions  of  bound  states  discussed  in  subsection  B above.  Such 
information  at  each  g vector  is  contained  in  local  densities  of  states  which  are  defined  for  each 
layer  by 

N^(q,  E)  - - -^  Im  Tr  G//(q,  E)  . 
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The  calculation  of  these  quantities  needs  the  explicit  determination  of  the  new  Green’s 
function  G in  a layer-orbital  representation.  For  this  purpose,  Dyson’s  equation  leq.  (13)1 
must  be  solved.  The  procedure  is  given  in  detail  in  Appendix  D.  Typical  results  of  such 

calculations  are  shown  in  Fig.  11,  where  9 was  chosen  to  be  (0.5,  0.5),  as  was  the  case  in  Fig. 

i 

9.  In  Fig.  11a  we  display  N/>(q,  E)  beginning  with  the  surface  layer  and  continuing  into 
I the  bulk  to  /»7,  for  which  N^(q,  E)  becomes  almost  identical  with  the  layer  density  of  states 

[ in  the  bulk  infinite  crystal.  The  curves  are  broadened  by  0.3  eV  and  contain  the  bound  states 

[ 

[ in  the  gaps  as  well.  In  Fig.  lib  we  display  the  changes  in  the  local  density  of  states  for  each 

layer  which  illustrate  more  clearly  the  decay  characteristics  of  the  bound  states,  the  resonances 
and  antiresonances.  In  particular  note  that  the  bound  state  at  1.54  eV,  whose  wavefunction 

I 

I was  shown  in  Fig.  7a,  is  once  more  seen  in  Fig.  11  to  be  localized  in  the  first  layer.  In 

contrast,  the  bound  state  at  —0.9  eV  (Fig.  7b),  is  seen  in  Fig.  11b  to  extend  over  several 
layers.  We  also  see  in  Fig.  11b  that  the  prominent  antiresonances  at  about  —10  eV  is 
localized  within  about  three  layers,  whereas  the  broad  feature  between  —6  and  —8  eV,  which 
we  discussed  in  connection  with  Fig.  9,  extends  over  more  than  seven  layers.  Note  that  for  the 
first  two  layers,  this  feature  appears  to  be  a single  broad  resonance  at  about  —7  eV,  and  was 
so  identified  in  the  slab  calculation  of  Ref.  58.  Fig.  1 1 confirms  the  interpretation  of  the 
broad  feature  as  a sum  of  two  resonances  and  one  antiresonance,  as  discussed  in  subsection  C 
above. 

In  conclusion,  we  find  that  the  scattering-theoretic  method  provides  a much  more  direct 
and  accurate  description  of  the  changes  in  the  continuum  density  of  states  than  the  procedures 
used  in  slab  methods. 

The  Koster-Slater-type  scattering-theoretic  approach  also  has  advantages  over  the 
transfer-matrix  method.^^  Both  these  methods  describe  truly  semi-infinite  solids  and  the  labour 
involved  may  be  comparable.  However  in  the  transfer-matrix  method  one  calculates  the  new 
Green’s  function  directly,  without  making  use  of  the  unperturbed-crystal  solutions.  $0  that 
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changes  in  the  electronic  structure  must  be  calculated  by  subtracting  two  large  quantities  as  in 
slab  techniques.  Furthermore  the  calculation  of  the  surface  energy  bandstructure  in  the 
transfer-matrix  method  is  significantly  more  laborious  as  compared  to  the  scattering  theoretic 
approach. 


E.  ^ Alternate  Calculation  for  Ge  ( 100) 

In  addition  to  our  calculation  based  on  the  Pandey,  Phillips  Hamiltonian  for  Ge,  we  have 
also  performed  surface  state  calculations  for  Ge  (100)  using  the  tight  binding  parameters  of 
Chadi^^.  Chadi  obtained  a good  fit  to  the  Ge  valence  bands  with  only  first-nearest-neighbour 
interactions  (Table  1).  By  comparing  our  results  of  the  two  calculations  we  are  able  to  arrive 
at  several  conclusions  about  the  origins  and  the  nature  of  the  surface  states. 

The  complete  PBS  and  the  surface  states  obtained  with  Chadi’s  parameters  are  given  in 
Fig.  12,  which  should  be  compared  with  Fig.  6.  It  is  clear  that  the  valence-band  projections 
are  very  similar  in  the  two  cases  whereas  the  conduction-band  projections  are  quite  different. 
The  first-nearest-neighbour  parameters  yield  a conduction  band  with  a width  of  7 eV  and  its 
minimum  at  P.  The  second-nearest-neighbour  Hamiltonian  yields  a conduction  band  with  a 
width  of  3.5  eV  and  its  minimum  at  the  L-point  (which  maps  onto  J and  I*).  It  is  also  clear 
from  Figs.  6 and  12  that  the  surface  states  obtained  in  the  two  cases  are  qualitatively  similar. 
In  the  valence  bands  there  is  even  quantitative  agreement  in  both  the  energy  position  and 
dispersion.  The  dangling-bond  band  calculated  with  the  Chadi  Hamiltonian  is  shifted  upward 
by  only  0.3  eV,  whereas  the  bridge-bond  band  is  shifted  upward  by  about  1.3  eV  with  respect 
to  Fig.  6.  The  width  of  the  bridge-bond  band  calculated  with  the  first-nearest-neighbour 
Hamiltonian  is  smaller  (.3  eV)  than  the  width  resulting  from  the  second-nearest  neighbour 
Hamiltonian  (1.45  eV).  These  results  suggest  that  the  dangling-bond  band  is  indeed  derived 
mainly  from  the  valence  bands,  whereas  the  bridge-bond  band  is  derived  mainly  from  the 
conduction  bands.  It  should  be  noted,  that  in  the  first-nearest-neighbour  approximation  the 
bridge  bond  band  lies  at  the  atomic  p-level  energy.  This  is  a direct  consequence  of  the  high 
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localization  of  the  bridge  bond  and  the  fact,  that  neighbouring  atoms  in  the  (100)  plane  are 
second  nearest  neighbours. 

1 

IV.  APPLICATIONS  TO  OTHER  MATERIALS 

A.  The  (100)  Surfaces  of  Ge,  GaAs,  and  ZnSe 

In  this  section  we  apply  the  scattering-theoretic  method  to  the  (100)  surfaces  of  the 
isoelectronic  series  Ge,  GaAs,  and  ZnSe.  This  study  will  allow  us  to  address  a number  of 
issues  that  arise  in  the  application  of  the  method  to  compound  semiconductors.  The  materials 
were  chosen  in  order  to  carry  out  a study  of  the  variation  of  the  surface  electronic  structure 
with  increasing  ionicity.  Similar  studies  of  the  (110)  and  the  (111)  surfaces  of  this  isoelec- 
tronic series  have  been  done  by  others.^** 

We  again  use  ETBM  Hamiltonians  as  in  Sec.  III.  For  this  study  of  qualitative  trends,  we 
use  the  first-nearest  neighbour  Hamiltonians  given  by  Chadi^^  (Table  1). 

The  geometry  for  the  zincblende  semiconductors  is  the  same  as  for  Si  or  Ge  (Figs.  1-4) 
except  that  we  have  alternate  layers  of  anions  and  cations.  For  example,  to  create  a Ga-  or  an 
As-  terminated  surface  in  GaAs,  it  is  sufficient  to  remove  one  As  or  one  Ga  layer,  respectively. 
This  removal  yields  a 4x4  Green’s  function  matrix  and  creates  two  inequivalent  surfaces.^^ 
The  left-hand  semi-infinite  solid  is  rotated  about  the  x-axis  by  90°  compared  to  the  right-hand 
semi-infinite  solid  with  respect  to  their  common  y-z  coordinate-system.  As  a result,  the  zeros 
of  det  I G*  I yield  singly  degenerate  states  (except  at  T and  K)  corresponding  to  the  two 
surfaces.  The  states  on  the  two  opposite  surfaces  are  related  by  a simple  90°  rotation,  so  that, 
e.g.,  the  J point  of  the  one  surface  corresponds  to  the  J’  point  of  the  other,  and  vice  versa. 
This  complication  may  be  avoided  by  removing  two  layers,  namely  a Ga  and  an  As  layer.  The 
size  of  the  Green’s  function  matrix  is  increased  to  8x8,  but  the  energies  of  the  surface  states 
of  both  Ga-  and  As-  terminated  (100)  surfaces  are  obtained  from  the  zeros  of  one  determi- 
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nant.  The  states  of  the  two  surfaces  can  then  be  distinguised  by  either  examining  the  wave- 
functions  or,  even  more  simply,  by  performing  a complementary  calculation  with  4x4  matrices 
(only  one  layer  removed). 

The  surface  band  structure  (SBS)  for  Ge  (100),  calculated  in  the  first-nearest-neighbour 
approximation  has  already  been  given  in  Fig.  12  and  was  discussed  in  Section  111  D.  Our 
results  for  the  GaAs  and  the  ZnSe  (100)  surfaces  are  shown  in  Figs.  13  and  14,  where 
separate  plots  for  the  anion-  and  the  cation-terminated  surfaces  are  given. 

The  first  result  that  is  evident  from  Figs.  13  and  14  is  that  surface  states  of  an  anion- 
terminated  surface  tend  to  be  within  the  valence  band  region.  Conversely,  surface  states  of  a 
cation-terminated  surface  tend  to  be  within  the  conduction  bands.  This  result  correlates  nicely 
with  the  well  known  fact  that  bulk  valence  bands  are  dominated  by  anion-like  states  while 
conduction  bands  are  dominated  by  cation-like  states  (see  also  Fig.  IS).  It  also  demonstrates 
that  the  character  and  energetic  position  of  most  surface  states  depend  primarily  on  the  nature 
of  atoms  in  the  first  layer. 

If  we  recall  that  our  ETBM  Hamiltonians  describe  bulk  valence  bands  quite  accurately 
and  bulk  conduction  bands  rather  poorly,  it  follows  immediately  that  our  calculations  of  the 
SBS  of  the  ideal  anion-terminated  surfaces  are  quite  reliable,  whereas  the  corresponding 
calculations  of  the  SBS  of  the  ideal  cation-terminated  surfaces  are  qualitative  at  best.  We  will 
therefore  first  discuss  anion-terminated  surfaces  in  detail  and  then  make  a few  remarks  about 
cation-terminated  surfaces. 

Both  GaAs  and  ZnSe  anion-terminated  surfaces  have  a backbond-type  band  that  lies 
entirely  within  the  heteropolar  gap.  These  states  are  highly  localized  and  are  predominantly 
anion  s-like,  very  much  like  the  bulk  band  immediately  below  them.  These  characteristics 
become  more  prominent  as  one  goes  from  GaAs  to  ZnSe.  Similarly,  the  dispersion  of  the 


bulk  band  immediately  below  it. 


Both  GaAs  and  ZnSe  have  another  backbond-type  surface  state  that  appears  in  the  pocket 
I around  —4  to  —6  eV.  These  states  turn  out  to  be  more  spread  out  and  their  dispersion 

I decreases  much  more  slowly  along  the  isoelectronic  series. 

1 

We  turn  now  to  the  two  surface  bands  in  the  fundamental  gap.  As  in  the  case  of  Ge,  the 
I lower  band  has  a dangling-bond  character  formed  by  s and  p^  orbitals,  .whereas  the  upper  band 

has  a bridge-bond  Py-Pj  character.  In  fact,  the  separation  between  the  bridge  bond  and  the 
dangling  bond  may  be  largely  due  to  the  presence  of  the  s admixture  in  the  latter.  This  is 
consistent  with  the  fact  that  the  average  separation  between  the  two  types  of  states  is  smaller 
in  ZnSe  where  the  Se  atomic  s states  is  lower  in  energy  than  in  As  and  thus  does  not  admix  as 
much.  As  for  the  dispersion  of  the  two  types  of  states,  we  note  that  the  bridge-bond  band  is 
almost  dispersionless  and  lies  at  the  atomic  p level,  whereas  the  dangling-bond  band  has  more 
dispersion  through  interactions  mediated  by  the  second  layer. 

i 

Turning  to  the  cation-terminated  surfaces.  Figs.  13  and  14  reveal  that,  qualitatively,  a 
I similar  set  of  surface  bands  is  obtained.  As  noted  earlier,  the  bulk  conduction  bands  are  given 

I rather  poorly  by  the  ETBM  Hamiltonian  used  in  our  calculations.  The  two  surface  bands 

within  the  range  of  the  projected  conduction  bands  are  therefore  only  qualitatively  meaningful. 

I 

I 

The  dangling-bond  band  of  GaAs,  however,  which  lies  in  the  fundamental  gap  appears  to  agree 
rather  well  with  the  self-consistent  results  of  Appelbaum,  Baraff,  and  Hamann.^^  We  investi- 
gated the  origins  of  this  state  at  sample  3 points  by  calculating  AN^(E)  and  integrating  it  over 
the  projected  valence  bands.  The  resulting  integral  was  — I indicating  that  the  state  is  mainly 
derived  from  the  valence  bands,  very  much  like  the  dangling-bond  band  of  the  (100)  surfaces 
of  Si  and  Ge  and  the  anion-terminated  (100)  surfaces  of  GaAs  and  ZnSe.  This  finding 
suggests  that  the  dangling-bond  bands  even  of  cation-terminated  surfaces  can  be  calculated 
fairly  reliable  in  a first-nearest-neighbour  approximation.  Finally,  the  back-bond  states 
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appearing  in  the  valence  band  pocket  of  cation-terminated  surfaces  (Figs.  13  and  14)  are 
probably  derived  totally  from  valence  states  and  are  therefore  even  more  reliable. 

B.  The  (100)  Surface  of  SiOi 

In  this  section  we  study  the  Si-  and  the  O-  terminated  (100)  surfaces  of  the  ideal  cubic 
form  of  Si02.  i-e.  /3-cristobalite.  Si02  is  an  important  technological  material  and  is  widely  used 
as  a substrate  for  depositing  other  substances.  Knowledge  of  the  surface  electronic  structure 
of  this  material  might  help  understand  the  way  it  bonds  with  other  substances.  The  study  we 
present  here  is  somewhat  crude,  since  we  calculate  only  surface  states  with  the  ^-cristobalite 
geometry.  On  the  other  hand,  studies  based  on  /3-cristobalite-like  geometries  have  helped 
understand  most  electronic  properties  of  bulk  crystalline  and  amorphous  Si02.  We  therefore 
expect  that  the  present  study  of  Si02  surfaces,  which  is  the  first  to  be  attempted,  will  serve  as 
a guide  for  more  detailed  understanding  in  the  future. 

The  first-nearest-neighbour  parameters  (Table  1)  used  in  the  bulk  ETBM  Hamiltonian, 
were  found  by  fitting  to  the  Si02  bulk  valence  bands  given  in  Ref.  (68)  and  to  a band  gap  of  9 
eV.^^  The  valence  bands  and  the  band  gap  are  reproduced  accurately  by  retaining  only  the 
Si-O  and  the  0-0  first-nearest-neighbour  matrix  elements.  The  Si-Si  interactions  are  left  out. 
The  geometrical  arrangement  of  the  atoms  for  a Si-terminated  (100)  surface  of  Si02  is  shown 
in  Fig.  16,  where  we  have  depicted  the  first  8 layers,  including  the  surface  layer.  The  two- 
dimensional  unit  cell  is  again  marked  by  a dashed  line.  It  should  be  noted,  that  the  unit  cell 
contains  one  atom  in  each  Si  layer  but  two  atoms  in  each  O-layer.  The  point-group  symmetry 
(C2v)  and  the  SBZ  are  the  same  as  in  the  zincblende  materials  (Fig.  3). 

The  creation  of  two  adjacent  Si-terminated  surfaces  is  accomplished  by  removing  one  O 
layer,  yielding  an  8x8  Green’s  function  matrix  (2  atoms  per  unit  cell).  The  creation  of  two 
adjacent  O-terminated  surfaces  cannot,  however,  be  accomplished  by  removing  only  one  Si 
layer,  because  the  resulting  two  semi-infinite  solids  are  still  coupled  via  the  first-nearest- 


neighbour  0-0  interactions.  In  order  to  accomplish  the  desired  decoupling,  at  least  two  layers 
(one  Si  and  one  O layer)  must  be  removed.  The  result  is  a 12x12  Green’s  function  matrix 
describing  one  Si-terminated  surface  and  one  O-terminated  surface.  The  calculated  states  can 
then  be  sorted  by  comparing  with  the  previously  obtained  states  for  the  Si-terminated  surface. 

The  projected  valence  bands  of  Si02  are  shown  in  Fig.  17.  The  upper  part  of  the 
projection,  with  a width  of  about  3 eV,  results  mainly  from  the  nonbonding  or  "lone  pair"  O 
2p-orbitals  which  lie  perpendicular  to  the  Si-O-Si  chain.  The  lower  part,  with  a width  of  about 
7 eV,  results  mainly  from  the  bonding  O 2p-orbitals,  which  lie  along  the  Si-O-Si  chain.  The 
bonding  bands  have  the  same  overall  structure  as  the  valence  bands  of  Si.  The  similarity  in  the 
corresponding  parts  of  the  projected  band  structures  may  be  seen  by  compairing  Figs.  17  and 
5. 


The  surface  states  for  the  Si-  and  O-terminated  (100)  surfaces  of  cubic  SiOj  are  shown  in 
Figs.  17a  and  b.  No  bound  surface  states  were  found  in  the  optical  gap  for  either  surface. 
The  nonbonding  bands  are  affected  very  little  by  the  creation  of  the  surface.  No  dangling- 
bond  bands  are  split  off,  as  expected,  because  there  are  no  bonds  to  begin  with.  Instead,  the 
states  in  the  non-bonding  bands  are  very  much  atomic  in  character,  so  that  the  creation  of  the 
surface  is  only  a weak  perturbation,  as  far  as  the  lone  pair  orbitals  are  concerned.  The 
bonding  bands,  on  the  other  hand,  show  pronounced  surface  states,  as  one  might  expect,  since 
bonds  are  broken  by  the  creation  of  the  surface. 

For  the  Si-terminated  surface  (Fig.  1 7a)  we  find  a number  of  surface  states  which  are 
predominantly  localized  on  the  O 2p-orbitals  in  the  O layer  immediately  below  the  surface  Si 
layer.  They  are  therefore  best  described  as  backbond  states.  The  Si  dangling-bond  states  are 
high  in  the  conduction  bands.  Most  of  the  backbond  states  lie  very  close  to  the  bulk  band 
projections  and  show  similar  dispersion.  This  behaviour  may  be  explained  by  the  fact  that  the 
O atoms  in  the  strengthened  backbonds  retain  all  their  Si  neighbours  and  four  out  of  their  six 
O neighbours. 
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In  the  case  of  the  O-terminated  (100)  surface  (Fig.  17b),  the  perturbation  has  an  even 
stronger  effect.  Two  of  the  six  O neighbours  and  one  of  the  two  Si  neighbours  of  each 
surface-layer  O atom  are  removed.  The  2p-orbitals  on  the  surface-layer  oxygen  atom  are, 
therefore,  left  dangling,  giving  rise  to  two  dangling-bond  like  bands  at  around  —4  eV  in  the 
gap  of  the  PBS  between  the  projection  of  the  non-bonding  and  the  bonding  bands.  These  two 
dangling-bond  bands  (resulting  from  the  two  atoms  in  an  O layer)  are  almost  degenerate  at  T, 
but  split  slightly  at  other  g-points,  because  the  two  2p  orbitals  point  in  different  directions. 
Finally,  the  lower  lying  bands  have  again  back-bond  character. 

The  most  significant  of  the  above  results  is  the  finding  that  oxygen  atoms  at  the  free 
surface  give  rise  to  occupied  dangling-bond  states  which  lie  in  the  gap  between  the  bonding 
and  non-bonding  valence  bands.  This  finding  is  independent  of  our  Hamiltonian  parametriza- 
tion  and  special  crystal  structure,  since  the  dangling  bond  states  must  lie  in  the  gap^^  above 
the  bonding  bands  and  below  the  non-bonding  bands  (the  energy  of  a Si-O  bond  is  lower  than 
the  energy  of  an  O 2p-orbital,  and  the  energy  of  a Si-O-Si  bond  is  even  lower).  These  states 
should  be  observable  and  can  probably  be  detected  by  high  resolution  experiments. 

Conclusions 

We  have  demonstrated  that  the  Koster-Slater  method  for  dealing  with  localized  perturba- 
tions can  be  formulated  in  a way  that  allows  the  study  of  the  electronic  structure  of  surfaces  of 
real  materials.  We  have  pointed  out  the  various  advantages  of  the  method,  when  used  in 
conjunction  with  tight-binding  Hamiltonians.  We  illustrated  applications  of  the  method  to  the 
free  (100)  surfaces  of  Si,  Ge,  GaAs,  ZnSe  and  cubic  Si02.  The  method  would  be  particularly 
suited  for  the  study  of  reconstructed  surfaces,  where  large  unit  cells  are  involved.  Such  studies 
have  not  been  persued  yet.  Instead,  we  have  modified  the  method  to  treat  interfaces  between 
two  semi-infinite  crystalline  solids.  This  work  will  be  reported  in  a separate  paper.  Finally  the 
method  has  the  potential  of  being  a very  efficient  way  of  carrying  out  self-consistent  calcula- 


tions of  the  electronic  structure  of  surfaces  and  interfaces. 
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Appendix  A 


Evaluation  of  the  Green’s  Function  Matrix  Elements 


The  matrix  elements  of  the  Green's  function  G®(E)  (equation  (3))  in  the  layer-orbital 
representation  (equation  (23))  are  given  by 


G)/  (E.  q) 


2 

nk 


</g|  nkXnk  I /'£> 


(A.l) 


E+-E„(k) 

A 

where  | nk>  are  the  Bloch  functions  [eq.  (22)]  and  1 /g>  are  the  layer-orbitals  defined  by  eq. 
(23). 


In  order  to  evaluate  the  matrix  elements  in  (A.l),  we  must  introduce  further  notation. 
We  first  rewrite  the  bulk  atomic  position  vectors  in  terms  of  two-dimensional  Bravais  lattice 
vectors  £j  plus  appropriate  basis  vectors  in  the  form 


Rj  I,  - £j  + 2' 

Furthermore  we  need  to  decompose  the  basis  vector  X™  into  a surface-parallel  component 
and  a perpendicular  component  k"': 

X^  - a™  + k"!  (A.3) 

Note  that  for  any  bulk-atom,  £j  locates  a two-dimensional  Bravais  lattice  point,  k"*  locates  a 
particular  plane,  and  finally  a™  defines  the  atomic  position  in  the  two-dimensional  unit  cell  in 
that  plane. 

AkMi|  the  same  lines,  the  bulk  k vector  must  be  decomposed  into  surface-parallel  and 
perpeadiculsr  components: 


k - (g-Hg)  + kj. 


(A.4) 
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Note  that  the  surface-parallel  component  has  been  written  as  the  sum  of  a £ vector  in  the 
surface  Brillouin  zone  (SBZ)  and  a two-dimensional  reciprocal  lattice  vector  g,  because  the 
projection  of  k onto  the  surface  does  not  necessarily  lie  within  the  two-dimensional  Brillouin 
zone. 

With  the  above  decompositions,  the  matrix  elements  appearing  in  (A.l)  can  be  evaluated. 
The  result  is 

G?/'(E.£)-^  2 2 2'  — — I PV(ki.a+g)-  (A.5) 

n kj,  g E+-En(kj_,£-)-s) 

The  auxiliary  function  P is  given  by 

P?/"  (•'i.  a+g)  = a+g)  X 

X exp  )1  exp  [ik^^fx"’-*'"  )1  . (A.6) 

% 

Note  that  the  Bloch-function  coefficients  C now  have  the  double  index  m,  n in  place  of  the 
single  index  v appearing  in  eq.  (22)  because  of  the  decomposition  given  in  (A2).  Finally,  the 
prime  on  the  g-sum  in  eq.  (AS)  means  that  the  sum  is  carried  out  only  over  those  g vectors 
neccessary  to  cover  the  projection  of  the  bulk  Brillouin  zone  onto  the  surface  as  required  by 
the  decomposition  (A4). 

For  the  (100)  surface  of  a diamond-  or  zince-blende-type  material,  the  projection  of  the 
bulk  Brillouin  zone  is  shown  in  Fig.  A1  together  with  the  SBZ.  Also  shown  is  the  irreducible 
segment  of  the  SBZ,  marked  by  the  four  high-symmetry  points  F,  J,  J',  and  K.  For  £ vectors 
inside  the  cross-hatched  region  of  the  SBZ,  no  g vectors  contribute  to  the  sum  in  eq.  (AS). 
For  £ vectors  in  the  remainder  of  the  irreducible  square  there  exist  non-trivial  contributions  to 
the  g-sum  in  eq.  (AS).  This  means  that  one  must  sum  over  kj^^k^  both  at  the  £ point  itself 
and  at  the  corresponding  £-t-g.  as  shown  by  the  hatched  regions  in  Fig.  AI. 
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For  g points  inside  the  cross-hatched  region  of  the  SBZ,  the  sum  over  extends 

from  -2ir/a  to  -t-2ir/a,  where  a is  the  lattice  constant,  as  illustrated  for  g points  along  the  FJ 
liki'  in  Fig.  A2.  The  figure  is  simply  the  cut  of  the  bulk  Brillouin  zone  by  the  plane  which  is 
perpendicular  to  the  surface  and  contains  the  line  FJ.  For  g-points  outside  the  cross-hatched 
region  of  the  irreducible  square,  the  limits  at  g and  g-fg  are  different  and  g-dependent,  as 
illustrated  in  Fig.  A2  for  one  such  point  q.  By  making  use  of  symmetry  considerations,  it  can 
be  shown  that  instead  of  summing  over  the  appropriate  ranges  of  k,  at  both  g^  and  go't'S*  the 
same  result  may  be  obtained  by  simply  summing  only  at  g^  over  the  range  — 2ir/a  to  2v/a 
(Fig.  A2).  This  alternative  way  of  evaluating  (he  sums  in  eq.  (AS)  is  a substantial  simplifica- 
tion especially  for  g points  away  from  high-symmetry  lines  in  the  SBZ. 

For  the  (100)  surfaces  of  diamond-  and  zinc-blende-type  materials,  k"'  > am/4  and  eqs. 
(AS)  - (A6)  simplify  to 

, -1  C?Mxq)C?(xq) 

G)(>  (q,E)  ■ 2 I dx  exp  [iffx(m-m')/2]  , (A.7) 

“ ^ n ri  E+-En(x.g) 

where  x^k,  a/2tr.  Before  this  expression  can  be  evaluated,  one  must  deal  with  the  fact  that 
E*^  stands  for  E-fie,  and  the  limit  e->0  is  to  be  taken.  In  the  gaps  and  pockets  of  the  PBS 
one  can  simply  set  < -0  and  evaluate  G}f>  (q,  E)  without  any  complication,  because  the 
denominator  in  (A7)  is  never  zero.  In  addition,  Gff  is  Hermitian,  which  reduces  the  compu- 
tational effort  even  further.  Within  the  PBS  continua,  however,  e cannot  be  set  to  zero. 
There  are  two  ways  to  proceed;  The  most  direct  way  is  to  set  e to  a small  but  finite  energy 
and  evaluate  eq.  (A7)  as  it  stands.  The  result  will  be  a Green’s  function  which  is  broadened 
by  the  amount  e . Alternatively,  one  introduces  the  spectral  function 


SfC  (g.  E)  - 2 </g  I nkxnk  | f'g>i(E-EnJ 
nk  ~ “ 


(A.8) 


and  takes  the  limit  c -*0  analytically  resulting  in 

S//'(g,E) 


G?^M3.E)-  f 


E-E' 


dE'  - i»  S//'(g,E) 


(A.9) 
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where  the  bar  on  the  integral  sign  denotes  principal  value.  can  be  evaluated  directly.  The 
result  is 


S(C  (q.  E)  - ^^2  2 2 2' 
Nj  n i 2 

where  kj  are  defined  by 


P?/Mki.a+g)  ^ 


aEn(kj^.q+g) 

bk , 


c^-kj  ] 


-1 


(A.  10) 


E - E„  (ki.  9+a)  . 


(A.11) 


For  our  particular  case,  the  final  result  for  SfC  is  (Xj  - k^  2e/a) 


S//'(q,  E) 


exp  fiirxj(m-m')/21  2 
n 


C?*  (xj.a)  C?  (Xj.g)  [ 


&En(x,9+g) 

b. 


(A.12) 


where  Xj  is  defined  by 


and  -l<xj<l. 


E - E„  (Xj.  g) 


(A.13) 
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Appendix  B 

Search  for  the  Bound  State  Energies 

The  energies  of  bound  surface  states  are  determined  by  the  zeros  of 

D(E)  » det  a G}f  (E.  g)  I - 0 . (B.l) 

where  ( and  l'  run  over  the  layers  which  have  to  be  removed  for  the  creation  of  the  surface. 
(See  Appendix  A for  the  evaluation  of  the  matrix  Gf  j’).  The  determinant  in  (B.  may  be 
evaluated  numerically  in  a variety  of  ways  but  the  search  for  the  zeros  could  be  a time- 
consuming  task.  Instead,  the  zeros  can  be  located  much  more  efficiently  by  first  converting 
the  matrix  G°i(',  to  triangular  form  and  counting  the  number  of  negative  diagonal  elements.  If 
this  number  is  different  at  the  two  ends  of  a given  energy  range,  the  difference  corresponds  to 
the  number  of  zeros  of  D(E)  in  that  range.  The  chosen  range  is  then  systematically  bisected 
until  the  positions  of  zeros  are  determined  to  the  desired  degree  of  accuracy.  In  our  calcula- 
tions the  surface  states  of  the  chosen  Hamiltonians  were  determined  with  an  accuracy  of  I 
meV. 
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Appendix  C 


Calculation  of  the  Wavefunctions 

The  eigenvectors  of  the  Hamiltonian  H are  given  by  the  Lippmann-Schwinger  equation 


m ip  + G + (C.l) 

where  9>  is  the  solution  of  H",  i.e.  <i>  stands  in  our  problem  for  a bulk  Bloch  function.  In  the 
pockets  and  gaps  of  the  PBS,  where  the  true  surface  states  are  encquntered,  there  exists  no 
nontrivial  solution  of  h”.  Therefore,  the  wavefunction  for  the  surface  states  is  given  as 

(C.2) 

In  the  surface  problems  discussed  in  this  paper,  the  perturbation  operator  acts  only  on  a very 
limited  part  of  the  Hilbert  space  of  H*,  as  was  extensively  discussed  in  Section  II.  We, 
therefore,  partition  the  perturbation  operator,  the  Green’s  operator  and  the  wavefunction  in 
the  two  subspaces  of  H*,  which  we  label  A and  B: 


U - 


' ^BA  ^BB  ' ' V'B  ' 


(C.3) 


where  is  the  unit  operator  in  the  subspace  A.  With  (C.3)  inserted  in  (C.2)  we  find 


'^'A  * “ "^A  (C.4a) 

'/'B  “ u V'A  • (C.4b) 

Eq.  (C4a)  may  also  be  written  as 

(-i-lA-G;y,)(u^A)-0-  ^C5a) 

This  is  a homogeneous  set  of  linear  equations  and  can  be  solved  for  (u^^)  arbitrary 

value  of  u.  In  the  limit  u-*ae  one  simply  diagonalizes  and  the  resulting  eigenvectors 


.1 
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are  finite  so  that  ^^>0,  as  expected.  The  quantities  are  then  inserted  in  (C4b)  and  i/'u 
can  be  directly  calculated  by  simply  multiplying  the  matrix  Gba  by  the  column  vector  (uV'^)- 
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Appendix  D 


The  Solution  of  The  Dyson  Equation 
The  formal  solution  of  the  Dyson  equation  (13) 

G-G°  + G°UG  (D.l) 


can  conveniently  be  written  in  terms  of  the  scattering  matrix 

T-U(1-G*U)-'  ‘ (D.2) 


so  that 


G-G”  + G°TG”.  (D.3) 

Equation  (D.3)  shows,  that  the  complete  information  about  the  new  features  of  the  perturbed 
system  H a U is  contained  in  the  scattering  matrix  T.  Due  to  the  limited  range  of  U,  the 
scattering  matrix  can  be  calculated  exactly.  To  calculate  T,  we  have  to  evaluate  (1-G°U)~*. 
We  introduce  the  notation 


Q - I -G  U 


(D.4) 


and  employ  the  same  partitioning  convention  as  in  (C.3).  The  matrix  Q is  of  the  form 


/Qaa 

'^Qba 


' ) 


(D.5) 


and  can  be  inverted  analytically,  even  though  Qg^  and  Ig  are  infinite  submatrices.  The  result 
is 


Q-i 


/IQaaI-' 

' “*B  Qba  WaaI"' 


(D.6) 


Using  this  result  we  obtain: 
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Gaa  - g.Ik  + g;^  Uaa  iqaaI-'  g;u 

Gab  “ G\b  + [Qaa1~*  G.^b  (D.7) 

Gba  “ Gba  + Gba  Gaa  IQaaI”’  Gaa 
Gbb  “ Gbb  + Gba  ^aa  IGaaI”'  Gab 

This  result  holds  for  any  perturbation  confined  to  subspace  A.  Note  that  is  precisely  the 
matrix  one  has  to  evaluate  in  order  to  search  for  bound  states.  The  complete  evaluation  of  the 
Green’s  function  matrix  for  the  perturbed  system  can  then  be  done  directly  by  inverting  the 
small  matrix  and  multiplying  matrices  as  required  by  (D.7). 

The  above  results  simplify  further  in  the  case  of  the  perturbation  that  creates  surfaces 
used  in  this  paper.  The  perturbation  is  of  the  form 

Uaa  - u 1a  (D.8) 

where  u is  a constant  that  is  eventually  let  go  to  infinity.  Eq.  (D7)  then  becomes 

Gaa-0 

Gab  - 0 (D.9) 

Gba  “ 0 

Gbb  “ Gbb  ~ Gba  IGaaI"'  Qab  • 

These  equations  are  valid  for  all  finite  energies^**.  Note  that  Gaa-0  corresponds  to  the  fact 
that  there  are  no  states  at  finite  energies  in  the  space  occupied  by  the  "removed"  layers. 
Gab-0  and  Gba-0  correspond  to  the  fact  that  the  "removed"  layers  are  completely  decou- 
pled from  the  two  semi-infinite  solids  described  by  Gbb-  particular  interest  that  Gbb 
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given  entirely  in  terms  of  G°,  which  describes  the  infinite  bulk  solid.  This  result  is  of  course 
true  because  no  electronic  or  lattice  relaxation  is  included  at  the  surface. 

Gqb  contains  a wealth  of  information.  In  particular,  one  can  evaluate  the  local  density  of 
states  at  the  layer  (with  / in  subspace  B)  as  follows: 

N/ (a.  E)  - • Im  Tr  jG?/-(a.  E)-2  G^j  [G^aIT/g]/'!  (D.IO) 

ir  t ij  I 

where  i and  j run  over  the  layer  orbitals  in  the  "removed"  layers  in  the  subspace  A.  The 
second  term  in  the  curly  brackets  clearly  yields  directly  the  change  in  the  local  density  of  states 
for  the  /'*'  layer  (The  sum  of  these  quantities  over  all  layers  is  equal  to  AN(q,  E)  discussed  in 
Sec.  III). 
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Table  1 


Tight  binding  matrix  elements  in  eV  defining  the  bulk  Hamiltonians  of  the  various 
materials  in  standard  notation.  The  superscripts  (c)  and  (a)  in  the  headings  of  the  first  four 
columns  denote  cation  and  anion,  respectively,  and  the  superscripts  1 and  2 refer  to  first-  and 
second-nearest-neighbours. 


E»»> 

V' 

ss 

V' 

V' 

''pa 

v' 

PPO 

V' 

PP* 

PP» 

PP* 

Si“ 

-4.203 

0.187 

/ 

/ 

-2.08 

2.12 

2.12 

2.32 

-0.52 

0.58 

-O.IO 

Ge* 

-5.830 

0.610 

/ 

/ 

-1.69 

2.03 

2.03 

2.55 

-0.67 

0.41 

-0.08 

Ge^* 

-5.63 

1.72 

/ 

/ 

-1.69 

2.76 

2.76 

3.13' 

-0.92 

/ 

/ 

GaAs** 

-3.07 

3.47 

-8.09 

1.2« 

-1.69 

2.37 

2.06 

3.51 

-0.96 

/ 

/ 

ZnSe*» 

+0.01 

6.20 

-12.03 

1.10 

-1.69 

2.59 

1.07 

3.46 

-0.75 

/ 

/ 

SiO, 

13.86 

18.36 

16.36 

-1.77 

-1.5 

3.76 

3.5 

5.71 

-0.64 

/ 

/ 

0-0  Interactions 

/ 

/ 

/ 

-0.6 

0.8 

0.8 

1.29 

-0.16 

/ 

/ 

a.  Ref.  58 


b Ref.  63 
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^AA  poles  at  infinity  corresponding  to  the  states  on  the  removed  layers. 


! 

\ 

[ 

i 

\ 

I 


2 


o 

o 

© 

0 

o 

o 

o 

o 

© 

© o 

o 

o 

o 

V 

1 

o 

o © 
o © 
o © 


u-»fl0 


o o 
o o 
© o 


Schematic  graphs  showing  the  creation  of  the  twin  surfaces  in  the  two  different  methods;  (a) 
"switching  off"  interactions,  (b)  "removal"  of  layers.  In  both  cases  the  resulting  surface  layer 


atoms  are  shaded. 
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Fig.  2 


Fig.  3 
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The  topmost  (100)  surface  layer  and  the  following  three  layers  for  diamond  or  zincblende 
lattices.  In  the  case  of  a zincblende  material  even  and  odd  layers  contain  either  anions  or 
cations,  respectively.  We  retain  the  crystal  y-  and  z-axis.  The  bulk  lattice  constant,  denoted 
by  a. 


Reciprocal  lattice  (dots)  for  a diamond  or  zincblende  (100)  surface  together  with  the  surface 
Brillouin  zone.  The  irreducible  square  (Ciy  symmetry)  is  cross-hatched. 
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Fig.  4 Geometry  in  a plane  prependicuiar  to  the  surface  containing  the  y-axis.  The  two  dotted  layers 

are  "removed"  to  create  the  adjacent  surfaces. 


PBS  a SURFACE  STATES  OF  Si  (100) 


Projected  bulk  band  structure  and  surface  band  structure  of  Si  (100)  calculated  with  the 
sccond>nearest*neighbour  tight-binding  Hamiltonian  given  in  Ref.  S8.  The  letters  b.-d,  and  br 
denote  the  back  bond,  dangling  bond  and  bridge  bond  bands,  respectively.  See  text. 
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Wavefunction  amplitudes  as  a function  of  the  layer  number  for  some  surface  states  in  Si  (100) 
(summed  over  the  s and  p orbitals).  The  corresponding  3 and  E values  are  explicitly  given. 
The  energy  values  resulting  from  a 20-layer  slab  calculation  (Ref.  38)  are  also  given  in 
parenthesis  for  comparison. 
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(«V) 


Fig.  1 1 (a)  Local  densities  of  states  (in  units  of  states/eV/unit  cell)  at  each  one  of  the  first  seven 

layers  of  the  Si  (100)  surface  at  the  symmetry  point  J of  the  surface  Brillouin  zone.  The  local 
density  of  states  at  any  bulk  layer  of  the  infinite  crystal  is  shown  by  the  heavy  line,  (b)  The 
changes  in  the  local  density  of  states  at  each  of  the  first  seven  layers.  The  total  change  is 
given  by  the  heavy  line  and  is  the  same  as  that  shown  in  Fig.  9,  except  for  different  broaden- 
ing. All  the  curves  in  this  figure  are  broadened  by  0.3  eV.  The  vertical  scale  applies  in  (a) 
and  (b)  for  the  first  layer  only.  The  other  curves  are  all  shifted  down  by  .4,  .8  etc. 


Fig.  12  Projected  bulk  band  structure  and  surface  band  structure  of  Ge  (100)  calculated  with  the 

first-nearest-neighbour  tight-binding  Hamiltonian  as  given  in  Ref.  63.  The  letters  b.  d and  hr 
again  denote  the  back  bond,  dangling  bond  and  bridgebond  bands,  respectively.  The 
"atomic"  energy  levels  are  shown  for  convenience.  Note  that  the  bridge-bond  band  is  at  the 
3p  atomic  level  in  a first-nearest-neighbour  approximation. 
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Fig.  13 
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Surface  band  structure  for  the  anion-  and  cation-terminated  (100)  surfaces  of  GaAs.  For 
details  see  Fig.  caption  12.  The  dots  and  stars  in  Fig.  13b  are  the  selfconsistent  results,  given 
in  Ref.  28,  for  the  dangling  bond  and  bridge  bond,  respectively. 
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Fig.  14  Surface  band  structure  for  the  anion-  and  cation-terminated  (100)  surfaces  of  ZnSe.  For 


details  see  Fig.  caption  12. 


Projected  valence  bands  and  surface  band  structure  for  the  Si-  and  O-  terminated  (100) 
surfaces  of  cubic  Si02.  The  energetic  position  of  the  "atomic"  02p  orbitals  is  indicated. 


The  dots  at  F and  at  1,  2,  3 and  4 are  reciprocal  lattice  points  for  the  Si  (100)  surface.  The 
heavy  line  shows  the  Surface  Brillouin  zone.  The  square  denoted  by  F,  J,  K,  J*  is  the 
irreducible  part  of  the  surface  Brillouin  zone.  When  g runs  over  this  irreducible  part,  the 
vectors  g -f  gj  for  i«l.  2,  3,  4 (g-  being  surface  reciprocal-lattice  vectors)  run  over  the  four 
squares  marked  by  the  corresponding  g-f  gj-  The  dashed  octagon  is  the  projection  of  the  bulk 
Brillouin  zone  onto  the  (100)  surface  and  the  dashed-dotted  square  is  the  projection  of  that 
part  of  the  bulk  Brillouin  zone,  for  which  bolds.  The  details  concerning  the 

cross-hatched  and  the  shaded  areas  are  discussed  in  Appendix  A. 


' 

Appendix  K 

RC  7200  (#30973)  6/26/78 
Solid  State  Physics  7 pages 

Research  Report 

NEW  METHOD  FOR  THE  ELECTRONIC  STRUCTURE  OF  HETEROJUNCTIONS  - 
APPLICATION  TO  THE  (100)  Ge-GaAs  INTERFACES* 

Sokrates  T.  Panteiides  and  J.  Pollmann^ 

IBM  Thomas  J.  Watson  Research  Center 
Yorktown  Heights,  New  York  10598 

Typed  by  Linda  P.  Rubin  (SP.2162) 


1 


^twifch  Oiviiioo 

S«n  JoM  Yorktown  Zurich 


RC  7200  (#30973)  6/26/78 
Solid  State  Physics  7 pages 

NEW  METHOD  FOR  THE  ELECTRONIC  STRUCTURE  OF  HETEROJUNCTIONS  - 
APPLICATION  TO  THE  (100)  Ge-GaAs  INTERFACES* 

Sokraies  T.  Panielides  and  J.  Pollmann^ 

IBM  Thomas  J.  Watson  Research  Center 
Yorktown  Heights,  New  York  I0S98 

Typed  by  Linda  P.  Rubin  (SP.2 162) 

ABSTRACT:  A new  method  which  for  the  first  time  can  treat  two  truly  semi-infinite  semicon- 
ductors in  contact  is  introduced  and  used  to  study  the  electronic  structure  of  the  (100) 
Ge-GaAs  interfaces.  Both  the  Ge-Ga  and  the  Ge-As  interfaces  exhibit  essentially  three 
interface  bands.  The  nature  and  origin  of  these  bands  are  discussed  in  detail  in  terms  of  local 
densities  of  states.  The  results  are  used  to  obtain  a new  interpretation  of  the  experimental 
data. 


* Work  supported  in  part  by  the  office  of  Naval  Research  under  contract  No.  N00014-76- 
C-0934. 

t Permanent  address:  Institut  fUr  Physik  der  Universitat  Dortmund,  Dortmund,  W. 
Germany. 


LIMITED  DISTRIBUTION  NOTICE 


This  report  hss  been  subnlcted  for  publication  elsewhere  and 
has  been  Issued  as  a Research  Report  for  early  dissemination 
of  Its  contents.  As  a courtesy  to  the  intended  publisher,  it 
should  not  be  widely  distributed  until  after  the  date  of  outside 
publication. 


Copies  may  be  requested  from: 

IBM  Thomas  J.  Watson  Research  Center 

Post  Office  Box  218 

Torkcown  Heights,  New  York  10S98 


Page  I 


r 

}• 

I 


Interfaces  between  two  different  semiconductors,  often  referred  to  as  heterojunctions, 
have  been  investigated  extensively  in  connection  with  device  fabrication.'  Theoretical  studies 
of  the  electronic  structure  of  realistic  models  of  heterojunctions,  however,  have  been  carried 
out  only  during  the  last  year.  Methods  used  thus  far  simulate  the  solids  in  contact  with  either  a 
finite-thickness  slab  on  a semi-infinite  substrate,^  or  with  two  finite-thickness  slabs  in  contact^, 
or  with  a superlattice  of  such  finite-thickness  slabs.'* 

i 

In  this  paper,  we  report  the  development  and  first  applications  of  a new  method  which  for  ^ 

I 

the  first  time  describes  the  electronic  structure  of  the  interface  between  two  truly  semi-infinite  ] 

I 

semiconductors.  We  use  the  method  to  give  a detailed  description  of  the  electronic  structure  ! 

of  the  (100)  Ge-GaAs  interfaces.  We  compare  our  results  for  the  Ge-Ga  interface  with 
previous  theoretical  work  (no  other  calculations  have  been  reported  for  the  Ge-As  interface). 

Finally  we  use  our  results  to  obtain  a new  interpretation  of  available  experimental  data. 

The  present  method  is  based  on  the  old  Koster-Slater^  idea  for  the  description  of  localized 
perturbations  of  otherwise  infinite  periodic  solids.  That  idea,  originally  proposed  for  point  j 

defects  has  since  been  extended  to  describe  free  surfaces,^^  which  may  be  viewed  as  arising  | 

from  a perturbation  that  has  two-dimensional  periodicity  and  is  localized  in  the  third  direction. 

( 

We  will  demonstrate  that  the  creation  of  an  interface  between  two  crystalline  semiconductors  i 

can  rtlso  be  described  by  a localized  perturbation  if  the  unperturbed  system  is  taken  to  be  two  j 

] 

infinite  non-interacting  solids.  j 

i 

I 

In  order  to  describe  the  local  perturbation  that  creates  the  interface,  we  use  the  Ge-GaAs  j 

■j 

(100)  interface  as  an  example  and  describe  the  two  unperturbed  bulk  materials  by  nearest- 
neighbor  tight-binding  Hamiltonians.^"  Each  (100)  layer  of  atoms  in  GaAs  has  either  Ga  or 
Ai  atoms,  so  that  we  can  schematically  denote  an  entire  layer  of  atoms  with  the  proper 
geometry  by  a single  circle.  Fig.  1 shows  the  local  perturbation  that  creates  the  interface: 

One  As  and  one  Ge  layer  is  removed '2  (u  discussed  extensively  in  Ref.  8 for  the  case  of  free 
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surfaces),  while  at  the  same  time  one  of  the  two  resulting  "free-surface"  Ge  layers  is 
"connected"  with  the  appropriate  "free-surface"  Ga  layer,  as  shown.  The  net  result  of  this 
perturbation  is  a system  which  consists  of  three  decoupled  subsystems*^:  a Ge-GaAs  interface 
between  two  semi-infinite  solids,  a semi-infinite  Ga-terminated  GaAs  crystal  and  a semi- 
infinite Ge  crystal. 

Once  the  interface  is  described  by  a localized  planar  perturbation,  the  full  machinery  of 
the  Koster-Slater  Green's  function  formulation  can  be  used^.  The  main  advantages  of  this 
approach  are:  (a)  it  deals  with  semi-infinite  solids;  (b)  it  exploits  the  fact  that  there  exists  an 
unperturbed  system  which  has  full  three-dimensional  translation  periodicity  whose  solutions  are 
easily  obtainable  and  whose  Green’s  functions  can  be  calculated  accurately  using  standard 
band-theoretic  techniques.  Thus,  all  properties  of  the  bulk  solids  are  built  in  from  the  start 
and  preserved,  and  changes  in  the  electronic  structure  are  calculated  directly,  rather  than  by 
subtracting  two  large  quantities,  as  in  the  case  of  slab  approaches.  The  method  yields  the 
positions,  the  dispersion  and  the  wavef unctions  of  the  bound  states  in  the  gaps  and  pockets  of 
the  joint  projected  band  structure  (PBS),  the  changes  in  the  density  of  states  within  the 
projected  band  continua  and  local  densities  of  states  (LDOS)  at  each  layer  of  atoms. 

As  a first  application  of  the  new  method  for  interfaces  we  have  studied  the  electronic 
structure  of  the  (100)  Ge-GaAs  and  Ge-ZnSe  interfaces.  We  have  used  the  tight-binding 
Hamiltonians  of  Chadi^  which  use  only  s and  p orbitals  on  every  atom  and  retain  only 
nearest-neighbor  Hamiltonian  interactions**^.  The  interface  matrix  elements  (Ge-Ga  hopping 
integrals)  were  chosen  to  be  the  averages  of  the  corresponding  Ge  and  GaAs  bulk  matrix 
elements.  In  this  communication,  we  discuss  only*^  the  results  for  Ge-GaAs  for  which  other 
calculations  are  available  for  immediate  comparison.^-'* 

In  Fig.  2 we  show  the  joint*^  PBS  of  Ge  and  GaAs  together  with  the  bound  states  for  the 
Ge-Oa  and  Oc-Aa  interfaces.  It  is  seen  that  very  few  bound  states  appear  in  the  gaps  and 
pockeu  of  the  joint  PBS  and  that  the  bound  states  are  in  fact  very  near  band  edges.  Such  a 
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situation  is  consistent  with  the  fact  that  the  two  solids  are  connected  without  lattice  mismatch 
and  have  very  similar  electronic  structures,  so  that  the  junction  is  a weak  perturbation.  Figure 
2 by  itself,  however,  does  not  provide  a complete  description  of  all  the  changes  in  the  electron- 
ic structure  produced  by  the  formation  of  the  heterojunction.  Changes  also  occur  within  the 
band  continua,  the  dominant  features  of  which  are  usually  referred  to  as  resonances  or 
antiresonances.  In  Fig.  3 we  show  the  y-integrated  (i.e.,  summed  over  the  two-dimensional 
Brillouin  zone)  LDOS  on  the  two  layers  next  to  each  interface  and  compare  them  with 

A 

9-integrated  LDOS  on  the  corresponding  bulk  layers.  (The  ^-integrated  LOOS  on  other  layers 
is  virtually  identical  with  the  corresponding  bulk  LDOS  indicating  strong  localization  of  the 
interface  effects).  Fig.  3 reveals  that  both  heterojunctions  give  rise  to  essentially  three 
interface  bands,  (see  also  Fig.  2) 

In  the  case  of  the  Ge-Ga  interface  (Figs.  2a  and  3a),  the  three  bands  lie  immediately 
aboye  bulk  bands,  and  their  dispersion  follows  closely  the  nearest  band  edge.  In  fact.  Fig.  3 
shows  that  each  interface  band  is  actually  derived  from  bulk  states  immediately  below  it. 
Physically,  what  happens  is  that  bulk  propagating  states  decay  at  the  interface  and  have  little 
or  no  amplitude  on  the  interface  layer.  (The  situation  corresponds  to  localized  antiresonanc- 
es). The  missing  amplitude  is  shifted  up  in  energy  into  new  states  that  are  localized  at  the 
interface  layers  and  are  either  bound  states  (if  in  gaps  or  pockets)  or  resonances  (if  within  the 
band  continua).  A more  detailed  illustration  of  this  situation,  is  given  by  Fig.  4 where  layer 
densities  of  states  at  the  J*  symmetry  point  are  shown.  Note,  for  example,  how  the  bulk  states 
at  -1 1 eV  on  the  Ge  side  and  at  -12  eV  at  the  GaAs  side  decay  at  the  interface  layers  and  give 
rise  to  a bound  state  at  -10  eV.  The  reason  that  amplitude  at  the  interface  is  shifted  up  in 
energy  and  interface  states  appear  above  the  bulk  bands  is  the  fact  that  Ge-Ga  bonds  are 
weaker  than  either  Ge-Ge  or  Ga-As  bonds,  since  the  strength  of  the  atomic  potentials^  is  in 
the  order  Oa  < Ge  < As.  This  behavior  is  also  consistent  with  the  fact  that  Ge-Ga  bonds  are 
acceptor  bonds,  in  the  sense  that  simple  valence  counting  shows  that  the  two  atoms  participat- 
ing in  such  a bond  contribute  a total  of  1.7S  electrons. 
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In  the  case  of  the  Ge-As  interface  (Figs.  2b  and  3b),  all  the  effects  are  of  the  same 
nature  except  that  the  three  interface  bands  now  appear  below  the  band  edges,  because  Ge-As 
bonds  are  stronger  than  either  Ge-Ge  or  Ga-As  bonds.  Furthermore,  Ge-As  bonds  are  "donor 
bonds",  in  the  sense  that  the  two  atoms  participating  in  Ge-As  bonds  contribute  a total  of 
2.2S  electrons. 

e 

The  results  described  above  for  the  Ge-Ga  interface  are  in  agreement  with  and  comple- 
ment the  results  obtained  by  Baraff,  Appelbaum,  and  Hamann  (BAH)^  who  carried  out  a 
self-consistent  pseudopotential  calculation.  These  authors  used  as  a model  geometry  three  Ge 
overlayers  on  a semi-infinite  GaAs  substrate  and  saturated  the  dangling  bonds  at  the  outermost 
Ge  layer  with  hydrogen  atoms.  Calculations  were  carried  out  at  high  symmetry  points  in  the 
two-dimensional  Brillouin  zone  and  interface  states  were  sought  only  near  the  top  of  the 
projected  valence  bands.  An  interface  band  was  identified  going  through  K and  J'  at  the  points 
shown  as  heavy  dots  in  Fig.  2a  (Because  of  a different  choice  of  coordinate  systems,  the  J 
point  in  Ref.  2 corresponds  to  the  J'  point  in  our  case).  In  addition  a state  at  F at  -0.1  eV 
was  found  which  is  resonant  with  Ge  bulk  states.  This  state  is  also  present  in  our  results,  but 
at  -0.2  eV.  It  behaves  like  a resonance  on  the  Ge  side  and  decays  over  four  layers  on  the 
GaAs  side.  BAH  assumed  that  in  addition  to  the  above-mentioned  interface  band,  which  for 
local  neutrality  must  be  3/4  filled^,  there  exists  anpther  fully  occupied  interface  band  at  lower 
energies  as  is  required  by  local  charge  neutrality.  Our  results  show  that  the  situation  is 
actually  somewhat  more  complicated  in  that  the  interface  charge  is  built  up  from  more  than 
one  interface  band  as  well  as  from  some  bulk  states  (Fig.  3a). 

The  results  for  both  the  Ge-Ga  and  Ge-As  (100)  interfaces  allow  for  an  interesting 
qualitative  comparison  with  the  results  obtained  for  the  (110)  Ge-GaAs  heterojunction  by 
Picket,  Louie  and  Cohen^.  It  should  be  noted  that  the  Ge-GaAs  (110)  interface  contains  both 
types  of  bonds,  Ge-Ga  and  Ge-As  bonds.  One  would  therefore  expect  to  find  interface  bands 
both  above  and  below  band  edges,  very  similar  to  a superposition  of  our  Figs.  2a  and  2b.  This 
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qualitative  conclusion  is  born  out  in  the  results  of  Pickett  et  al.*^ 

Finally,  we  turn  to  a comparison  of  the  theoretical  results  with  available  experiments. 
Esaki  et  alJ^  concluded  from  transport  measurements  that  no  empty  interface  states  (<S  x 
lO'^  cm'^)  are  present  in  the  heterojunction.  Experimentally,  however,  it  is  not  possible  to 
decide  whether  the  sample  contains  a Ge-Ga  interface,  a Ge-As  interface  or  some  other  mixed 
interface.  In  fact,  since  according  to  our  calculations  both  the  ideal  abrupt  interfaces  are 
metallic'^  in  disagreement  with  experiment,  we  believe  that  the  most  jikely  configuration  is  a 
stoichiometricaUy  mixed  interface,  i.e.  there  exists  at  least  one  layer  which  contains  both  Ge 
and  Ga  or  both  Ge  and  As  atoms  so  that  there  are  on  the  average  an  equal  number  of  Ge-Ga 
and  Ge-As  bonds.  In  the  simplest  such  case,  one  can  have  the  last  Ge  layer  followed  by  a layer 
in  which  Ge  and  Ga  atoms  occupy  alternate  sites,  followed  in  turn  by  an  As  layer,  etc.  Such  a 
configuration  would  have  interface  states  very  much  like  those  of  the  (110)  interface^  and  no 
empty  states  would  be  present  in  the  gap  in  agreement  with  experiment. 

In  conclusion,  we  have  presented  a novel  method  to  calculate  the  electronic  structure  of 
two  semi-infinite  semiconductors  in  contact.  We  have  shown  that  the  method  can  provide  a 
detailed  description  of  bound  states,  resonances,  antiresonances  and  local  densities  of  states  for 
semiconductor  interfaces.  The  good  agreement  with  the  results  of  self-consistent  calculations 
provides  further  support  to  the  notion  that  empirical  tight-binding  Hamiltonians  can  provide  a 
meaningful  description  of  the  electronic  structure  of  semiconductors,  their  surfaces  and 
interfaces  in  the  valence  band  region. 

We  would  like  to  thank  S.  G.  Louie,  P.  Maldague  and  A.  R.  Williams  for  valuable 
discussions. 
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For  the  Ge-As  interface,  electron  counting  reveals  that  all  valence  states  and  interface 
states  within  the  valence  band  energy  region  (see  Figs.  2b  and  3b)  are  full  and  an 
excess  of  O.S  electrons  per  two-dimensional  unit  cell  would  have  to  partially  occupy  an 
interface  band  in  the  conduction  bands. 


Fig.  I Interface  creation  by  a local  perturbation.  Each  circle  represents  a two- 
dimensional  plane  of  atoms.  The  perturbation  consists  of  the  "removal"  of  an  As  and  a 
Ce  plane  (arrows)  and  the  appropriate  link  between  a Ge  and  a Ga  plane  as  shown 
(solid  line). 
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Fig.  2 Joint  PBS  and  bound  states  for  the  two  Ge-GaAs  interfaces. 
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Fig.  3 f-incegrated  LOOS  on  the  Ge,  Ga  and  Ge,  As  layers  next  to  the  interfaces  of 
Ge-GaAs  in  comparison  with  the  LOOS  on  the  corresponding  bulk  layers 
(dashed  lines).  The  shaded  areas  show  the  excess  density  due  to  the  interface. 
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Fig.  4 LDOS  on  the  first  five  layers  on  both  sides  of  the  interface  (J*  point). 
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Temperature  und  intensity  dependence  of  rate  of  introduction  of  neutral 
traps  in  Si02  by  e-beam  irradiation* 

J.  A.  Van  Vechten,  C.  M.  Serrano,  and  J.  M.  Aitken 
IBM  Thomas  J.  Watson  Research  Center 
Yorktown  Heights,  New  York  10598 

Recent  experiments  have  shown  that  uncharged  electron  traps  are  created  in  Si02  films 
during  the  course  of  irradiation  with  x*rays*  and  with  electron-beams^  at  energies  suffi- 
ciently low  that  Frenkel  pair  formation  is  excluded.  These  defects  are  in  addition  to  the 
well  known  positively  charged  electron  traps  which  are  also  introduced  by  irradiation.  The 
neutral  defects  are  not  removed  by  the  usual  ~ 400  C anneals^  which  remove  the  positively 
charged  defects.  Because  any  electron  traps  that  remain  in  the  Si02  of  the  final  device  will 
aggravate  hot-electron  related  threshold  instabilities,  and  because  many  devices  are  exposed 
to  ionizing  radiation  at  some  point  in  their  fabrication  (during  e-beam  or  x-ray  lithography, 
during  reactive  ion  etching,  RIE,  etc.),  it  is  important  to  understand  the  nature  and  origin 
of  the  trapping  centers. 

Here  we  report  the  dependence  of  the  radiation  induced  neutral  trap  density  upon 
e-beam  intensity  and  upon  sample  temperature  during  irradiation.  A model  is  proposed  to 
account  for  these  results 

The  experiments  were  performed  on  50  nm  Si02  films  grown  on  0.2  Q cm,  p-type 
<100>  substrates  at  1000  C in  dry  O2.  500  nm  thick  A1  dots  were  rf  evaporated  through 
a shadow  mask  to  form  capacitor  structures  which  were  subsequently  annealed  at  400  C for 
20  min.  in  forming  gas.  Individual  capacitors  were  diced  from  the  wafers,  placed  in  the 
column  of  a JEOL  200  scanning  transmission  electron  microscope,  STEM,  on  a liquid  N2 
cold  stage.  Samples  were  irradiated  with  6 nA  of  50  KV  electrons  at  both  room  temper- 
ature and  at  the  temperature  of  the  cold  stage  when  filled  with  liquid  N2,  which  is  about  90 
K.  (The  same  arrangement  was  used  at  both  temperatures  to  insure  that  the  exposure 
conditions  were  otherwise  the  same.)  Exposure  intensity  was  controlled  by  varying  the 
diameter  of  the  beam.  The  entire  area  of  each  sample  was  scanned  and  the  total  time  the 
beam  scanned  the  sample  was  the  same  for  all  intensities.  Thus,  the  total  dose  given  all 
samples  was  the  same  and  all  but  the  highest  intensity  samples  were  multiply  scanned. 

After  irradiation,  the  samples  were  annealed  in  forming  gas  at  400  C for  30  min.  This 
removed  all  the  positively  charged  and  some  of  the  neutral  centers  that  were  introduced  by 
the  radiation.  The  remaining  neutral  traps  were  studied  by  injecting  hot-electrons  from  the 
substrate  into  the  oxide  where  a small  fraction  were  trapped  at  the  defects.  The 
resulting  shift  in  flat-band  voltage  is  proportional  to  the  number  of  electrons  trapped  in  the 
oxide.  Fig.  1 shows  total  flat-band  shift  for  samples  irradiated  at  1x10^  and  at  0.1  W/cm^ 
for  5 min.  plus  an  unirradiated  control  sample  plotted  against  the  time  that  a constant  d.c. 
current  of  2 x 10*^  A is  passed  through  the  oxide.  We  see  that  the  flat-band  shift  in  the 
samples  irradiated  at  the  higher  intensity  is  less  for  low  temperature  than  for  room  tempera- 
ture irradiation,  that  the  shift  in  samples  irradiated  at  room  temperature  is  less  for  low 
intensity  irradiation  although  the  total  dose  is  the  same,  and  that  for  low  temperature,  low 
intensity  irradiation,  the  shift  is  essentially  the  same  as  that  in  the  unirradiated  control. 
Analysis  of  the  trapping  cross  sections  of  the  trapping  centers  showed  that  the  were  two 
types  of  defect  present  in  the  oxide.  One  defect  had  a cross  section  of  about  10*'^  cm^  and 
was  not  present  in  samples  irradiated  at  90  K.  The  other  defect  had  a cross  section  of 
about  lO*'^  cm‘  and  was  not  present  in  samples  irradiated  at  intensities  below  10  W/cm‘. 
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We  propose  the  following  explanation  of  these  results.  From  the  fact  that  the  neutral 
defects  in  question  do  not  anneal  out  we  conclude  that  they  are  more  stable  configurations 
of  the  material  than  that  in  the  as  grown  oxide.  The  energy  of  the  incident  radiation  is 
much  lower  than  the  threshold  required  for  direct  lattice  displacements.  However,  point 
defects  (vacancies  and  impurities)  may  migrate  by  means  of  recombination  enhanced 
diffusion^,  RED,  during  irradiation  to  form  bound  complexes.  These  processes  are  normally 
thermally  activated  but  depend  only  on  dosage  and  not  on  intensity.  An  intensity  depend- 
ent mechanisim  is  required  to  explain  the  data  in  Fig.  1.  Dislocations  are  known  to  exist  in 
amorphous  materials^.  Misfit  dislocations  have  been  induced  to  glide  in  GaAlAs  heteros- 
tructures by  ionizing  radiation^.  This  process  has  a sharp  threshold^  at  about  10^  W/cm^. 
Such  motion  is  reasonable  in  the  present  case  since  migration  to  the  interface  would  remove 
the  misfit  strain  due  to  the  difference  in  thermal  expansion  coefficients  between  Si  and 
Si02.  Both  types  of  process  are  required  to  explain  the  data  in  Fig.  1.  We  have  no 
explanation  for  the  apparent  discrepancy  with  the  results  of  EerNisse  and  Norris^  who 
conclude  interfacial  strain  is  increased  by  irradiation. 

*This  work  was  supported  in  part  by  the  Office  of  Naval  Research  under  contract  No. 
N00014-76-C-0934. 
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ABSTRACT 

We  solve  the  effective  mass  equations  describing  donor  impurities  localized  at 
semiconductor-insulator  interfaces  using  a method  which  involves  an  expansion  of  the  impurity 
wavefunction  in  terms  of  spherical  harmonics.  We  also  investigate  the  behaviour  of  the  | 

impurity  energies  as  a function  of  the  donor  distance  from  the  interface  both  in  the  insulator 
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and  in  the  semiconductor.  This  allows,  for  the  first  time,  the  connection  between  surface  and  i 
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bulk  donor  properties.  Results  for  the  Si-Si02  case  are  presented.  i 
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I.  Introduction 


A number  of  properties  of  MOS  devices  are  strongly  influenced  by  the  presence  of 
impurities  at  the  semiconductor-oxide  interfaces.  Karpushin'  and  Bell  et  al~  have  investigated 
the  behavior  of  donor  localized  at  semiconductor  surfaces  without  inversion  layers.  The  effect 
of  inversion  layers  was  studied  by  Stern  and  Howard^,  but  their  numerical  results  are  restricted 
to  the  two-dimensional  limit.  Very  recently*.  Martin  and  Wallis  have  given  a theoretical 
treatment  of  wavefunctions  and  binding  energies  of  donors  localized  at  semiconductor  surfaces 
in  inversion  layers,  using  simple  variational  wavefunctions  which  are  adequate  only  for 
extremely  high  electric  fields.  They  show  that  the  two-dimensional  limit  is  approached  only 
for  unrealistically  high  electric  fields. 

The  purpose  of  the  present  paper  is  twofold.  First,  to  solve  more  accurately  the  same 
hamiltonian  as  that  used  by  Martin  and  Wallis  for  all  realistic  values  of  the  electric  field  in  the 
inversion  layer.  Second,  to  extend  the  treatment  to  include  the  case  of  donors  situated  not  just 
at  the  semiconductor-oxide  interface  but  also  near  it.  i.e.  as  a function  of  the  donor  distance 
from  the  interface  both  in  the  insulator  and  in  the  semiconductor. 

In  section  11.  we  present  the  equation  for  the  case  of  a donor  impurity  situated  at  the 
interface  and  define  the  basic  approximations  using  in  deriving  it.  The  method  of  solution  is 
also  presented.  In  section  III  the  case  of  donors  situated  near,  and  not  just  at.  the  interface  is 
discussed  and  solved.  Numerical  results  are  presented  in  section  IV  and  discussed  with 
particular  reference  to  the  recent  experimental  data  of  Hartstein  and  Fowler^. 

II.  Theoretical  framework  and  method  of  solution. 

We  consider  two  contiguous,  semi-infinite  half  spaces  and  assume  that  the  potential 
energy  of  an  electron  undergoes  a discontinuous  jump  as  the  electron  passes  from  the  semicon- 
ductor to  the  insulator.  For  the  St-Si02  case  this  discontinuity  is  about  3 eV^.  Here  we  shall 


assume  that  this  discontinuity  is  infinite  and  that  the  electron  wavefunction  vanishes  at  the 
boundary.  The  effect  of  this  approximation  on  the  binding  energy  has  been  shown  to  be 
negligible^.  In  the  region  of  the  inversion  layer  there  is  an  electric  field  present  which  is  in 
general  a complicated  function  of  the  distance  from  the  interface.  We  shall  assume  the  electric 
field  to  be  constant,  which  has  been  shown  to  be  a good  approximation^.  For  simplicity,  we 
assume  the  <I00>  orientation  for  Si  (for  Ge  it  would  be  the  <11 1>)  even  though  the  present 
method  could  be  applied  to  any  other  orientation.  For  this  situation,  Karpushin*  has  shown 
that  the  lowest  lying  bound  states  are  primarily  derived  from  the  energy  ellipsoids  whose  major 
axis  are  perpendicular  to  the  interface.  Therefore  only  these  ellipsoids  will  be  considered  here. 
We  also  neglect  intervalley  interaction  between  the  two  ellipsoids  perpendicular  to  the  surface. 
With  the  above  assumptions  the  Hamiltonian  for  an  impurity  of  charge  e can  be  written  as 
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where  E is  the  external  electric  field,  m,  and  m^  are  the  transverse  and  longitudinal  effective 
masses,  respectively;  cj  and  e2  are  the  dielectric  constants  of  the  insulator  and  semiconductor 
respectively,  as  illustrated  in  figure  la.  Equation  (I)  has  been  written  using  as  unit  of  energy 
and  length,  respectively,  the  following  quantities 


(3a) 

(3b) 

Finally,  the  potential  energy  U (7)  is  approximated  by 

U (7)  « - 1 (4) 

i.e.  the  screening  of  the  impurity  ion  potential  by  the  electrons  in  the  inversion  layer  is 
neglected.  Hamiltonian  ( 1 ) is  exactly  the  same  hamiltonian  as  that  considered  by  Martin  and 
Wallis-». 

To  solve  equation  (1).  we  note  that,  in  the  absence  of  the  electric  field,  the  image 
potential  and  the  anisotropy  (i.e.  y a 1),  Levine^  has  shown  that  a state  cannot  exist  unless 
f-fm  is  an  odd  integer,  with  the  ground  state  being  the  2p„  level,  in  the  usual  hydrogen  atom 
nomenclature.  When  the  anisotropy,  the  image  potential  and  the  electric  field  are  switched  on 
I is  no  longer  a good  quantum  number  but  m still  is.  Therefore,  one  can  write  the  solution  of 
I ( I ) as 

(7)  - 2 fi  (r)  (e.g,)  (5) 

i 

where  the  angles  0£9s«r/2  and  0<<f<2v.  As  a result  of  the  angle  B varing  only  in  half  space, 
the  spherical  harmonics  obey  different  orthogonality  rules^.  In  particular,  spherical  harmonics 
I with  different  m are  still  orthogonal,  but  those  with  different  f are  orthogonal  only  if  is 

even.  The  matrix  elements  of  all  operators  appearing  in  (1)  in  the  basic  function  (5)  can  be 
calculated  analytically,  and  they  will  be  given  elsewhere*.  As  a result,  one  is  left  with  a system 
of  radial  differential  equations,  whose  order  is  given  by  the  number  of  terms  in  expression  (5). 
This  system  can  be  solved  by  using  the  following  expansion 


2m,  e"^ 
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This  allows  one  to  transform  the  problem  from  that  of  solving  a system  of  differential  equa- 
tions to  that  of  diagonalizing  a secular  matrix.  The  size  of  the  secular  determinant  is  given  by 
the  product  of  the  number  of  terms  (N)  in  expansion  (S)  times  the  number  of  terms  (M)  in 
(6).  The  convergence  of  the  results  with  respect  to  the  number  of  these  terms  and  the  choice 
of  the  parameters  aj  is  discussed  in  section  IV. 

In  order  to  calculate  binding  energies'*,  we  must  also  solve  Hamiltonian  ( 1 ) in  the  absence 
of  an  impurity  potential.  As  shown  by  Martin  and  Wallis'*,  it  is  sufficient  to  solve  the  Hamilto- 
nian. 


Ho  - - Y + - + aEz 
dZ“  z 


(7) 


This  one-dimensional  Hamiltonian  can  be  solved  very  accurately  by  using  an  expansion  of 
the  type  (6).  i.e..  by  assuming 


4'(z) 
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and  diagonalizing  the  corresponding  secular  matrix,  as  discussed  in  section  IV.  The  binding 
energies  are  therefore  given  by 


Eh  - E(Ho)  - E(H) 


(9) 


III.  Impurity  located  near  the  interface. 


In  the  previous  section  we  have  considered  the  situation  in  which  the  impurity  is  located 
exactly  at  the  interface.  We  now  consider  the  situation  in  which  the  impurity  is  located  either 
in  the  semiconductor  or  the  oxide  (fig.  Ib).  In  the  first  case,  the  potential  experienced  by  the 
electron  is  given  by 
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where 


r,  - v' x“+y^+(d-z)^ 
r2  - ^ x-+y-+(d+z)“ 


(10) 
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d being  the  distance  of  the  donor  ion  from  the  interface.  In  the  second  case.  i.e.  when  the 
donor  is  located  in  the  oxide. 


V,  (r)  - - i- 
'’2 
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In  order  to  calculate  the  matrix  elements  of  (10)  and  (12)  in  the  basis  (S)  we  use  the  well 
known  expression 
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Pf  (cos  y) 


(13) 


where  r^  (r^)  is  the  smaller  (large)  of  |T|  and  |^|.  and  y is  the  angle  between  T and 


Since  d is  chosen  along  the  z axis,  we  can  write 


r 


The  matrix  elements  of  (14)  and  (IS)  in  the  basis  states  (S)  can  be  also  calculated 

I analytically^.  In  actual  calculations  the  number  of  terms  in  (14)  and  (IS)  is  kept  finite  and  the 
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convergence  of  the  results  with  respect  to  the  number  of  terms  is  discussed  in  the  next  section. 

IV.  Numerical  results  and  discussions 

All  the  calculations  reported  here  are  for  the  case  of  the  Si-SiOt  interface  using  the 

* j 

parameters  given  by  Stern  and  Howard^.  We  first  address  the  problem  of  convergence  of  the  \ 

calculated  binding  energies  as  a function  of  the  basis  functions  in  (6).  All  the  results  reported 
here  are  for  the  ground  state  but  the  analysis  can  be  extended  in  a straightforward  manner  to 
treat  other  excited  states  as  well.  The  exponents  aj  should  in  principle  be  treated  as  variation- 
al parameters.  In  practice,  however,  as  we  have  previously  seen^,  the  aj  can  be  fixed  and  not 
considered  as  variational  parameters  if  one  chooses  a sufficiently  large  number  of  terms  in  (6) 
and  if  the  a/s  are  selected  in  a sensible  way.  We  have  seen  that  it  suffices  1)  to  include  IS 
terms  in  (6),  and  2)  to  choose  the  largest  and  the  smallest  aj  to  be  50  and  O.OS  respectively, 
with  the  rest  of  them  in  geometrical  progression,  to  insure  a convergence  in  the  eigenvalues  to 
better  than  1 part  in  10**.  The  next  question  regards  how  many  terms  should  be  included  in 
(S).  This  number  obviously  should  increase  with  the  strength  of  the  electric  field.  In  table  I 
we  show  the  contribution  of  the  various  terms  in  expression  (S)  for  fields  up  to  10**  esu,  which 
corresponds  to  3x10'*  volt/cm.  Since  10^  esu  corresponds  to  the  physically  highest  achievable 
fields,  we  see  that  usually  less  than  S terms  provide  very  accurate  solutions.  Also  Hamiltonian 
(7)  can  be  solved  very  accurately  using  IS  terms  in  expansion  (8).  In  fig.  2 we  compare  our 


results  with  the  previous  results  of  Martin  and  Wallis  and  we  see  that  our  analysis  is  more 
accurate  and  that  the  two  methods  agrees  for  very  large  values  of  the  electric  field.  In  figure  3 
we  illustrate  the  individual  effect  of  the  anisotropy,  and  of  the  image  potential  on  the  binding 
energy.  We  see  that,  for  small  and  intermediate  fields,  all  these  effects  are  important,  whereas, 
as  shown  by  Martin  and  Wallis,  the  binding  energy  tends  always  to  4 in  the  limit  of  infinitely 
large  electric  fields,  which  corresponds  to  the  two-dimensional  limit.  In  fig.  4 we  plot  the 
mean  radius  for  the  ground  state  as  a function  of  the  electric  field. 

t 

We  now  discuss  the  case  in  which  the  donor  ion  is  not  located  at  the  interface  but  near  it. 
In  fig.  5 we  plot  the  binding  energy  of  the  ground  state  as  a function  of  the  distance  of  the 
impurity  ion  from  the  interface  in  the  absence  of  the  electric  field.  In  the  case  in  which  the 
donor  moves  in  the  semiconductor,  the  binding  energy  increases  and  reaches  a value  which 
corresponds  to  the  bulk  donor  binding  energy  for  the  ground  state.  It  is  important  to  note  that 
the  value  of  ~.7  is  in  units  of  the  Rydberg  (3a)  appropriate  for  the  interface  while  for  the  bulk 
one  would  define  a Rydberg  in  terms  of  63  only;  this  means  that,  in  terms  of  the  bulk  effective 
rydberg,  one  has  for  the  binding  energy  the  value  ~1.6  which  is  in  exact  correspondence  to  the 
value  calculated  using  bulk  analysis*^.  From  fig.  Sa  we  see  also  that  the  bulk  value  is  reached 
after  2 effective  radii  (i.e.  after  about  40  A),  and  that  for  distances  from  the  interface  larger 
than  that  the  binding  energy  does  not  change  since  the  radius  for  the  ground  state  is  smaller 
than  the  distance  of  the  donor  ion  from  the  interface.  In  fig.  5b  we  show  the  dependence  of 
the  binding  energy  of  the  ground  state  as  a function  of  the  impurity  distance  from  the  interface 
in  the  oxide.  We  see  that  the  binding  energy  decreases  as  the  distance  increases,  as  expected. 

We  now  briefly  compare  our  theoretical  results  with  the  experimental  findings  of  Hart- 
stein  and  Fowler^,  who  have  measured  the  conductivity  of  n-channel  silicon  MOSFET  devices 
in  which  Na*^  ions  were  diffused  through  the  oxide  to  the  semiconductor-oxide  interface.  The 
temperature  dependence  of  the  peak  conductivity  can  be  described  by  an  activation  energy 
which  Hartstein  and  Fowler  identify  as  the  binding  energy  of  an  impurity  band.  By  varying 


the  substrate  bias  Hartstein  and  Fowler  were  able  to  study  the  activation  energy  as  a function 
of  the  electric  field.  For  electric  fields  of  19.9  esu  and  60.4  esu  they  found  activation  energies 
of  18  meV  and  25  meV  respectively.  In  fig.  6 we  show  the  calculated  binding  energy  as  a 
function  of  the  electric  field  for  the  cases  in  which  the  impurity  ion  is  at  the  interface  and  at  .2 
a.u  (~4  A)  in  the  oxide.  We  see  that  the  observed  binding  energies  agree  much  better  with 
the  case  in  which  the  impurity  ion  is  located  not  exactly  at  the  interface  but  very  close  to  it.  It 
is,  however,  important  to  note  that  the  present  model  neglects  important  effects  as  discussed  in 

t 

the  beginning  of  section  II.  Therefore  one  should  take  the  result,  presented  in  fig.  6.  as  only  a 
qualitative  indication  rather  than  quantitative  determination  of  distance  of  the  impurity  ion 
from  the  interface.  Fig.  6 shows  that  the  solution  is  sensitive  to  the  distance  of  the  donor  ion 
from  the  interface. 

In  conclusion,  we  have  presented  a method  which  provides  accurate  binding  energies  of 
donors  located  at  semiconductor-insulator  interfaces  for  any  value  of  the  anisotropy  and  of  the 
electric  field.  We  found  that  the  binding  energy  depends  strongly  on  the  distance  of  the 
impurity  from  the  interface.  Clearly,  one  should  now  improve  the  theoretical  model  before  an 
accurate  description  of  the  experimental  data  is  attempted.  In  particular,  the  screening  of  the 
impurity  ion  potential  due  lo  the  electrons  in  the  inversion  layer  has  to  be  included.  Further- 
more. intervalley  coupling  has  to  be  considered  and  a better  model  for  the  interface  should  be 
used.  This  aspect  is  currently  under  investigation. 
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Table  I 


Contribution  of  various  terms  in  the  expansion  (S)  given  in  text  for  different  values  of  the 
electric  field  E.  Only  the  I * odd  terms  contribute  tot  he  ground  state.  Contributions  are 
given  in  percent. 


0 91.2  8.3  .5 

1 91.1  8.4  .5 

10  90.3  9.1  .6 


102  85.2  13.1  1.5  .2 

10^  70.7  21.9  5.5  1.4  .4  .1 

lO-*  50.7  20.5  12.5  5.2  2.1  .8  .3  .1 


Figure  Captions 


Fig.  1.  Schematic  structure  of  semiconductor  insulator  structure.  In  fig.  (la)  the 

impurity  ion  is  located  at  the  surface  z— 0.  In  fig.  (lb)  the  impurity  ion  is 
located  in  the  semiconductor  (1)  or  in  the  oxide  (2).  at  a distance  d from  the 
interface.  The  coordinate  axes  are  also  shown. 

Fig.  2.  Binding  energy  for  the  ground  slate  as  a function  of  the  electric  field  for  the 

Si-SiOi  case.  Comparison  with  the  results  of  Martin  and  Wallis,  ref.  4 in  text. 

I 

is  shown. 

Fig.  3.  Impurity  ground  state  binding  energy  as  a function  of  the  electric  field  for  the 

Si-SiO->  case.  The  effects  of  the  image  potential  and  of  the  conduction  band 
anisotropy  are  illustrated. 

Fig.  4.  Effective  radius  for  the  impurity  ground  state  as  a function  of  the  electric  field. 

The  symbols  are  the  same  as  in  fig.  3. 

Fig.  5.  Impurity  ground  state  binding  energy  as  a function  of  the  impurity  ion  from 

the  interface.  Fig.  (5a)  show,  the  case  in  which  the  donor  ion  is  located  in  the 
semiconductor,  while  fig.  (5b)  represents  the  case  in  which  the  donor  ion  is 
located  in  the  insulator.  The  numerical  values  correspond  to  the  Si-SiO^ 
situation. 

Fig.  6.  Impurity  ground  state  binding  energy  as  a function  of  the  electric  field  for  two 

values  of  the  distance  of  the  impurity  ion  from  the  interface.  The  first  value 
(daO)  corresponds  to  the  impurity  ion  located  at  the  interface;  the  value 
daO.2  eff.  units  corresponds  to  the  impurity  ion  located  in  the  oxide. 
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